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1. Introduction 

Jupiter's bands are a clear indication of atmospheric 
differences at different lattitudes on the planet. The aim of 
the present study was to attempt infrared observations at dif- 
ferent lattitudes in order to obtain spectra xn the 10y_ region. 
These spectra might then help one to understand differences 
in chemical composition or physical structure of the well known 
optical features. 

In order to obtain such spectra of a rotating planet, it 
seemed well to attempt simultaneous observations at different 
lattitudes. We planned to use a Hadamard transform spectrometer 
with 15 entrance slits in order 'to obtain 15 simultaneous spectra, 
at a resolution of O.Oly. The spectral band covered contained 
255 spectral elements. 

2 . Results 

The results obtained during the course of the study are 
covered in ten publications which include seven ;]ournal articles, 
two general review papers, and a Ph.D. thesis. These publications 
are listed in section 3 below. The thesis which is the best 
single summary of the work done under the present contract appears 
as the Appendix to this report. 

The publications fall into two main groups; those that deal 
with technical aspects of the instrumentation we had to develop 
in order to make observations; and those that deal with obser- 


vational work. 
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The technical difficulties encountered turned out to be 
rather more severe than had been initially anticipated. Many 
of the problems we had to solve are of a very general nature that 
applies to other varieties of multiplex instrumentation, as well 
as to other classes of Hadamard transform spectrometers. As a 
result, all the technically innovative work was of sufficiently 
general interest to warrant publication. 

The astronomical observations that were undertaken ,are 
described in the thesis of Ming-Hing Tai, the graduate student 
who worked on the project (see Appendix) . The results cited 
there represent the first observations undertaken with our ins- 
trument. The instrument did in fact perform with the spectral 
resolving power and the number of spatial and spectral elements 
we had anticipated in our design. Because of the above mentioned 
technical difficulties, we suffered , considerable delays in 
in initiating observations. Ultimately we therefore had to content 
ourselves with far less observing time than we had originally 
hoped to obtain, and the observational side of the program had 
to be carried out at far more modest levels than anticipated 
because important technical requirements had drained a dispro- 
portionately large fraction of our available funds. 

Despite these setbacks, however, we feel that the results 
obtained — while not precisely in line with what we had set out 
to do — are interesting, led to a large number of publications, 
and constitute a worthwhile contribution to observational planetary 


science. 
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3. Publications 

1. Two asymmetric Hadamard transform spectrometers, Martin 

Harwxt, Perry G. Phillips, Leon W. King, and Daniel A. Briotta, Jr 
Applied Optics, 2669 (1974) • 

2. Infrared astronomy and the shuttle, Martin Harwit, Astro- 
nautics and Aeronautics 1 ^, 23 October (1974) 1. 

3. Masks for Hadamard transform optics, and weighing designs, 

Neil J. A. Sloane and Martin Harwit, Applied Optics 207 (1976) 

4. A practical multi-spectrum Hadamard transform spectrometer, 

Ming Hing Tai, D. A. Briotta, Jr., N. Kamath and Martin Harwit, 
Applied Optics 2533 (1975) . 

5. Errors in Hadamard spectroscopy or imaging caused by imperfect 
masks, Ming Hing Tai, Martin Harwit and Meil J. A. Sloane, 

Applied Optics 3^, 2678 (1975) . 

6. Fourier and Hadamard transform spectrometers: a limited 

comparison, Ming Hing Tai and Martin Harwit, Applied-Optics 15 , 
2664 (1976) . 

7. "Hadamard Transform Analytical Systems", Martin Harwit, 
to be published in "Transformations in Analytical Chemistry," 

P. Griffiths, Ed., Plenum Press (1978). 

8. Fourier and Hadamard transform spectrometers: a limited 

comparison II. M. Harwit and Ming Hing Tai, Applied Optics 

3070 (1977) . 

9. Distortion in Hadamard Transform Optics, N.J.A. Sloane, 

Martin Harwit and Ming-Hing Tai, accepted for publication. 

Applied Optics (1978) . 



4 


10. A practxcal Hadanard transform spectrometer for astronomical 
application, Ming-Hing Tai, Cornell University Ph.D. thesis and 
Center for Radiophysics and Space Research Report CRSR 655 (1977) . 

4 . Personnel Participating on the Grant 

Prof. M. Harwit, principal investigator 
Mr. Ming Ring Tai, graduate student 
Mr. Gary Melnick, graduate student 

Mr. Frank W. Dain, undergraduate student (part time) 

Mr. Gerald Stasavage, research technician (part time) 

Mr. Daniel A. Briotta, Jr. , graduate student (part time) 

Mr. Leon W. King, graduate student (part time) 

Mr. Narayana Kameth, programmer (part time) 

5, Financial Status 

Al-1 funds provided under this contract have been spent 


by this time. 



CENTER FOR RADIOPHYSICS AND SPACE RESEARCH 


CORNELL UNIVERSITY 
ITHACA, NEW YORK 


January, 1977 


A PRACTICAL HADAMARD TRANSFORM SPECTROMETER 

FOR 

ASTRONOMICAL APPLICATION 


by 

Ming-Hing Tal 



DEDICATION 


To my parents 




ill 


ACKNOWLEDGEMENTS 


This project is by no means one man’s product. Many 
people have devoted their talents to the completion of this 
project. I am just the one to make a summary of it. 

My special thanks go to Prof. Martin Harwitj my thesis 
advisor. Martin set himself up as an admirable^ example: he 
always moved the heaviest equipment himself; tackled the most 
difficult jobs in the experiment so that he himself was res- 
ponsible for any misfortunes; was earnest in discussions 
about problems with his student; and was modest in sharing 
the credit for the results. Above all, he truly expressed 
concern about his student, as a respectable advisor and as 
a warm friend. 

My second thanks go to Dr. Daniel A. Briotta_._ Dan helped 
to automlze the whole process and taught me the mini computer 
technique. Dan has shown me how to run the computer as though 
he were showing me how the lines in his palm run. 

My thanks also go to the infrared astronomy group at 
Cornell. Prof. J. Houck has served on my committee and has 
always been available for discussion of astronomical and ex- 
perimental problems. I especially thank him for his help in 
finishing up my degree while Martin was on sabbatical leave. 
Dr. Phillips built the spectrometer. The late L. King modi- 
fied it. G. Stasavage built some electronic parts and H. 
Kondracki built the dewer. G. Melick helped with the obser- 


Iv 



vatlon, data processing and electronic parts, G, Gull with 
the field trip, and Drs. W. Forrest, Diward, and D. Shacck 
were involved in many valuable discussions. Thanks also go 
to Westy Dain for not only being always willing to help on 
chilly night's observations and tedious laboratory work, 
but also for proof reading my thesis. 

I thank Prof. E.E. Salpeter for serving on my committee. 
Prof. A. Siever and Allen Chin for helping the He^ cooled 
dewer, and Prof. P. Glerasch for always being available for 
discussing astronomical problems. I also thank G. Vincent 
for building many mechanical parts, B. Boettcher for figures 
and Eckelmaun for photographs. 

The author thanks the staffs of the Kltt Peak National 
Observatory and McDonald Observatory for generous allotments 


of observing time and for;_ as si S;tanc^ during- the- -ob serving-runs ; 
Finally, my thanks go to Divan Sun, Adeline Wu, and W.¥. 


Chiu, for typing the thesis for me, and to N.H. Cheung, 
K.K. Mon, K.M. Lo and many of ray Chinese friends for their 


continual help during these years. 


Portions of this work were supported by the NGR 33-010-210 


and NSG 1263- 


V 



TABLE OF CONTENTS 


PAGE 

LIST OP FIGURES vii 

LIST OP TABLES xl 

ABSTRACT Xii 

I. HISTORICAL INTRODUCTION TO HADAMARD TRANSFORM 1 

SPECTROMETRY 

II. HADAMARD MATRICES 10 

(A) . Weighing Designs 10 

(B) . General Mathematical Formulation 13 

(C) . Hadamard Matrix 16 

(D) . Optical Realization of Hadamard Encoding 20 

(E) . Errors in Hadamard Spectroscopy 34 

III. INSTRUMENTATION 37 

(A) , Multislit Spectrometer 37 

(B) . The Optics 42 

(C) . The Electronics , 5^ 

(D) . Laboratory Calibration 60 

IV. COMPARISONS BETWEEN FOURIER TRANSFORM AND 65 

HADAMARD TRANSFORM SPECTROSCOPY 

(A) . Mathematically 65 

(B) . Computer Requirements 67 

(C) . Optics 68 

(D) . Mechanical Requirements 73 

vi 



V. PROGRAMMING FOR HADAMARD TRANSFORM SPECTRAL 75 

DATA REDUCTION 

(A). HTS Program ' 75 

'(&).' DHTS Program 8l 

(C). Correction Program 89 

VT. ASTRONOMICAL OBSERVATION 90 

(A) . Correction Procedure 90 

(B) . Observation of a-Orionis 91 

(C) . Observation of Jupiter 100 

(D) . Observation of Mercury 121 

APPENDIX A. ESTIMATE OF CODING ERROR FOR FOURIER 136 

TRANSFORM SPECTROMETRY 

APPENDIX B. SINGLY HADAMARD TRANSFORM' PROGRAM 139 

APPENDIX C. DOUBLY HADAMARD TRANSFORM PROGRAM l68 

APPENDIX D. CORRECTION PROGRAM 201 

•“BI-BLIOGRAPHY- “"2^ 


vil 



LIST OF FIGURES 


FIGURE 

1-1 

2-la 

2-lb 

2-lc 

2-2 

2-3 

2-4 

2-5 

2-6 

2-7 


TITLE 

An 8x8 Hadamard matrix and two 7x7 cyclic 
matrices that can be derived from it, 

Hadamard transformation of a Sine wave 
input . 

Hadamard transformation of a square wave 
input . 

Hadamard transformation of a straight 
line input. 

Schematic representation of HTS using 
the S code. 

The 509 exit slit mask and the 255- 
elemenu ^ code . 

Comparison spectrum of the mercury emi- 
ssion lines in the 1.4-1.8iim region: (a, 
bottom) as obtained in the Hadamard- 
transform mode^ (b, top) as obtained 
under identical conditions using the 
same optical system as a conventional 
monochromator . This figure is taken 
from Decker, 1971b. 

Schematic representation of DHTS. 

The 29 entrance slit mask and the 15- 
element S_ code. 

(a) Spectrum of the 1.7m mercury vapor 
doublet showing negative peaks to the left 
at the emission peaks; (b) Shows the res- 
ponse we would obtain to a single spectral 
line with a perfect mask; (c) Shows the 
response for a single line with the ra- 
diation simulated as passing through a 
mask with slits too narrow because each 
opaque mask element protrudes into the 
adjacent transparent slot by a tenth of a 
slot width: (d) Shows the effect of simu- 
lating slits that are systematically too 
wide. Note that the main spectral line 
has been placed in different positions 
for the synthetic runs (b),(c), and (d) . 


PAGE 

3 

21 

22 

23 

25 

26 

28 

29 

33 

36 


vlii 



FIGURE 

3-1 

3-2 

3-3 

3-il 

3-5 

3-6 

3-7 

3-8 

3-9 

3-10 


5-1 

5-2 

5-3 


TITLE 

The flow chart of the data taking pro- 
cess, starting with radiation from the 
telescope and ending with the output of 
the computers' 

Optical path through the spectrometer. 

The dedisperser and exit mask are shown 
rotated by 90°. 

Grating diffraction efficiency as a func- 
tion of slit position at 8 and l4vim. 

Liquid helium cooled post optics. 

HTS alignment circuit . 

Logic diagram for mask motion and data 
processing. 

HTS drive imlt circuit. 

Basic detector bias and preamplifier 
circuit . 

Spectrum of the mercury vapor 1.7ym 
line using the 1x255 mode. 

Calibration spectra obtained for the 
- -me-r c ury— v-ap o r— e mi-s s 1 on— 1 i nes— at— l-.-6-9-Vim — 
and l.Tlym: (a) The eighth of a series 
of fifteen individual spectra obtained; 

(b) An average of all fifteen spectra; 

(c) A diaplay of the fifteen spectra, 
each spectrum being displaced vertically 
from, the next. The diagonal pattern 
near the right-hand edge represents a 
displacement of the peak between suc- 
cessive spectra. This represents the 
actual shift in spectral range between 
adjacent spatial elements. 

The flow chart of the 1x255 Inverse 
transform program. 

The flow chart of the 15x255 inverse 
- transform-program. 

Matrices generated by the computer 
during the reduction of the spectral 
data. 


PAGE 

kl 

^3 

45 

47 

55 

56 

58 

59 

62 

64 


79 

83 

86 


lx 



FIGURE 


TITLE 


PAGE 


\ 

6-1 The spectral emissivity of various re- 92 

gions as calculated from the flight 
data (Murcray ^ 1970). 

6-2 The raw spectra of a-orionis and the 9 4 

Moon. 

6-3 The ratio spectrum of a-orionis to 95 

the Moon corrected for lunar tem- 
perature. 

6-'^ (a) Low resolution spectrum taken by 99 

Gillett et al (1969). Different sym- 
bols represent spectra taken on dif- 
ferent nights, (b) High resolution 
soectrum uaken by Treffers and Cohen 
(1973). 

6-5 The raw spectrum of Jupiter and the 101 

Moon. 

6-6 The ratio spectrum of Jupiter to the 102 

Moon corrected for lunar temperature. 

6-7 The atmospheric profile of Jupiter for 103 

a solar-composition model. 

6-8 (a) Schematic representation of the 108 

molecule. _ The components of an- 
gular momentum and the motion in the 
V 2 vibrational mode are also shown. 

(b) Energy levels of the V 2 vibrational 
mode of ammonia. Superscripts a and s 
refer to the antisymmetric and symme- 
tric levels which arise due to inver- 
sion splitting. 

6-9 Low resolution spectrum taken by Gillett 109 

( 1969 ). Different symbols repre- 
sent spectrum taken on different nights. 

This figure is taken from their paper. 

5-10 (a) Room temperature absorption spec- 111 

trum of ammonia, p= 0.06 atmos, w=0.6 cm 
atmos. (b) Brightness temperature; 

(c) Surface brightness of the central 
region of Jupiter from 8 to 13.5vm. 

(Taken from Aitken and Jones). 



FIGURE 


TITLE 


PAGE 


6-lla Spectrum of the N and S polar regions 113 

of Jupiter at 3—4 cm**^ resolution di- 
vided by the spectrum of the Moon. 

Data points are shovm as solid circles, 
and the solid line represents the best 
fitting synthetic spectrum calculated 
from the model calculated by Lacy ^ 
al (1975) (The graph is taken from Lacy 
et al ) . The dotted curve is our ob- 
served spectrum' by matching Lacy’s 
spectrum at points A and B, 

6-llb Same as 6-lla except matching our spec- ll4 

trum with Lacy’s at A' and B’. 

6-12 (a) Thermal emission spectrum of Jupiter 117 

corrected for absorption in the earth’s 
atmosphere observed by Ridgway (1973). 

The dashed line is the predicted form of 
the H 2 continuum, (b) The ratio of the 

Jovian spectrum to the atmospheric ab- 
sorption spectrum observed by Combes et 
al ( 1974 ). The solid and dashed lines 
are the blackbody curves at 135 °K and 
120°K respectively. 

6-13 The raw spectra of Mercury and the Sun. 123 

-6—1-4 • - -Di-u-rna-l— path~of- the~Sun about Mer ctiry, T26~ 

draira to scale. The relative positions 
of the Sun are marked at 11 day intervals 
with the planet held as a fixed reference. 
Planeto-graphic longitude are indicated 
for Mercury. (Taken from Soter and 
Ulrichs, 1967) 

6-15 Two coordinate systems on the surface of 129 

Mercury. The unprimed system is the 
"solar system” with the Z-axis pointing 
towards the Sun. A is the subsolar point. 

The primed system is the "earth system" 
with the Z’-axis pointing towards the 
earth. A’ is the subearth point. 

6-16 The final Mercury spectrum, corrected 135 

for solar temperature, with a number of 
blackbody slopes shown to match. 



LIST OP TABLES 


TABLE TITLE PAGE 

2-1 The Value of t for Different Matrices 18 

2- 2 Comparison of Three Different Grating 34 

Spectrometer 

3- 1 A Brief Description of Each Filter 48 

6-1 The Main Constituents of The Jovian 105 

Atmosphere 


xii 





W4^#A:^M. 

Aiwr~ 




xlii 



CHAPTER I 


HISTORICAL INTRODUCTION TO HADAMARD 
TRANSFORM SPECTROMETRY (HTS) 

The idea of modulating or encoding the optical output of 
a spectrometer goes hack to the original work of Golay (19^9) 
and Pellgett (1951). Its purpose is to allow many different 
wavelengths of radiation to fall on a, detector simultaneously, 
and thereby to increase the signal-to-noise ratio (SNR) of the 
resulting spectrum. This improvement comes about because each 
element of the spectrum is effectively viewed a larger fraction 
of the total available observing time. One idea is to encode 
or modulate each spectral wavelength exiting the spectrometer 
output with an audio frequency that contains the optical wave- 
length information. The use of a conventional wave analyzer 
then allows recovery of the original optical spectrum. There 
are many variations of this technique. 

In 1968, Ibbett et a2 and Decker ^ al independently 
suggested the use of sequentially stepped multiplex spectro- 
meters. In both systems radiation enters a dispersive instru- 
ment through a single slit and is analyzed at a number of exit 
slits. Decker ^ al pointed out that two constraints should 
be imposed on the encoding scheme: 

(1) To obtain the optimum signal to noise ratio, each spectral 
element should be viewed during exactly half the step positions. 

(2) To impose the smallest dynamic range requirements on the 
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detector amplifier system, each step position should pass 
light from exactly half the spectral elements. They also 
worked out a scheme fhat satisfied the two constraints for 
masks having elements m=4n+2, where n is an arbitrary integer 
or zero. Ibbett e^ ad introduced the Hadamard pattern for the 
mask. As discussed below, this is a pattern based on a set of 
binary orthogonal matrices first studied by the French mathe- 
matician Jacques Hadamard, Ibbett e^ ad also described the 
application of their scheme to a real time computer aided 
measurement. 

In 1969, Sloane e;t ad worked put a number of binary cyc- 
lic coding schemes for multiplex spectrometry and evaluated 
the performance of each scheme in terms of a linear, least 
mean square, unbiased estimate. These schemes include a Hada- 
mard matrix H, and various modified Hadamard matrices , -which 
these authors refer to as G matrix and S matrix (Fig. 1-1). 

A Hadamard matrix H of order N is an N x N matrix H„ 
of +l's and -I’s which satisfies: 


where is an N x N unit matrix 

A modified Hadamard matrix G of order M is a partitioned 
matrix from the H matrix: 


H 




(M = N-1) 


1 



3 


H 


+ + + + 4 - + + + 
+ + +- -- +■• 
+ +- “- +- + 
+ *"““+~ + + 
+ -- +- 4-4-- 

4 -- 4 -“ 4 - 4 --- 
4 ' 4 -~ 4 " 4 -~ ~ “ 
+ - 4 - 4 ---- 4 ' 


0®^ pace s 

^oit quality 


G 


+ 4 - ~ “ 4 " ~ 

4 " “ ~ “ 4 - ^ 4 " 

- -~ 4 -— 4 ’ 4 - 

--+- 4 - 4 -- 
-+- + +-- 
+ - + + — "- 

— + + — — — + 


0 0 1110 1 
0 1110 10 
1110 10 0 
110 10 0 1 
1 0 l‘o 0 1 1 
0 10 0 111 
10 0 1110 


Figure 1-1. An 8 x 8 Hadamard matrix and two 7x7 cyclic 
matrices that can be derived from it. 
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where the first row and first colunm of H are all +i’s. A 
feature of the G matrix is that it can be written in cyclic 
form — a factor which we will show to be o.f considerable 
practical importance. 

A modified Hadamard matrix S_ is a matrix obtained from 
G by replacing +l’s by O's and -I’s by +l^s. 

The properties of H, G, and S will be discussed in sec- 
tion 2-6. 

VJhen we talk of encoding by means of a Hadamard matrix 
we have the following in mind. A mask is used to modulate - 
open or close - a series of entrance and exit slits in a spec 
trometer. If a certain slit location is open, we can desig- 
nate it by a +1,; if it is closed we can designate it by a 0; 
if it can be used to subtract from the signal incident on the 
detector, we designate it by -1. The sequence of +l's and 
“ -I’s characterizing a mask in a given modulating position 
corresponds to a row of a matrix. The whole set of mask posi 
tions corresponds to the set of rows of the matrix. If the 
sequence of mask patterns corresponds to the rows of a Hada- 
mard matrix we say we are encoding with a Hadamard pattern. 

Sloane e;b ^ also introduced the idea of using a cyclic 
matrix for coding masks. This greatly decreases the experi- 
mental cost and facilitates operation, since any N slits of 

a single mask— ^2N-1 -slits long can“be used to provide ~bne of 

the required mask patterns. 

The first single entrance Hadamard spectrometer (HTS) 
was .built by Decker and Harwit (1969). The spectrometer had 
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a single entrance slit and 19 exit slits. The exit mask was 
stepped manually. The authors used this spectrometer to take 
the spectrum of the mercury vapor 1.7u band to demonstrate the 
Hadamard transformed spectrum’s fidelity and freedom from sys- 
tematic errors. 

With 19 exit slits, the HTS had a theoretical signal-to- 
noise advantage of 2.18 over the conventional spectrometer, 
which is rather hard to verify experimentally, Decker ( 1971 ) 
therefore preceded to build a 255-slit HTS. In this spectro- 
meter, the radiation, after being decoded by the exit mask, 
exits along the same path it comes in. This reverse pass de- 
disperses the beam and allows it to be brought to a focus at 
the entrance plane. Thus the dimensions of the focused image 
are roughly the same as the dimensions of the entrance aperture, 
and the detector size can be minimized. This is important since 
sufficiently large detectors sometimes do not exist, and if 
available tend to be noisy . Decker experimentally verified 
the theoretically predicted multiplex advantage of an HTS. 

DeGraauw and Veltman (1970) were the first to use an HTS 
for astronomical work during the 1970 solar eclipse. Houck 
et al (1973) subsequently used an HTS to obtain near infrared 
spectra of Mars from airplane altitudes. 

Besides putting an encoding mask at the exit plane, one 
can also put another encoding mask at the entrance plane of 
a spectrometer. In that way the radiation is modulated at 

' N 

both the entrance and exit apertures. Harwit ^ al (1970) 
worked out this scheme of doubly multiplexed dispersive spec- 
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trometry. The double multiplexing scheme allows one to in- 
crease the total amount of radiation that can be transmitted 
through a spectrometer. Furthermore,, by a proper reduction 
of the data, one can also obtain a one dimensional picture 
of the source at the entrance plane. For a spectrometer of 
m entrance slits and n exit slits, one needs m x n data points 
to recover m spatial spectra, with each spatial spectrum con- 
taining n spectral elements. For a homogeneous source one 
does not need the spatial information, so (n + m - 1) data 
points will be enough to recover the spectra. Karwit ^ ^ 

( 19 7 a) describe tv?o schemes for recovering the spectrum with 
(n + m - 1) data points. 

In 1975 Tal et al (1975a) finished the construction of a 
doubly coded Hadamard transform spectrometer. The spectro- 
meter -has 15 entrance slits and 255 exit slits, whic h ca n 
simultaneously obtain 15 spatial spectra, each having 255 
specrral elements. Tai ^ (1975 d) went on to give an ana- 
lysis of the errors in Hadamard spectrometry caused by imper- 
fect masks. 

Besides coding the radiation at both the entrance and the 
exit aperture, one can go one step further and use a two di- 
mensional mask at the entrance aperture (Harwit, 1971). 

This yields a two dimensional picture at the entrance aperture 
- where each spatial point -at the entrance ha‘s its' own~'speictrvun; 
To put it a different way, one obtains a two dimensional pic- 
ture of the source at the entrance aperture for each color of 
the spectral elements. 
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Harwit (1973) experimentally verified the operation of 
imaging spectrometry, and Swift et ^ (1976) constructed the 
first Hadamard imaging spectrometer. 

There are other discussions of Hadamard transform spec- 
trometry in the literature, mostly of theoretical aspects. 
Nelson and Predman (1970) give a more complete theoretical 
treatment of Hadamard matrix encoding. They also rediscovered 
a theorem due initially to Hotelling (19^^) showing that the 
Hadamard matrix is the best design for a singly coding mask. 
Sloane and Harwit (1976) show the connection between Hadamard 
spectrometry and the mathematics of weighing designs in sta- 
tistics. 

There have been various comparisons of Hadamard transform 
spectrometry with other spectrometry, Larson £t ^ (197^) 
makes a theoretical comparison of singly multiplexed Hadamard 
transform spectrometers and scanning spectrometers. They pre- 
sent a general mathematical framework for the comparison of 
relative performance and also verify their prediction by com- 
puter simulation of various characteristic spectra. Their re- 
sults show that where the noise level is constant and indepen- 
dent of the Incident photon flux, the determined multiplex 
advantage is /N/2 , as predicted by Pellgett (195D . This is 
usually the case in a low energy region, such as the infrared. 
For a noise level that is signal-dependent, such as in the DV 
energy region, the detector is characterized by an output with 
statistics approaching a Poisson distribution and variance 
therefore proportional to the input signal. In that case the 
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HTS technique will be advantageous only for spectra that are 
characterized by a few well-defined and intense peaks on a 
very low intensity background. For spectra with high back- 
ground, for dense spectra, or for spectra having very weak 
spectral features, the HTS will have no advantage over the 
conventional single slit (SS) technique. 

Hirschfeld and Wyntjes (1973) compare Fourier transform 
and Hadamard transform spectrometry. They also describe 
various limitations of Hadamard -transform spectrometry. This 
paper was followed by an exchange of notes between Decker 
(1973) and Hirschfeld and Wyntjes (1973) in the Journal Applied 
Optics in v;hich some of these limitations are disputed. 

These papers concern themselves with a number of practical 
matters on which opinions can vary. Here we mention these con- 
troversial papers mainly for completeness. Their contends will 
— ^b'e 'dl'S'cus*sie'd~farth^“ b^Ibw^i 

V.^yatt and Esplin (197^) analyzed the effect of band width 
on noise equivalent power (NEP) for multiplex spectrometry 
with cryogenically cooled, cooled-background extrinsic long 
wavelength infrared detectors. They find that the NEP is 
directly proportional to band width, so multiplex schemes that 
require increased band width are not of real advantage. They 
further conclude that doubly encoded systems that are based on 
m--+ n 1 measurements would have a- real, throughput a'dvantage 
‘ Various other aspects of Hadamard matrices and Hadamard 
transform spectrometry which have not been mentioned above are 
covered in articles by: Baumert, Pratt ^ (1969)» Hirschy 
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et ^ (1971), Allen et ^ (1972,1973), Kowalski ^ al (1973), 
Planky ^ ^ (197^), Oliver et ^ (1974). 

In this thesis Chapter II will describe the mathematical 
properties of Hadamard matrices and their application to spec- 
troscopy. Chapter III describes the Hadamard transform spec- 
trometer, and gives results on laboratory performance. 

Chapter IV gives a comparison of Hadamard transform and Fourier 
transform encoding in spectrometry. The output of an HTS is 
fed into a mini computer. The computer performs a real time 
inverse Hadamard transform to recover the spectrum. Chapter V 
describes the algorithm and programming of inverse Hadamard 
transform* Chapter VI discusses observational results and 
their interpretation. 



CHAPTER II 


HADAMARD MATRICES 


(A) Welp;hlng Designs 

In order to understand the mathematical advantage of 

Hadamard transform encoding, let us look at the following 

examples (Sloane et 1976 ). 

Suppose four objects are to be weighed, using a spring 

balance which makes an error e each time it is used. Assume 

2 

that e is a random variable with mean zero and variance a . 

First suppose the objects are weighed separately. If the 
unknown weights are measurements are 

n2> and the errors made by the balance are ej, 62 > 

64, then the four weighings give four equations: 

ni = ’f'l* + ei ri2 = 'f’z + ®2 

TI 3 = i|)3 + ea ri4 = 'f'4 + ®4 

The best estimate of the unknown weights are the measurements 
themselves : 

A 

-J'l ~ 

^2 = H2 = ’^'2 + ®2 


T r\ 
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These are unbiased estimates : 



Ei{/i = 

’/'i 







Ei|»2 “ 



(E denotes 

expected value) 


with variance 

or mean 

s quar e 

error 





A 

2 

2 

2 





E(i|ji - 

^1 ) = 

EcJi 

= ff 




On 

the other hand 

, suppose the balance is 

a chemical 

balance 

with two pans, 

and the four weighings ; 

are made 

as 

follows 

• 

• 








m = 

+ 

ij)2 

+ '/'3 + 

11*4 + 

ei 



t2 = 

4 'i - 


- i |'3 

+ 

©2 



ts = 



- i |'3 





TUt = 

iPi - 

4>2 

- ij ?3 + 

1|)4 + 

64 ( 2 

-1) 

This means that in the 

first 

weighing 

all four 

objects 

are 

placed in the 

left hand pan. 

and in the other weighings 

two 

objects 

are in 

the left pan 

and two in 

the right. (Note 

that 


the e are independent of the weights on the balance. This 
point is crucial). It is easy to solve for 4>2 9 'l'3» 
as long as the coefficient matrix for is not singular. 
Thus the best estimate for t/;i is 
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'i’l 



+ ri2 

+ 

+ TI4) 


hiei 

+ 62 

+ 


©If ) 


p 

The variance of Ce, here C is a constant, is' C times the 
variaince of e, and the variance of a sum of independent ran- 
dom variables is the sum of the individuals variances. 


Therefore the variance of i|>i(and also of 
a 


16 ii 


’f’4) is 


Weighing the objects together has reduced the mean square 
error by a factor of 4 . In effect the signal to noise ratio 
(SNR), which is given by the root mean square (rms) error is 
reduced by a factor of 2. 

Finally, suppose the balance is a spring balance with 


only 

one pan, so 

only coefficients 0 and 1 can 

be used 

. A 

good 

method of weighing 

the four ob.lects is: 




Til = 


i|»2 + t |»3 + ip 4 + 

ei 



T12 = 



62 



t 3 = 

’I'l 

+ ’f '3 + 

63 



Titi = 


+ 1|;4 + 

64 

(2-2) 

In this case the 

9 0 

2 

lift 

variances of i/ji, ^3^ ^4, are 

2 

7 o 

Q » 


7 ct 7 ct 


9 ’ 9 

vious case. 


respectively, a smaller improvement than in the pre- 


The theory of weighing designs is of immediate interest to 


multiplex optics, -since the simultaneous measurement of the 
intensities of different bundles of rays is completely ana- 
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logons to the simultaneous weighing of different groups of 
weights. In measuring the intensity 'of radiation passed ' 
through slits in a mask, we are effectively 'weighing’ that 
radiation. 


(B) General Mathematical Formulation 

One can put the problem into a more general form. 

Let be the i^ unknown, be the measurement, e^ be 
the error associated with the J measurement. Let w^^ be 

4* Vi VVi 

the weighing coefficient of the J measurement with the i^ 
unknown . Then 


ji^i 


i=l' 

. 3 = 1 ' 


•n 

•n 


(2-3) 


In matrix notation: 

a = W + e (2-4) 

With the notation < > for ensemble averages, the error e^ has 
the following properties : 

(1) <e^> = 0 

(2) e^ is independent of jj|_ 

(3) = 0 if errors are assumed to be 

uncorrelated. 

= if i=^ 

The problem now is the following; (!) For a particular coding 
matrix W, what should be the decoding matrix A? i.e. What is 

•ys 

A such that £~A ti_ where is an unbiased estimate of 



Ik 


(ii) What is the best choice of W (or A) that will minimize 
the error of measurement i.e. What is W such that 

2 = < S > 

j=l ^ ^ 


is a minimum. 

In the absence of noise, i.e. n=0, it is clear from (2-^) 


that 


1 

and therefore 


A 

and 

A 

i 



i i = i 


In the presence of noise, 

= A W + A _n 

and with the assumed properties of the noise <ii)^>= one ob- 
tains 

A 

<ji_> = AW + A<£> 

= A W i 

Assuming no prior knowledge of the unknowns, one may use the 

A 

unbiased condition jb. This again implies 

A = W”^ 

So with the assumed properties of noise and unbiased con- 
dition, the decoding matrix is just the inverse of the coding 
matrix. 


One still has to find a coding matrix which will minimize 
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the uncertainty, e , 

The second question can be solved in the following way: 

A •4» Vi 

Let Ti^ be the i measurement in the absence of noise, then 
for each measurement 


"i = * 

2 


= Hi + 




2^2 

+ ••• + ^3n"n 

= Ajl (Hi + 

e^_) + 

*32 ^'>2 * ^2'> * ' 


= (AjiOi + ••• + + (Aj2 ^2 +•••+ 

= i|j. + noise 


The mean square of the noise term corresponding to the j 
unknown is therefore 


th 


■J 


= (A 


n 




= i? (j2 (2-7) 

■ 3 

where 

and Aj represents the improvement in the SNR for the weighing 
design, compared to the SNR for individual weighings. 

Hence, the problem of maximizing the signal to noise 
ratio becomes the problem of minimizing 5 or A^ (Nelson and 
Predman (1970)). 

Sloane ^ al^ (1969) 'independently developed an expression 
for c/a^ where 
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e/a^ = Trace 

T 

= Trace A A (2-9) 

Equation (2-8) and (2-9) amounts to the same thing be- 
cause it can be seen very easily that 

m n 2 

Trace A A = z . Af 
0=1 j 

The question of minimizing e had been answered by Hotelling 

(Hotelling, 19^^) and rediscovered by Nelson and Fredman. 

Hotelling has shown that for any choice of mask W with |W. j<l, 

2 ^ 
the are bounded by e^> and that it is possible to have 
_ 2 - 

e.= — - for i=l,... N if and only if *a Hadamard Matrix of the 

IN N 

order N exists (by taking This leads to the discussion 

of the Hadamard Matrix. 


(C) Hadamard Matrix 

* 

A Hadamard matrix of order N is an NxN matrix Hj^ of +l’s 
and -I's v;hich satisfies: 




(2-10) 


where is an NxN unit matrix. 

A. Hadamard matrix has following properties: (Golomb(1964) ) 

(1) Its row vectors (or equivalently, its column vectors) 
are mutually orthogonal. 

(-2-) The Hadamard' p'roperties wiIl~not be disturbed by: 

a. Interchanging rows, 

b. Interchanging columns, 

c. Changing the sign of every element in a row, or 
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d. Changing the sign of every element in column. 
These properties enable the first row and column of 
every Hadamard matrix to be normalized to contain only +l’s. 
If _G represents the remaining M x M matrix (M= N-1), then H 
can be partitioned into 


1111 

1 


1 


(M= N-1) 


It is conjectured that Hadamard matrices exist for all 
multiples of four. Purther, if one of the following conditions 
is also satisfied, • 

(1) N = P + 1 P prime 

(2) N=P(P+2)+l P and P + 2 prime 

(3) n = 2^ m an integer 

then G can be made cyclic. That is, the (J + 1) row can be 
generated by shifting the J^ row one position to the left. 

For example, when N=8, we have matrices of the form shown in 
Pig. 1-1. Note that H and G are symmetrlx. 

Another choice for W is the matrix S obtained from G by 
replacing +l^s by O’s and -I’s by +l*s. 

The properties of 'the H, G and S are discussed by Sloane 
et al . If rows 1 and J are any two rows of H, G or S, it can 
be shown that their dot product is: 
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In : row i * row J 


In Gj^ : row 1 * row J 


In Sj^j ; row i • row j 


The inverse of each matrix is: 


0 

N i=j 


-1 

M 


M/ij i^j 
M/2 i=J 



“IT • G"*^ 
N % " % 


M+1^^ ^ 


f5+l(2S^ 


where J is a M x M matrix consisting entirely of -l*s and 
N= M+1. 

Table 2-1 gives the value of A for different matrices. 

The matrix ^ represents the weighing scheme v^eighing each 
obje ct s e para tely... This corresponds to-a convent ional^s'lngre” 
slit spectrometer or to a wedge filter monochromator. • 


MATRIX A 

I 

H 

G 

S 


Table 2-1 


ELEMENTS OF A 

1,-0 

1,-1 

1,-1 

1 , 0 

cosine squared 
functions 


1 

/N 

/S(2-|)-^ 

N 2 
¥ N+l 


(FOR LARGE N) 
1 

/N 

if 

/jT 


i! 


* F is the Fourier Transform case which will be disscussed in 
Chapter jy. 
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If the number of measurements N is a multiple of 4, and 

the matrix coefficients are +1, the best weighing scheme is 

the Hadamard matrix H. t will be reduced by a factor — 

/N 

compared to weighing the unknown separately. This is the 
maximum advantage a weighing scheme can obtain with weighing 
coefficient |W^^j<l. If N is- not a multiple of. 4, or if the 
weighing coefficients are O’s and I’s, it is not possible to 
simultaneously minimize and some other criterion 

must be used (Sloane 1976) ♦ Also the errors are uni- 

formly larger than for the H-matrix, as shown for the G and 
^ matrices in Table 2-1 above. 

It is' interesting to see that a spectrometer using the 
Fourier Transform, such as a Michelson interferometer* has a 
multiplex advantage a factor of lower than the H-matrix 
and a factor of /2 lov;er than the S-matrix encoding instru- 
ment. 

Following are some computer simulations for S-matrix 
transformations with various inputs (See Pig. 2-l(a) to (c)). 


INPUT 

1. Constant 

2. Hadamard code: representing 

single line emission 

3. Single line: 1 at the 1^”^ 

element and 0 for the ‘rest. 
This represents an unknown 
impulse coming in during the 


OUTPUT 
Constant 
Single line 

Hadamard code. 

Note, unlike the mono 
chromater, the error 
propagates to other 
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INPUT 

observation. 

4 . Sine wave . 

5 . Square wave . 

The input is not a perfect wave 
because we have an odd number of 
elements. The input values have 
amplitudes 0 or 1 for each element. 

6. Straight line at a slope 1/255* 

This may represent a shift in 
baseline. 

(D) Optical Realization of Hadamard Encoding 

We have made use- of (modified Hadamard) S-matrices in two 
optical instruments: One is _a Ha damar d tra nsform spe_ct-iiO-me_t-er — 

having an encoding mask at the exit aperture. The other is a 
doubly encoded HTS which has encoding masks at both the en- 
trance and exit apertures. 

S codes can be used for both the entrance and exit masks 
for the HTS, with +1 standing for an open slot through which 
radiation is transmitted, and with 0 standing for a closed 
slot where radiation is blocked. 

The cyclic property* of the S-raatrix is very desirable, 
for thSi only a single mask 2N-1 slots wide need be constructed. 
Successive encoding positions are generated by stepping the 
mask one slot width along its length. , This avoids the con- 
sturction of N masks with slots. 


OUTPUT 

elements. 

sine wave with dif- 
ferent phase. 

Not a perfect square 
wave . 


Refer to figure (2-lC). 




(a) Input 


(b) Output 


Figure 2-lb . Hadamard transfo 




Although H and G matrices have better coding efficiency 
than the S matrix, they introduce technical difficulties when 
used for spectrometer coding schemes. To utilize H and G ma- 
trices, one must measure the reflected as well as the tran- 
smitted radiation. In this case +l’s represent reflecting 
slots and -I’s represent transmitting slots. Therefore a mi- 
nimum of two detectors must then be used, one in a subtracting 
mode, the other in the normal mode. The use of two detectors, 
however, increases the noise. Furthermore, the H matrix, with 
all elements +1 in the first row, makes the dynamic range of 
the detector system change by about a factor of two. Also its 
lack of the cyclic property does not allow one to generate the 
rest of the masks by the simple stepping technique mentioned 
above . 

1* Hadamard Transform Spectrometer (HTS) 

For tlie singly encoded HTS, we use the following optical 
arrangement (figure 2-2). Radiation passing through the single 
entrance slot is rendered parallel and directed to the dis- 
persing element. The dispersed radiation is then de-collimated 
and focused upon the multi-slot mask at the-' exit aperture. 

The spectral elements transmitted by the mask pass through 
suitable post-optics and are collected onto a detector. One 
then makes N (in our case N=255) measurements by sequentially 
stepping the mask N times. The inversion procedure \{)=S ti, 
recovers the spectrum. Figure (2-3) gives a 255 cyclic S- 
matrix code for the exit mask. 

Theoretically, the mean square error in the spectrum 




single entronce »Ut 


Eni) mask 
In posilionj 


Figure 2-2. 


Schematic representation of HTS using the 
S code. 








00000 00101 
01001 11110 
01011 11110 
01001 01101 
00101 11010 
01101 01000 
01001 00011 


|ifl!ll;ll mil IIB 


10001 

11 101 

00001 

loiqi 

OHIO 

00001 

11110 

01101 

non 

OOOlil 

00111 

00111 

iniio 

non 

11100 

00100 

10001 

11101 

10010 


mil 11001 00001 
10001 01011 00110 
10010 10100 10100 
10001 10110 00010 
00110 10011 Oipil 

01001 . loood 00111 


Figure 2-3 . The 509 exit silt mask and the 255-element 
S code . 
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given by a monochromator is a^. For an HTS using the S code, 

this error is (sioane et al, 1969). Hence the rms 

gain in S/N for the HTS is G= , For N=255, 

N-’ ^ ^ N 

G'\'8.0 (figure gives Decker’s results. The experimental 
gain was measured as 8.0±0.3). Note that the scale in figure 
(2~^) are in arbitrary units, and the zero point appears to be 
shifted between parts (a) and (b). 

2 . Doubly Encoded Hadamard Transform Spectrometer (DHTS) 
Figure (2-5) is a schematic representation of an optical 
system which has a number of entrance as well as exit slits. 
Instead of passing radiation through only one entrance slit, 
a mask M slits wide is placed at the entrance aperture. 

Radiation passed into the spectrometer through different com- 
binations of open and closed slits. The dispersed radiation 
at the exit plane is analyzed in the same fashion as in the 
HTS. Encoding is accomplished by sequentially stepping one 
of the masks through its N different positions for each posi- 
tion of the other mask. 

In a DHTS the entrance aperture is modulated by a P x P 
S matrix. 

Let t>e the P x P matrix whose rows represent P 

different entrance masks, e. =1 for open slots and 0 for closed 

ir 

slots (l<i<P, l<r<P) . Similarly let x“Xj_j represent the exit 
mask. When the entrance mask is in position i and the exit 
mask is in position j , the detector measurement is 



Intensity in Arbitrory Units 
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Mercury Emission Spectrum - Monochromator 



Mercury Emission Spectrum- HTS Spectrometer 



Wavelength, (>*) 


Figure 2-1} . Comparison spectrum of the mercury emi- 
ssion lines in the 1.4-1. 8pm region: (a^ bottom) 
as .obtained in the- Hadamard-transf^orm -mode,— (b, -top) 
as obtained under identical conditions using the 
same optical system as a conventional monochroma- 
tor. This figure is taken from Decker, 197115 • 





Entrance mosk 
in position I 


Exit mask 
in position ] 


Figure 2-5. Schematic representation of DHTS. 
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"ij “ "ir <2-n) 

Where is the spectral element produced by radiation pas- 

X’S 

sing through the r^^ entrance slot and the exit slot, 
is the noise in the (i,J)^^ measurement j it has the following 
properties: 








If the instrument has no optical magnification, the spec- 

t h 

trum of radiation that passes solely through the r entrance 
slot to first order is shifted by r spectral channels from the 
spectrum passing solely through the first entrance slot. 

Hence, only P+N-1 distinct spectral elements exist: 


(P—l ) , « • . • S * 


'N-1 


where 


t 1) ~ ij; 

’^rs r-s 


(t=r-s) 


In matrix notation one may write 


£ i X X 


(2-12) 


Employing the same analysis as one does for the HTS, i.e, 

A 

using the unbiased condition <i^>=£ and the properties of 
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one obtains 

7 _ -1 , Ts-l 

where 

’^-N+1 

i = ’i'-N 

■ « 

’i'p-1 ip„2 •“: *♦ ’^'-N-p+2 

A. 

Each row 1 of represents a spectrum at the exit mask 

til 

for radiation that enters the instrument through the i en- 

+• Vi 

trance position. Hence the J diagonal gives a one-dimensional 
spatial picture across the entrance aperture for the spectral 
element j . 

A 

One may obtain an average spectrum ij«. by averaging all the 

- 1 

elements in each diagonal. 


As 




P 

s 

r=l 




I*, 


r-t 


t > 0 


- 1 P , 

N-l't r 


t < 0 


Harwit et ^ (1970) showed that if both the entrance and 
exit masks are S matrices and we define 


a 


2 

t 






> 
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Then 



16 W ^-1 -2 

(N+1) H-lt 1 


^ l6g^ 

(N-|t|)H2 for N large 

where 

t = -(N-1), (N-1) 

If the total mean square error for the unknown is 

t=^ 

2 

e = S 0-2 

t ^ 


(2-13) 


where one sums only the central element, then for the S-code 

- /t 2 22. 18 , ^/I'l 

e “ — + 0( — ) N large 

where = constant ^ . 

There are two points that should be made about the DHTS. 

(1) It has not been shown that Hadamard codes are the best 
codes for such an instrument. In fact some evidence suggests 
that Hadamard codes are not precisely optimum for this "two- 
ended” operation (Harwit ^ ^,197^b). (2) For PxN data points 

the spectrum yields only P+N-1 spectral elements, plus a one- 
dimensional image. It is also possible to reconstruct "the 
- (-P+N-1-) elements with (P+KT-l)" data" po'ints only THarwit eib al , 
197^b). Fig. (2-6) gives a 15 element cyclic S-matrix code for 


the entrance mask. 




1 i 1 1 O 1 O 1 1 O O 1 o o o 


Figure 2-6. 


The 29 entrance slit mask and the 15-element S code. 
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Table 2-2 compares three different grating spectrometers. 
The first column represents the conventional single entrance 
and exit slot instrument. N measurements are made in time T, 
with a mean square error in each. The second column is 
for a singly multiplexed instrument with an exit mask -S, and 
is taken from Sloane ^ ^ (1969)* The last column is for the 
doubly multiplexed system, using Equation (2-13) for c^,_and 
has been multiplied by a factor of N to allow for having to 
make measurements in time T. 

Table 2-2 

NO MASK SHT DHT 



1 


N 

^ /2272 


(E) Errors in Hada.roa-rd- S-peotroscnpy 

During the manufacture of the masks, whether by deposi- 
tion of metal or by removal of metal through an etching process, 
it is possible to obtain a systematic error that leaves each 
of the opaque portions of the mask either too wide or too 
narrow by a -fixed amount. This will cause a systematic vari- 
ation in signal passing through the slit. For example if the 
open slit is too narrow by a fixed amount e, the light passing 
through an open slit position will be 

when the open slit is bounded by two open 


slits . 
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I^(l-e) when the open slit is boimded by one open 
slit and one closed slit- 

I^(l-2e) when the open slit is bovinded by two closed 

slits. 

A similar analysis holds for open slits that are too wide, 
except that the minus sign in these expressions is replaced 
by a plus sign. 

The spectrum of a single (spectral line resulting from 
such imperfect masks) is remarkably simple (Tai e^ a^, 1975 b). 
Independent of the particular S-matrix mask to be used, there 
are always precisely four false blips present in the final 
spectrum. The amplitude of these blips is always the same for 
a fixed narrowing or widening of the transmitting slits. Two 
of the blips always surround the main spectral line, and a pair 
of adjacent blips always are some distance removed from that 
line. For the 255 element S-matrix, these two are located 24 
and 25 elements away to the left of the parent line. The am- 
plitude of the displaced blips is positive when the transparent 
slits are too wide and is negative when the slits are too nar- 
row. In contrast, the two 'blips surrounding the parent line 
always are positive. -Figure (2-7) shows the negative features 
accompanying the 1.7 ym mercury vapor doublet and the computer 
simulation of a pure spectral line input and its distorted spec- 
trum. For the general case the reader may refer to Tai et al 
(1975 b). 



36 


ORIGINAL PAGE IS 
• OP POOR QUALITT 



Figure 2-7 . (a) Spectrum of the 1.7vm mercury vapor doublet 

showing negative peaks to the left at the emission peaks • 

(bj Shows the response we would obtain to a single spectral 
line with a perfect mask; Cc) Shows the response for a single 
radiation simulated as passing through a mask 
with slits too narrow because each opaque mask element pro- 
trudes into the adjacent transparent slot by a tenth of a 
slot width; (d) Shows the effect of simulating slits that are 
systematically too wide. Note that the main spdctral line 
has been placed in different positions for the synthetic runs 
(b),(c), and (d). 





CHAPTER III 


INSTRUMEIWATION 
(A) Multislit Spectrometer 

A conventional spectrometer has four e^ssential elements, 
an entrance slit, a dispersive device such as a grating or 
prism, a set of imaging optics, and an exit slit. 

Such a spectrometer has two important parameters. The 
first is the "resolution” R, which is a measure of how well 
the spectrometer can separate two neighboring lines. The 
second parameter is the "throughput" E. This is a measure of 
the light gathering capability of the system. The tvro para- 
meters, R and E, are lumped together into what is called "lu- 
minosity" L, defined as (Vanasse, 197^ 

- L = E • R - 


For a conventional grating spectrometer having a grating 
of area A given by VJH, where W and H are the width and height 
respectively of the grating, the throughput is determined by 
the product of A with the solid angle n subtended by the slit 

s 

at the collimating mirror (or lens). The solid angle is given 


by 




w*l 



where w and 1 are the vridth and height respectively of the slit 
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and F Is the focal length of the collimating mirror 

E -x. W • H • — . _(3“1) _ 

- - - - p2 ~ - 

The resolution of this instrument is 


R 


X 

AX 


= nN 


n . ^ 


(3-2) 


where X is the wavelength, AX the closest wavelength that can 
be separated, n the order, N the total number of lines on the 
grating, ¥ the width of the grating and d the spacing between 
rulings. 

Since d sin a = nX (3-3) 
substitute (3-2) into (3-3) 


W 

^ sin a = 
it 


nX 


(3-^) 


if the slit width is limited by diffraction, which is the mi- 
nimum slit width, then 


w 


PX 

W cos a 


(3-5) 


substituting (3-5) into (3-^) one gets 
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R ~ tan ot (3“6) 

w 

Comparing equation (3-1) and (3-6) one sees immediately 
that, for a fixed grating area W-H, and fixed optical system, 

E is proportional to the slit width w and R is inversely pro- 
portional to the slit width w. This means that an increase in 
luminosity of the system by increasing the slit width is made 
at the sacrifice of resolution, and vice versa. 

Prom (3-1) and (3-6) one obtains 

L E * R 

W • H • i 
P 


Another feature of a conventional spectrometer is that it 
transmits only one narrow spectral* range of light to the de- 
tector, and all other spectral elements are wasted. As a re- 
sult, the instrument is inefficient. 

Within the past two decades there has been much research 
done in an effort to design new spectrometric systems with a 
view to maximizing the luminosity L, and to observe a number 
of spectral elements simultaneously. This can provide a mul- 
tiplex advantage, or a wide aperture advantage. Two quite dis- 
tinct modulation techniques have been employed in the past. 

The -first depends on the wave nature of radiation, and makes 
use of interferometry. The Pabry-Perot interferometer, 
Michelson interferometer, and Mach-Zehnder interferometer 
(Jacquinot, 195^ a I960; Vanasse and Sakai, 196?) are instru- 
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merits of this type. The other technique employs dispersing 
spectrometers in which entrance and exit slits are replaced by 
opaque or transmitting masks. Golay's multislit spectrometer, 
Girard’s Grill spectrometer and Hadamard spectrometers (Harwit 
et al,197i{a) are representative of these instruments. 

A spectrometer, whether interferometric os mask-multiplexed 
yields a multiplex advantage mainly for detector noise or ampli- 
fier noise limited applications. In these cases it can be Shown 
that for N spectral elements, one can achieve of the order 
improvement in the overall spectral signal-to-noise ratio, 

S/N, over a conventional spectrometer (Chapter II). 

For photon noise limited applications, the multiplexing 
advantage is cancelled by the K-fold increase in the photon 
noise attributed to the N-fold increase in the energy falling 
onto the detector. Nevertheless, for photon noise limitations, 
the throughput advantage can still be realized. The large 
throughput will become a disadvantage when the noise is back-!* 
ground noise which increases faster than the noise due to the 
source (Harwit ^ 1974a)* 

In this chapter we will describe the experimental study 
of a Hadamard transform spectrometer (HTS) and calibration in 

the laboratory. Figure (3-1) is the flow chart of the whole pro 
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cess, starting with radiation from the telescope and ending 
- -with the -output of the -computer.- Each' comp’onent 'will be des-" 
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Figure 3-1 . The flow chart of the data taking process, 
starting with radiation from the telescope 
and ending with the output of the compu- 
ter . 
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(B) The Optics 
1. Spectrometer 

Fig. (3-2) .shows the basic spatial design of the spectro-^ 
meter. It works in the Ebert-Pastie mode. Radiation passing 
through the entrance aperture S falls upon the spherical mirror 
M which, in turn, renders it parallel and directs it to the 
grating G. The dispersed radiation is collimated by the other 
half of the spheroid and focused upon the exit aperture S*. 

A 255-slot encoding mask is located at this position. The exit 
focal plane is positioned in such a way that it bisects a 90° 
corner reflection. The corner reflector returns the radiation 
through the spectrometer again and displaces the beam from 
the center of the principal plane to one side. This reverse 
process dedisperses the beam and allows it to be brought to a 
focus at the entrance plane. The dimensions of the focused 
— i-mage—a-re— roughiy— the-“same*~as~the~~dime‘nsTbns^'f“"fhe entrance 
aperture (Decker 1971). These procedures allow one to use a 
smaller detector. 

A diagonal mirror at the entrance directs' the dedispersed 
radiation to the liquid helium cooled post optics. 

The entrance mask can be a single entrance slit with any 
width between zero to 1.5 mm for the 1x255 program, or it can 
be a fifteen S-matrix code with each slit having width .0.1 mm 
for the 15x255 element program. In normal use the height of 
of entrance slit is 3.5 nim. It can be increased up to 10 mm. 

is a spherical mirror with a 49.5 cm focal length. 

On its back it is held in place with three teflon-tippled 



43 



I 


Post Optics 


Figure 3-2 . Optical path through the spectrometer. 

The dedisperser and exit mask are shown 
rotated by 90°. 





screws. Three teflon-tippled springs bear on the front edges 
of the mirror. This design allows one to make slight adjust- 
ments in the position of, M when aligning the- 1-nstrument . -The 
central part of the mirror is blocked off to reduce stray ra- 
diation. 

G is a 75 mm X 75 nun grating with 20 llnes/mm and blaze 
angle 5° 11*. The corresponding blaze wavelength at first order 
is 9.03u. Figure (3~3) is the calculated grating efficiency as 
a function of slit position, at two wavelengths. At 8y the 
energy imaged within one slit width is 80^ of the total. For 
l^u the energy within one slit is 52%. The grating is mounted 
on a yoke which allows it to be adjusted In three mutually per- 
pendicular directions. It is located at 0.82 focal length from 
the primary mirror M. 

All mirrors inside the spectrometer are silver coated with 
A_filLoi;_e_cjti.ve_ c oat i-ng— o-f - -S-iO^-^ — The— r efTeb‘tUVlir5r^f"5 il ver~c oat - 
ing at lOy is better than 97%. 

For a multiplexing spectrometer which has N entrance and 
N exit slits, one wishes to image entrance slits S^, 82* 

t t I 

onto exit slits S^, Sg, such that is imaged onto 

1 I 

, onto at a particular vravelength X . Let 5 be 

the angle subtanded by S. Where S and 6’ are measured from 
the center of M. The grating equation for imaging S onto S* 
is 


Sin a + Sin B 


mX 

a 


(3-7)- 



Diffraction Efficiency 


in 

-=r 



Linear Distance from Center of Central Slit (mm) 


Figure 3-3 . Grating diffraction efficiency as a function of slit position 
at 8 and l4uni. 
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Differentiating (3-7) with respect to a gives 

^ Cos ° (3-8) 

da -Cos- g 

The minus sign indicates that a and g change in opposite direc- 
tions. da is the width of the entrance slit and dg is the width 
of exit slit, so unless a-g. 

In our spectrometer, the exit slit width is 0.1024 mm, 2A% 
larger than the entrance slit width. This is the effect of 
anamorphic dispersion- a magnification produced by the grating 
(3-8), when the lower limit on the slit width is set by dif- 
fraction. The total number of spectral elements that can be 
observed simultaneously is limited by the optical aberrations 
of any particular optical system, which set a limit on the 
total useful width over which the spectrum can be displayed. 

2 . Post Optics 

The 'pbs^~ o^ics consist of a liquid helium cooled Arsenic- 
doped silicon (As:Si) detector, with appropriate optics for 
focusing the radiation onto the detector (Pig. 3-4). Radiation 
enters the evacuated dewar through a barium fluoride v/indow,' 
passing through a filter with pre-selected band-width. The 
filtered radiation then passes through a cooled barium-fluoride 
filter and is focused onto the detector inside the housing by 
a gold coated mirror. A light baffle is partitioned in front 
of the detector housing. 

Below is a brief discussion of each of the cryogenic 


components 
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(a) Barium Fluoride Window 

The barium fluoride window is 2 mm thick and 0,75" in dia- 
meter. - I-t is u^ed to cut off wavelengths- longer than l'4y. It 
has a transmission efficiency of around 90 % out to llv before 
it starts to cut off. At l4.29v (700 cm“^), its transmission 
efficiency is 50 %. 

(b) Filters 

A liquid helium cooled filter wheel with 8 filter posi- 
tions is housed inside the dewar. Table 3-1 gives a brief 
description of each filter. 

Table 3-1 

Peak transmission 

Filter position Range Band-width efficiency 


1 

no filter 



2 

closed 
_(-b looked- b-y- 
aluminum foil) 



3 

8y 

9.9v 

1.9v 

7155 

4 

8.7v - 

ll.lv 

2.4y 

83$ 

5 

llv 

12. 4y 

1.4y 

82 $ 

6 

11. ly - 

13. 8v 

2.7w 

86$ 

7 

8 . 4y - 

15w 

6 . Oy 

85$ 

8 

glass 


shorter than 2y 



Filter positions 1 and 2 are for testing purposes. 

The filter wheel is held^ in place by a spring-loaded screw.. 


The wheel was originally connected to^ the outside world through 
a stainless steel_ rod and could be changed to different posi- 
tions by turning the rod. It was found that the- stainless 
steel rod conducts too much heat from the outside into the 
helium containing can. When the stainless steel is replaced 
by a G-10 Glass Epoxy Lamitex rod, the holding time for the 
liquid helium of the dewar increases from 9 hours to 15 hovirs . 

( c ) Field Mirror 

A gold coated mirror with focal length 7 mm and f/o.^1 is 
used as a field mirror to focus the radiation onto the detector. 
The reflection efficiency for gold mirrors at lOy is over 99?. 
The mirror has 3 degrees of freedom of adjustment, one trans- 
lational and two rotational adjustments. 

(d) Liquid Helium Cooled Barium Fluoride Filter 

This barium fluoride window is also used to cut off the 
radiation longer than' l^y . Although the filter barium fluo- 
ride window cuts off radiation longer than l4>i from the out- 
side world, it will emit radiation of its own because it is 
at room temperature. Since all the interference filters have 
a long wavelength leak between 20]i to 26y, and since the de- 
tector will cut off radiation longer than only, there is 
still radiation from 20y to 2^y that gets into the detector 
as background radiation. The insertion of the liquid helium 
cooled barium fluoride filter eliminates this peak. It was 
found to cut down the background radiation by a factor of four. 

(e) Detector 

An arsenic doped silicon detector with dimension 1.2 mm x 
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3.2 mia is positioned inside a housing, which has a baffle at 
its entrance. 

Arsenic doped silicon is a N-type extrinsic semiconductor. 
When the detector absorbs radiation, free carriers are provided 
for 'the conduction band, thus changing the resistance of the 
detector. The following discussion follows the work of Put ley. 
For further details, one can refer to Putley (1964) and Kittel 
(1966). 

Let a be, the conductivity of the detector 
e be the electric charge of the carrier 
T be the life time of free carriers 
N be the density of free carriers, and 
y be the mobility of free carriers 
Then, 

= NePy (3~9) 

and 

AOp = AJeTpTi . 

AP 

= Ev (3-*l0) 

where 

is the change in conductivity of the detector 
AJ is the number of photons incident in unit time 
AP is the radiation power incident 
V is the frequency of the incoming photons, and 

T) is the quantum efficiency, i.e, nmuber of electrons 


freed per photon. • 
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Since 



-Where Rj^ is. the resistance of the detector 
z is' the length of the detector, and 
A is the area of the detector 

Then 


( 3 - 11 ) 



A 

& 

A 



Al T1 
hv N^eyT 


Let 




where is the detector current, and 

Vg is the bias voltage across the detector 

Then 





(3-12) 


( 3 - 13 ) 


Let 



V 

^ A A?. n 

Rp A hv N^evT 


"d^l 



( 3 - 14 ) 


( 3 - 15 ) 
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where V is the voltage across the load resistor 
o 

is the load resistance 

Then 

= «L.^B ^ 1 

Prom equation (3-15) one can calculate AP from 
(f) Procedures for Alignment of the Optics 

i) Place all components in their respective positions 
and line them up visually. Be sure there is no mechanical 
binding in the mask and in the driving mechanism. 

ii) Using a laser, put the spot from the entrance slot 
on the middle of the grating. It is suggested that only one 
entrance slot be used. 

iii ^ _A^i^t_J^e__grating3-llt_-until_-the_line--.of-di-spersed 

dots exits at the proper position at exit. As the grating is 
rotated, this line of spots should remain level, not displaced 
normal to itself. 

iv) Put one half of the dedispersing mirror combination 
in place. Use a T square to line it up .roughly. At this point, 
use a mercury emission lamp with proper f-number to simulate 
the beam coming from the telescope. A number of colored image 
of the single entrance slot will be seen at the exit. 

v) Adjust the spherical mirror for coarse adjustment 
and the position of the dedispersing mirror as a fine adjust- 
ment to bring the image to a focus on the exit plane. 

Vi) Adjust the grating position in its yoke until the 
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color image from the mercury lamp is parallel to the exit .slit 
length. 

vii) Adjust the angle of the dedisperser so that the 
color image reflected by it is perpendicular to the exit mask. 

viii) Put in the other half of the dedisperser and line 
it 'up at an angle so the radiation falls back upon the‘ grating. 
Adjust with fine adjustment screws so the grating is fully 
illuminated. The image of the grating will appear on itself 
when viewed from the diagonal, ^5° * mirror which diverts the 
radiation to the dewar. 

ix) Place the dewar on the spectrometer using h" spacers 
to represent the thickness of the' dewar bottom cover. Adjust 
the ^5° mirror so that the dedispersed image (with color) falls 
on the center of the filter on the filter wheel. 

x) Turn the filter wheel to position one (no filter), 
so that radiation can fall on the field mirror. Adjust the 
field mirror until the dedispersed radiation impinges upon the 
detector. Be sure that all the light falls onto the detector. 

Be sure that the mirrors accept all the radiation. 

xi) Insert the housing. Be sure that the incoming ra- 
diation is clear of the housing. 

xii) Put in the liquid helium shield, and the radiation 
shield. Put on the nose. Use GE varnish and aluminum foil 

to reduce openings in the baffles so that they will transmit 
only the bright white fringes. 

xiii) When the spectrometer is on the telescope, maximize 
the signal by tilting and rotating the dewar. 
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(C) The Electronics 

The data-taking process is controlled electronically to 
ensure a smooth process. The operator only needs to turn the 
switch on. The spectrometer will then automatically take date 
in, process it, and stop at the end of the transform indicated 
by the operator. All the operator has to do in the whole obser- 
vation is to keep the astronomical object in the beam. Figure 
3-5 shows the block diagram of electronic and computer set 
up. The following are brief descriptions of the electronic 
parts incorporated in the system. 

1 . Alignment Sensor 

The alignment sensor (figure 3-6) is used to synchronize a 
Monsanto electronic counter, and the computer with the spectro- 
meter. The circuit is shown in figure- (3-6). The exit mask is 
continuously moving. When the exit mask is at its starting 
position, i.e. the first 255 slots are at the exit aperture, 
a light pulse goes through an alignment slot on the exit mask 
and is detected by a photo-cell on the other side of the exit 
mask. Two transistors amplify the output light curve and two 
IC chips change the light curve into an alignment pulse. This 
alignment pulse, through the drive unit, readies the counter 
for counting, readies the computer to accept data, and to turn 
on an indicator light showing that the system is taking data. 

One can adjust the starting position of the exit mask by ad- 
justing the intensity of the light. After 255 data points have 
been obtained, the exit mask is at its other end, and another 
alignment pulse turns off the counter and the indicator light. 
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The computer after taking 25i5 readings checks that the indi- 
cator light is off. This step ensures the synchronization 
of computer and spectrometer. 

2 . Drive Unit 

A drive unit (Fig. 3-7) is used to drive the entrance 
and exit masks of the spectrometer. It can be set tp be 
adjusted by pulse streams at 100 Hz, 200 Hz, or 400 Hz. The 
unit performs the following functions; 

(a) It receives an alignment pulse from the sensor 
circuit (18) and generates a pulse to reset the counter for 
counting ( 17) . The pulse will also ready the computer for 
taking the data. 

(b) It drives the exit mask in a continuous mode at a 
displacement rate of one slot for every 4l pulses (l8). 

(c) After each set of 4l pulses the unit instructs the 
the mini-computer to read the integrated signal off the counter 
and then reset the counter for the next data integration. 

(d) After 255 reset pulse the unit advances the en- 
trances mask by one slot by sending the entrance mask advance 
motor 40 pulses (11). 

3. Preamplifier 

The circuit (Pig. 3-8) shows a transirapedance amplifier 
implemented with a Burr Brown operational amplifier. The cir- 
cuit has the advantage of high speed, low susceptibility to 
microphonics, and detector operation at constant voltage with 
a high resistance load resistor. 


Neglecting the voltage noise and current noise in the 
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first order approximation, the current through the detector, 
also goes through the load resistor Rp because the input 
impedance of the operational amplifier can be taken to be very 
large, so 



where is the output voltage, 

V_ is the constant bias voltage, and 
b . 

R is the load resistance 

(D) Laboratory Calibration 

1 . Calibration of Post Optics 

Liquid nitrogen is used to calibrate the efficiency and 
sensitivity of the post optics. Liquid nitrogen is in a dewar 
with black paper along the wall to simulate a blackbody. 
Through a chopping device, the post optics will alternately 
see radiation from the liquid nitrogen and from the room. The 
difference between these two radiations gives the A.C. signal. 
D.C. measurements are obtained by putting liquid nitrogen 
directly in front of the dewar window. The following are the 
results of the calibration: 

Dewar profile: f^ = 7.6 

fy = 10.3 

8.7u - 11. ly 


Wavelength region: 
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Band-width: 

2.4y 

Bias voltage: 

15y 

Load resistor at LH^ 

temperature: 1.2 MO 

A.C. power: 

9x10“^^ watt. 

A.C. signal.: 

272.7 mv 

A.C. noise: 

4 .^yv 

(NEP) A.C. detector: 

5.2x10“^3 

(NEP) A.C. system: 

1 . 25 x 10 ”^^ 

A.C. responslvity ; 

1.2 amp/watt . 

D.C. power": 

3. 05x10“^ watt. 

D.C. signal: 

11.32V 

D.C. noise: 

4yv 

(NEP) D.C. system: 

1.47x10“^^ 

z 

D.C. responsivity : 

2.58 amp/watt. 

Background noise: 

^^^Blip 


=3.5x10'^^ 


where is the D.C. power 

V is the frequency that is assumed to be 
The system is a factor of 3-57 away from background limited. 

2 . Test of the Spectrometer 

A mercury vapor lamp with emission at 1.7y is used with 
the spectrometer to test the computer program. A PbS detector 
and an appropriate blocking filter isolates the 1.7y doublet 
of mercury. The slit width and length for each mask is 0.15 mm 
and 3.5 mm. Fig. (3-9) is the mercury vapor spectrum at 1.7y 
for the 1 X 255 mode. Fig. (3-10) is the mercury vapor spectrum 



Relative Flux 



Figure 3-9 . Spectrum of the mercury vapor 1.7pm line using the 
1 x 255 mode. 
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for 15 X 255 mode. Figure (3-lOa) shows the eighth of a series 
of fifteen individual spectra. Figure (3-lOb) shows an ave- 
rage of all fifteen specrra, and (3-lOc) shoves all 15 spectra, 
each spectrum being displaced vertically from the next. The 
diagonal pattern near the right-hand edge represents a dis- 
placement of the peak between successive spectra. This re- 
presents the actual shift in spectral range between adjacent 
spatial elements. 

Once it v;as clear that the instrument with the computer 
program worked properly in the 1.7v region, the spectrometer 
v;as tested wirh the cryogenically cooled, arsenic-dop.ed sili- 
con detector. The transmission spectrum of polystyrene with 
a soldering iron as the source was obtained. The shap e of the 
filter profile and the transmission spectrum of polystyrene 
shovfed that the instrument worked in the lOv region. 

The wavelength calibration of the spectrometer is obtained 
by comparing the polystyrene transmission spectrum obtained by 
the spectrometer and the spectrum obtained by a Perkin-Elmen 
monochromator. The wavelength calibration is then checked 
against the moon spectra at 9.5v where the atmosphere has 
strong absorption features. 
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Figure 3-10 . 

Calibration spectra obtained 
for the mercury vapor emi- 
ssion lines at 1 . 69 vra and 
l.Tlum: (a) The eighth of a 
series of fifteen individual 
spectra obtained; 


(b) An average of all fifteen 
spectra; 



(c) A display of the fifteen 
spectra, each spectrum being 
displaced vertically from 
the next. The diagonal pat- 
tern near the right-hand 
edge represents a displace- 
ment of the peak between 
successive spectra. This re- 
presents the actual shift in 
spectral range between ad- 
jacent spatial elements. 




CHAPTER IV 


COMPARISONS BETWEEN FOURIER TRANSFORM AND 
HADAMARD TRANSFORM SPECTROSCOPY 


Since the Michelson interferometer spectrometer (MIS) and 
Hadamard transform spectrometer (HTS) both have the multiplex 
advantage and the advantage of large through-put j there are a 
number of comparisons between them in the literature. In this 
chapter the comparisons are carried out in four different as- 
pects: mathematically, computationally, optically and mecha- 

nically. 

(A) Mathematically 

Fourier transforms and Hadamard transforms can be viewed 
as using two different weighing schemes. Appendix A gives the 
mathematical analysis of the coding error for Fourier transform. 

Let 5 be the path difference in a two beam interferpmeter . 
P(v) is the power at wave number v (i.e, the spectral density 
function), v here is taken to be the inverse of wavelength, 
then S(c), the power received for path difference is i 

S(?) = /”P(v) cos^ (2‘ir^v)dv . A-1 

0 

= 1/2P + l/2/“P(v) cos(4ir5v)dv A-2 

o 0 

In the Fourier transform each spectral element can be 
viewed as having a phase modulation by a cosine squared term. 
Its argument depends on the stepped distance and the particular 
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wavelength. In the Eadamard transform each spectral element 
is modulated by a step function of 1 and 0 for a S~matrix 
coding. _ , - - „ - 

One can ask oneself whether these modulating or encoding • 
schemes are equally efficient? 

In table 2-1 we stated that for large N the multiplex ad- 
vantage, or the efficiency, of the H-matrix coding of elements 
1 and -1 is >^17, the S matrix coding of element 1 and 0 is ^ 
and for a single detector MIS the Fourier coding is /k/8. 

This last figure is based on calculations shown in Appendix A. 
The other two' values vjere described by Sloane ^ ^ (1988). 

The E-matrix is an orthornormal matrix. The element -1 means 
"subtract” the radiation, while the +1 element means "add" the 
radiation. No radiation is wasted and thus the coding scheme 
has the highest efficiency. For S-matrix coding half the 
slits are open, 1, letting light pass through; and half the 
slits are 0, blocking the light. Half the radiation is there- 
fore wasted each time and one can intuitively see v?hy the 
efficiency of the S-matrix is only half that of the K-matrix. 

Although the S-matrix is less- efficient, it has the im- 
portant advantage that it is cyclic; that is, the (i+1)^^ 
column of the S-matrix is obtained by shifting the i^^ column 

cyclically one place dovmwards. Instea'd of constructing a 
2 

mask of N slits for N spectral elements, one constructs only 
'one mask with 2N-1 slits. Such a mask has two advantages. 
First, the cost of mask construction is reduced by '^N/2 and 
the design' of the advancd mechanism is considerably simplified-. 
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since the total weight of the mask also decreases as> 'vN/2.-, 
Secondly, it can be self-supporting and therefore permits the 
construction of a spectrometer which requires no transmission 
materials. In operation the mask is stepped one slit width - 
along the length of the mask- for each successive encoding po- 
sition. 

In the Fourier case., equation (A-2) of Appendix A 

S(0 = 1/2P^ + l/2/“p(v)cos(i<TrCv)dv (A-2) 

shows that half the power goes into 1/2P^, the first term on 
the right hand side, which is not modulated at all. This re- 
duces the efficiency by a factor of two as in the S-matrix case. 
The other half of the power is modulated by a cosine term. 

Cosine modulation gives a factor 1//2 because cosines functions 
do not form an orthronormal set themselves. The total Fourier 
modulation efficiency is therefore l/2^/N/2j. Mathematically, 
the S-matrix is a factor of /2 better in SNR than the Fourier 
transform. The true Hadamard code, which cannot be realized 
experimentally, as yet, is a factor of /5 better that the 
Fourier code. 

(B) Computer Requirements 

Both Mis and HTS require a digital computer to decode the 
data. However, in the HTS case this reduces to nothing more 
than a series of additions and subtractions. Hence, as much 
as an order of magnitude in computer time can be gained over 
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the Fourier decoding procedure required by MIS (Decker, 1971). 
In addition, HTS do not have the large zero path-length spike 
which is characteristic of MIS, and hence can be operated 
with a substantially lower dynamic range. 

(C) Optics 

HTS attempts to "liberate" the grating instrument from 
its inferior position and offers a possibility to convert a 
conventional scanning spectrometer .into a multiplex instrument 
at a moderate cost. However, it is also the grating and optics 
that limit the capabilities of HTS. 

1 . Resolution 

The resolution of a grating instrument is 


X 


AX 


mN 

(4-1) 


(4-2) 


where m is the order of diffraction. N is the total number of 
rulings in the grating. W is the width of the grating, d is 
the separation of lines. 

The MIS introduces variable path differences between two 
interfering beams. The resolution is determined by the maxi- 
mum permissible path difference in the interfere gram. 


1 

2x 

(4-3) 

1 

(4-4) 
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where x is the maximum path difference between interfering 
beams. 

5 is the displacement of one of the two mirrors from 
the white light fringe positions and 
Av is the increment in wavenumber. 


Since 


R 


_A 

AX 


1 

XAv 


(Mertz a, P.5) 


X 


A MIS can have a resolutuion as high as '^10^. 


(4-5) 


2. Slit Width 

In the HIS, the minimum usable slit width is determined 
by the diffraction pattern 


w I F (4-6) 

where P is the focal length of the imaging mirror. 

It has been argued that a boxcar profile is a poor match 
to a sine diffraction pattern(Hirschfeld, ^ al, 1973, Mertz, 
1976 b). In the presence of diffraction, the transmission at 
each point of the mask will be a 'complex function of the spec- 
tral distribution, the mask position, and the relative width 
of the nearby transparent and opaque slits. 

One way to correct this is to make the slits, wider than 
the diffraction limit, allowing the mask’s transmission to 
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approach the geometric optical limit. This way, however, not 
only the resolution of the instrument is reduced, but the to- 
tal number N of spectral elements that can be observed simul- 
taneously, also decreases. 

There is another way to look at the same problem. The 
grating is an operator that changes the frequency domain into 
the spatial domain, so that intensity as a function of fre- 
quency, after passage by the grating, becomes a function of po- 
sition. With diffraction effects included, the intensity has 
a new functional dependence on distance. The Hadamard mask 
code and the subsequent decoding process just translate this 
spatial distribution function back into its spectral' domain. 
Therefore, an intensity pattern which is complicated by the 
diffraction pattern, after the Hadamard coding and decoding, 
should still show up the same intensity pattern. Since one 
knows the diffraction pattern for a given optical system, the 

diffraction effect on the spectral intensity distribution can 

2 

be computed and corrected, so that the sine diffraction effect 
would not make the slit width any larger than the diffraction 
limit. One point, however, should still be noted. The cor- 
rection that would have to be applied is wavelength dependent 
because diffraction is wavelength dependent. 

A similar problem exists in MIS (Stev/art P.296) because 
the, move able_ mirror mu^ be stoppe.d ,when its. maximum- displace — 
ment is reached. 

The sidelobes in the interferometric case is of the form 
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suo = ^ (ii-7) 

[D-.[lJ ’ 

(Stewart P.295) 


which has considerably stronger side lobes than the diffraction 
pattern. Here, T is the time for the mirror to travel from one 
end to the other, and ea is the central frequency. Various 
schemes of apodization have been introduced to compensate for 
the side lobes in the interferometric case. 

3 . Multiplex Number 

The total number of elements N that can be observed simul- 
taneously yields the multiplex advantages. In HTS, the total 
aperture size is limited by aberrations largely due to off axis 
radiations. The aperture width is (Hirschfeld, 1973) 

w = P • S (^-8) 

w 

where P is the instrument focal length and S /2 is a factor of 

w 

the order of 0.05'^ 0.1, that describes how far off axis one 
can go before the aberration pushes the individual slit width 
up to the point where no further gain in N is possible. 

Since the minimum slit width, determined by diffraction ef- 
fects, is 'v.i.22Xf, the total number N is 


N 


PS 


w 


1.221f 


and is of the order of 10*^. 

For our HTS at Cornell, we have 
F = 49-5 


(4-9) 



72 


S^'vO.08 (assume medium value) 
lO.Ovim 

f= 7.5 

N-V.250 

Actually, however, Mertz and Flamand (1976) suggested that 
values >>0.1 may be realized in practice. 

The MIS has a very large effective value of N. This is 
its main gain. The total wave number range observed in MIS is 

max min 4 A 

1 

TiS 

where A is the step size of the mirror, a is the resolution in 
wave number. N for the MIS can be 10^. 

^ ♦ Spectral Range 

The MIS also has a broad free spectral range. Its range 
is limited by the beam splitter efficiency which usually varies 
approximately as the cosine of the wavelength. HTS free spec- 
tral range is usually about one grating order. Its free spec- 
tral range can be increased by using order sorting, but this 
increases the technical difficulty. Although the HTS has a 
smaller spectral free range, it can be set to recover only 
those spectral bands of particular, interest throughout any • 
spectral regions (Decker, 1971). This is Impractical with a 


MIS. 
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5. Throughput 

The throughput is defined in section 1-A as the product 
of aperture A and angular acceptance n 
E = An 

For a MIS it can be shown that nR=2n and 

2 = A . ^2L 

^MIS ' R 

where R is the resolution of the instrument. A is the area 

m 

2 

of the interferometer mirror. Typically, A-'-l cm , and for 
R'^'IC^, X 10“^. 

For the HTS, from equation (3-1) one obtains 

F = A 

^HTS g^j.2 ^ 

where h and w are slit height and width and F is the focal 

length. For a double multiplexed spectrometer, the throughput 

-2 2 

of MIS and DHTS are about the same, of the order 10 cm . 

For the same throughput, the HTS may have a worse system trans- 
mission because the HTS requires a dedispersing process. 

DHTS has the additional advantage that one can construct 
a one dimensional picture of the source. 

(D) Mechanical Requirements 

The HTS mask can be made self-supporting hence no beam- 
splitters or transmission optics, are required. Furthermore, 
in the MIS, construction tolerances usually involve dimensions 
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and motions that have to be malntaine'd to within fractions of 
wavelengths. For a HTS, the corresponding tolerances are frac- 
tions of a slit widths and these tolerances are normally two 
orders of magnitude more relaxed, so that this instrument will 
be more suitable for rugged applications and less costly. 

It is clear from the above comparison that the MIS has 
the advantages of highest resolution, very large multiplex 
number and free spectral range . The HTS has the mechanical 
advantages and computational advantages for large N. The HTS 
can have on the order of 10 spectral elements and a resolu- 
tion sufficient to resolve the rotational lines of many mole- 
cules. For most IR astronomical observations this will be 
sufficient. Furthermore, its potential for modifying the 
existing sacnning spectrometer at a moderate cost make this a 
very worthwhile field for further study. 



CHAPTER V 


PROGRAMMING FOR HADAMARD TRANSFORM 
SPECTRAL DATA REDUCTION 


An 8 k Computer Automation minicomputer model L.S.I. or 
model Alpha-l6 can be used to interface v?ith the output of a 
Monsanto scalar counter which digitized the output of the de- 
tector used with the Hadamard transform spectrometer. The com- 
puter processes each data point as soon as it receives it and 
when the data gathering run is completed, the computed spec- 
trum is also ready within a fraction of a second. The final 
spectrum can be displayed on a cathode ray tube for quick vi- 
sualization, or printed on paper by a teletype machine for 
more detailed analysis. Also, it can be stored on paper tape 
for future use. 

There are two inverse transformation programs: 1 x 255 

for the single entrance slit and 255 exit slit Hadamard trans- 
form spectrometer, HTS, and 15 x 255 for the 15 entrance slit 
and 255 exit slit instrument, DHTS. 

(A) HTS Program 

The inverse HTS program processes the raw data obtained 
by the combination of a single entrance slit and 255 exit slit. 
This program is in double precision format. Two areas in the 
memory are reserved by the program to store the final spectrum. 
The final spectrum can be stored in either the plus beam or the 
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minus beam area. The plus and minus beams are arbitrarily 
named. 

Data can come in at a rate of 10 data points per sec., 

5 data points per sec., or 2.5 data points per sec., depending 
on how the clock driving the spectrometer is set. -Since a com- 
plete pass has 255 points, each pass takes 25.5 sec., 51 sec., 
or 102 sec., depending on the data rate.- Each pass yields one 
spectrum. One can take as many passes as one wants until one 
is satisfied with the SNR of the spectrum. 

The whole program (Appendix B) is linked by the following 
subprograms: COMMAND, TRANSFORM, INPUT/OUTPUT, READ, CLEAR, 

PUNCH, GRAPH, DISPLAY, MATHMATICAL PACKAGE and MASK. The func- 
tion of each subprogram is described briefly in the following 
sections . 

1. COMMAND; This subprogram performs two functions. 

(a) It commands the computer to do one of the following func- 
tions: TRANSFORMATION, CLEAR, READ, PUNCH, GRAPH or DISPLAY. 

(b) If two distinct spectra are stored in two different beam 
areas, the COMMAND program can take the difference and ratio 
of the two spectra. 

2. INVERSE TRANSFORM; This is the most important subprogram 
in the whole program. It will take data points from the counter, 
transform them and enter them into either the plus or minus beam 
areas, or it will read the data points from the paper tape into 
one beam area, transform them 'and enter them into the other 
beam area. This program is called the inverse transform, since 
it inverts the transformation performed by the coding mask, and 
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yields a spectrum. The program also performs the Hadamard 
transform, i.e. transforms the spectrum back to raw data, from 
either input . 

The algorithm is based on the following idea: Let i|;j be 

the j spectral element and w^^^ be the weight of the j ele- 
ment of the mask, w^^ equals 1 for transmitted radiation 
and 0 for blocked radiation. Each measurement then has a 
value 


255 
2 S 
i=l 


'■i 


where is the random detector noise satisfying the properties 
mentioned in Section II-B. S_ is the 255 x 255 matrix, is 
the 1 data point entered into the computer. The computer’s 
job is to decode to reconstruct the original spectral values 
t|).. Therefore 

cl 




255 

2 

1=1 


,-l 

^ji ^i 


'■ -1 

where is the unbiased estimate of and S is the inverse 

of the S matrix. 

According to the relation 



one obtains S 


-1 


by keeping all +l’s in the S-matrix and re 


placing all O’s by -1, The matrix obtained in this way is the 
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2 

inverse matrix of S except for a constant factor which 

■only gives a different normalization. 

To reconstruct the spectral values one needs to add or 

subtract each measured value to different bins, according 

to whether S“^ is plus or minus. Fig. (5-1) is a flow chart 

for the 1 x 255 transform program. 

By a Hadamard transform we mean a program that transforms 

the spectrum back to its raw data . This procedure is useful 

because by inspecting the raw data display which usually appears 

quite smooth, any bad data point can be easily identified, 

and for example, replaced by the average of its two adjacent 

data points. This procedure will improve the final spectrum. 

The Hadamard transform turns out to be extremely easy. 

All one has to do is to change one statement in the inverse 

“1 

transform program. When is -1, instead of negating the 
data, one just sets it to zero. 

Data points are taken both with the exit mask moving in 
a forward and in a reverse direction. The inverse transform 
program takes care that when the exit mask moves in the for- 
ward direction, the spectrum is transformed into the plus beam 
area. When the exit mask moves in the reverse direction, the 
-final spectrum is stored in the minus area. Not adding the 

* The notation here may be a bit confusing. We use the to 
transform the raw data into an intensity spectrum. The in- 
verse transform uses S to transform the intensity spectrum 
back into raw data. 
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Figure 5-1 . 


The flow chart of the 1x255 inverse trans- 
form program. 
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forward and backward spectrum eliminates a degradation of the 
final spectrum due to any asymmetry between the data taking for 
different directions of motion of the mask. 

3« INPUT/DUTPUT ; This program links up the computer, the 
teletype and high speed reader. It consists of the following 
functions: Keyboard Input, Paper Tape Input, Output to Tele- 

type, Output Text from Buffer, Output Floating Point Number, 
Wait for Execute Signal, Command Error Exit, Carriage Return- 
Line Feed. 

READ, CLEAR, PUNCH ; This program reads the spectrum 
from the paper tape into either plus or minus beam areas for 
further manipution, or punches the spectrum out from the beam 
area; it also can clear the beam area. The speed of teletype 
for reading is 100 words per sec. Hence it takes about 12 min. 
to read the spectrum. Punching has the same rate. 

5. GRAPH : This program, plots the graph on teletype paper, 

with its numerical value in floating point format. This pro- 
cedure offers one the chance to inspect the spectrum in detail 
if it is needed. It takes about fifteen to twenty minutes to 
finish a spectrum, depending on the complexity of the spectrum. 

6. CRT ; A cathode ray display is interfaced with the out- 
put of the computer. It takes about 2 sec. to display the 
spectrum, with a factor of 5 higher resolution than the graph 
printed by the teletype on paper. One can display the spec- 
trum at the end of any pass to see how good it is. 

7. MATHMATICAL PACKAGE : This package is supplied by the 

Computer Automation library tape, with a little modification 
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from our own on its double precision part. 

^SK ; This part contains the matrix, the 255 ele- 
ments exhibited in the first mask position plus an additional 
25 ^ elements representing -the further- cycling of this mask. 

(B) DHTS 

The doubly encoded Hadamard transform program prpcesses 
the data obtained by the various combinations of fifteen en- 
trance slots and 255 exit slots. It is a single precision 
program. It can accept data at a rate of 5 data points per 
sec. and 2.5 data points per sec. only, because it takes a 
longer time to process each data point. The whole transform, 
takes about 1^^ minutes. The final spectrum consists of 15 se- 
parate spectra, representing a one-dimensional color picture 
across the spectrometer entrance aperture. Each separate spec- 
trum contains 255 spectral elements. The program can co-add 
all fifteen separate spectra yielding a sum spectrum with im- 
proved SNR. 

The- program consists of the following subprograms: 

COMMAND, TRANSFORM, DATA, INPUT/OUTPUT, READ, CLEAR, PUNCH, 
GRAPH, DISPLAY, ENTRANCE MASK, EXIT MASK. (Appendix C). ^ 

Since most of the subprograms are preforming the same function 
as their counterpart in the 1 x 255 casp, except written in 
single precision format, their description will not be repeated 
here. Only the TRANSFORM, DATA, and DISPLAY programs will be 
discussed because they are different from those in the 1 x 255 


scheme. 
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1. INVERSE TRANSFORM ; The program can accept data eigher 
from the counter or from paper tape. In order to eliminate 
any asymmetry due to the different -directions o'f' motl'on of ~the 
exit mask, the computer will accept data only when the mask is 
moving in a given direction, either forward or backvrard, de- 
pending on the operator. If one wants to save time, one can 
still choose the mode in which the computer will accept data 
in both directions of mask motion. Hence, the program provides 
six modes for operation: accepting data from the counter, in 

the forward, backward, or both directions, and accepting data 
stored on the paper tape, in the forward, backward and both 
directions. Figure (5-2) shows the flow chart for 15 x 255 
program. 

Let both the entrance and the exit masks be linear arrays 
encoded by Reed-Muller codes. Then the matrix of spatial- 
spectral elements ^ is related to the matrix of measurements 
n by 

^ JIL — ~ IL 

To obtain the spatial-spectral information about the viewed 
scene we solve this equation by premultiplying the data matrix 
by and postmultiplying by S ^ 

1 = s"\ S“^ 

Now consider the element multiplied only by 

ejemejits _pf the £irst„ column -of- .and -in t-urn it multiplies 
only the elements of the first row of S“^. To a given spectral- 
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Initialize the entrance MASK 


Store the first mitioiized entronce element into MASK 1 

I 


Initialize exit mask according to direction that the exit mosk is moving 


Store the first initiolized exit element into MSK1 
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FINISH TRANSFORM 


glKure 5-2 , The flow chart of the 15x255 inverse trans- 
form program. 
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spatial element , it therefore contributes an amount 

-1 -1 

But the elements and' all have values, 
either +1 or ~1, and the result is that each element of 
the matrix ij> receives a contribution or 

reading, 

This procedure is generally valid. Any reading will 
make additive contributions that can only have values or 

~^kA each element of the i{) matrix. 

For real time decoding we therefore need the following: 

(a) A memory that consists of bins containing the con- 
tributions to the elements acciunulated up to any given 
time t in the cycle of measurements. For a device that can 
resolve m spatial and n spectral elements, this memory re- 
ruires mn bins and of the order of mn memory words . 

(b) For each acquired reading we perform a series of 

additions of values either or one to each of the 

memory bins. But before that can be done, we need to de- 
cide on the assignment of + or - needed for a given bin. 

This is done in the following way. 


We store the sequence of + and - signs in one column of 
and in one row of S“^ and in one cycled permutation of 
each of these vectors. Let us designate these signs by their 
positions in these two vectors, as and , respectively, 
i=l,... m; j=l,... n. (Since each of these sequences' is cyclic 
it can, respectively, be brought into its kth and fi-th cycling 
position after a measurement The elements of the two 

vectors then are multiplied in all possible combinations to 
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give a matrix having mn elements. 




s"V^ 


i«l- 


m; j=l* • • • n 


Each element is either + or - depending only on whether the 
-1 -1 

signs s and are similar or dissimilar for a particular 
combination of i and j values. 

The additions or to the bins are made as 

successive elements, are computed, so, that the elements 
Z^^ need never be stored. Figure (5-“3) shows the relation of 
to a superarray containign the set of all elements that 
are constructed at various stages of the computation. 

When only a restricted number of spectral elements are 
of interest, we need to compute elements representing only 
selected J values. This might be useful, for example, if only 
certain atmospheric CO 2 absorption lines needed to be studied, 
and the spectral elements between were of lesser interest. In 
that case only s,+j-l corresponding to given j 

values need to be used, and the computing time decreases as 
p/n, where n is the total number of available spectral elements 
and p is the actual number of interest. 

—1 -1 

One starts with the inverse of the codes, s and S , for 
the entrance and exit masks stored in the computer. One stores 
509 elements of the exit mask, i.e., the 255 elements exhi- 
bited in the first mask position plus an additional 25^ ele- 
ments representing the further cycling of this mask. Simi- 
larly one stores twenty-nine elements for the entrance mask. 



2m-l 


column L 



2n - i 


Suptr-orroy of + and — tignt. 

Row m**- 1 to 2ffi- 1 on o npatitton 
of rovt I to m-l. Column* n + 1 
to ^-1 on a ropotlllonof 
column* lton~l 


PlKure 5-3. Matrices generated by the computer during the 
reduction of the spectral data. 




representing the first fifteen elements used, plus the further 
cycling of fourteen elements. 

For each reading one essentially makes use of the ma- 
trix (figure 5-3) making use of elements k to k+l4 of the stored 

entrance code and elements t to il+ 25 ^ of the stored exit code. 

“1 -1 

This matrix consists of + and - signs. When s^j^ and have 
the same sign, both being + or both being the matrix posi- 
tion ij is assigned a + sign, and the reading is added to 
the accumulatively stored value of . If the elements s~^ 
and S”^’ have dissimilar signs, a - sign is assigned to ij and 
the reading is subtracted from the stored values. This 
whole process takes vlOO msec, and is carried out while the 
succeeding intensity measurement is being made. 

In pratice, we start with the first spatial element, i=l 
and add or subtract the contributions to all the values, 
successively going from j=l to j=255. We then repeat this pro- 
cedure for i values going from 2 to 15. This whole procedure 
is carried out while the exit mask is moving from posiiton Z 
to A+1 depending on whether the exit mask is moving forward 
or back. The entire process is then repeated for the next 
reading . When the exit mask reaches its 255th posi- 

tion, l remains unchanged, but the entrance mask moves from 
the position k to k+1. 

For odd values of k, the exit mask moves in the direction 


of increasing £ values, and for even values of k, it moves to- 
ward decreasing values. In short, the exit mask moves back 
and forth as readings are taken. After the entrance mask has 
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moved. through all its fifteen positions ^ .and the total- 3825 
readings have been taken and added onto or subtracted from the 
elements, the run is completed. 

2. DATA ; Instead of processing a data point immediately 
as it comes in, this program stores the data in the memory, so 
one can display the raw data points first, correct them if there 
are any obvious bad points, and then transform them. This pro- 
gram serves the same function as the inverse transform in the 
1 X 255 system. 

3- DISPLAY The display program allows one to display the 
information in a number of different ways: 

(a) One can call for the spectrum corresponding to any 
one of the entrance slit positions and display it individually. 

(b) One can display the sum of the different spectra. 

In order to do this, one has to take into account that the 
spectrum' for a given entrance slit is displaced by one spec- 
tral position from adjacent entrance slit position. In other 
words, the wavelength for element corresponds to the wave- 
length for element because of the slightly displaced 

light paths through the spectrometer. 

(c) Finally one can display all fifteen of these spectra 
simultaneously, with the zero baseline of each spectrum verti- 
cally displaced from the next one. While this format is some- 
what crowded, _ it _do^s allow a_quick ^comparison- of—the indivi— — 
dual spectra. 
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(C) Correction Program ; 

This program is shovm in Appendix D vrritten in BASIC lan- 
guage. It corrects the error introduced by the imperfect mask. 
Instead of correcting the mask which in practice is not pos- 
sible, the program corrects the final spectrum. That is much 
easier. 

As seen in section II-D, for any spectral line I^, the 
distorted spectrum shows a line 1^=1^ (1-e), two positive blips 
with amplitude (1/2 ) (sI^/1-e ) adjacent to the line on both 
sides, and two negative blips with same amplitude, l.e. (e/2) 
(I^/l-e) at 2^, 25 elements to the left. The correction pro- 

t 

gram takes the intensity of every element, I^, multiplies it 
by e/2, adds it to the elements 24 and 25 positions to the 
left, and subtracts it from the two adjacent elements, one on 
each side of the line. The final spectrum is complete except 
for a different normalization factor. This is a linearized 
correction procedure valid only for small values of e, e<<1. 



CHAPTER VI 


ASTRONOMICAL OBSERVATION 
(A) Correction Procedure 

The correction of the spectra for telluric absorption and 
for Instrumental response is a critical procedure. The correc- 
tion is carried out by comparing the source spectrum (either 
stellar or planetary) with a lunar or solar spectrum which is 
taken on the same day at an airmass as close as possible to 
the star. The following procedures are used in the data re- 
duction: 

(1) Correct the raw star spectrum and lunar spectrum (or 
sun) for negative dip due to the Imperfect mask as described 
in section II-D. 

(2) Correct for different entrance slit width if neces- 
sary because different slit widths will give different reso- 
lution. 

( 3 ) Most of the astronomical infrared sources to be ob- 
served are weak sources, hence a positive offset is always 
added to the signal to prevent the signal becoming negative. 
(The electronics are confused by negative signals.) The cor- 
rection procedure shown previously also take out this offset. 
One can take out the offset from- the raw spectrum if one knows 
how large the offset is. Another way to correct this is by 
substracting a constant intensity from the stellar spectrxim 
imtll the ratio of the stellar spectrum over the Moon spectrum. 


90 


91 


around any large telluric absorption reigon, such as the ozone 
band, is optimally corrected. By optimal correction we mean 
that the atmospheric feature appears as neither a positive, nor 
a negative band. The Moon is a strong infrared signal and does 
not require an offset, so the result can be used as a- calibra- 
tion for base line. 

(4) Align the stellar spectrum and the lunar spectrum by the 
telluric absorption feature. The final stellar spectrum is ob- 
tained by taking the ratio of the stellar spectrum and the lu- 
nar spectrum and multiplying it by the black body temperature 
of the Moon. The lunar temperature is obtained by noting the 
phase angle of the lunar east limb where it has usually been 
observed, and extrapolating the temperature' from the value gi- 
ven by Linsky (1973)* This procedure assumes the lunar in- 
frared emissivity at 8-l4vm as unity, which is not true. Murcray 
et al’s ( 1970 ) results for the lunar enilssivity of 8-l^\im re- 
gion are shown in figure (6-1). The observation "was done from 
a balloon. The strong feature centered at 9.6vm is a result 
of telluric ozone absorption. Murcray ’s result has not been 
used in our analysis because in the region to be discussed, 8.5vm' 
l4ym, the Moon’s emissivity is about constant except for the 
ozone absorption. 

(B) Observation of a-Orlonis 

The- observations of ct-orionis were carried out with the 50" 
infrared telescope at Kitt Peak National Observatory, Arizona, 
in May 1974. The beam size is I 8 sec. x 47 sec. The Kitt Peak 
50" telescope bolometer was used with the Hadamard spectro- 

f 

meter. The dewer has a band pass of 8 - l4ym. The spec- 



Percent Emissivity 
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trometer operated in the 8-11 vm region with a resolution X/AX 
aroimd 500. The chopping frequency w&.s 10 cycles per second. 

Three runs were taken. Each run consisted of 10 scans of 
the sources. Two lunar spectra were taken on the same day 
for correction, purposes. For the lunar temperature we used 
383°K. Figure (6-2,) shows the raw spectra of a-orionis and 
the Moon. The cx-orionis spectrum is the sum of two indepen- 
dent runs. 

Figure (6-3) shows the ratio spectrum of a-orionis corrected 
for lunar temperature. Except for the region immediately aroimd 
the ozone band all the telluric absorption features are gone. 

The region between 9.35ym to 9 -Turn is unreliable because the 
ozone band has a very low transmission. 

a-orionis Is a late type super-giant with temperature 
around 3000°K. A stellar continuum corresponding to a 3000°K 
blackbody is also shown In figure(6-3), normalized to arbitrary 
units. The broad .emission feature around lOym, which is due 
to silicate emission, is clearly shovm in the spectrum. This 
lOym 'silicate’ emission feature of a-orionls has been previous- 
ly discussed by others. 

Woolf and Ney (1969) renormalized Gillett ^ aT's spectra 
( 1968 ) of a-orionis and interpreted the lOym emission as coming 
from circumstellar dust which absorbs starlight and reradiates 
at infrared wavelengths. Since the emission is far more sharply 
peaked than a black body, the wavelength dependence of the 
emission probably closely mimics the wavelength dependence of 
the opacity of the material. In the same paper, Woolf and Ney 







96 


also proposed that it is silicate grains which one is observing 

in the circiimstellar dust cloud. This suggestion has been 

strengthened by high resolution spectra of Gamown e^ ^ (1972) 

-1 —1 

with resolution X/AX=250 from 850 cm to 1100 cm . The emi- 
ssivity of silicate grains should have a second peak near 20 ym. 
This second emission feature was observed by Low and Swaray (1970) 
in narrow-band photometry of a-orionis. Another supporting 
piece of evidence for the silicate model of the a-orionis dust 
cloud comes from observations of silicon monoxide (SiO) . 

Silicon monoxide is expected to be among the most abundant 
molecules present in the atmospheres of cool stars of normal 
composition. It is also a reagent in the condensation mecha- 
nism thought to produce circumstellar silicate grains (Mass 
et al a 1970). The presence of silicon monoxide in a-orionis 
was confirmed by Cuduback e^ ^ (1971)* They observed SiO ab- 
sorption features around 4ym. 

The silicates are, expected to form from the material 
ejected by cool stars. Gilman (1969) calculated what solids 
would condense from gas of stellar composition as it moved 
away from' a star and cooled. This is critically dependent on 
the ratio of oxygen to carbon in the gas. These elements 
first combine to form carbon monoxide; the subsequent develop- 
ment depends on which of the two is left over when all the other 
has_.beejo us.e.d up_ in ihis-™way._ If-carbon predominates-, graphi-te 
will be the principal condensate, or under certain circumstances 
silicon carbide. If oxygen wins, the grains which should form 
are the silicates of calcium, magnesium, aluminiiam and iron; 



combinations of these are responsible for the 10 and 20 ym spec- 
tral features. Aluminium and calcium' silicate may be rare, be- 
cause of the low cosmic abundance of aluminium and calciiun. 
Woolf and Ney (1969) expected magnesium silicate, MgSlO^, with 
some ^ iron silicate, PeSiO^, to be most abundant. For recent 
woric on dust grains one can refer to Salpeter's paper (197^ a) 
on the theory of nucleation and dust grains in carbon-rich 
stellar atmosphere and his paper (1974 b) on formation and flow 
of dust grains in cool stellar atmospheres. 

Laboratory spectra exist for silicate absorption features 
(Day, 1974). Any fine feature of the astronomically observed 
silicate emission may be washed out by different particle sizes 
and shapes, uncertainty in temperature and mixture of composi- 
tion. Gammon et ^ (1972) examined the excess in XY Seg and 
0 Get and concluded that the type of silicates involved are 
basic rather than acidic. Day (1974) synthesized the' amorphous 
magnesium silicate and obtained an absorption band quite si- 
milar to the material causing the interstellar lOym absorption 
feature. He suggested that the existance of disordered struc- 
tures seems a more reasonable expectation than crystalline 
terrestrial-type minerals. For the momnet, the nature of the 
silicates is certainly at an unsettled stage. 

Penman (1976) had measured the middle infrared reflectivi- 
ties of five silicate minerals, and used Kramers-Konig analysis 
to obtain the optical constants of the samples. He then used 
the optical constants in Mie computations of the absorption 
properties of very small mineral grains. The final absorption 
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cross*-section spectra compared to the observed lOvm silicate 
of the source W3/IRS5. All the calculated spectra have sharper' 
features than the astronomical features due to the application 
of Mie theory. However, the hydrated silicate, Chloritite 
(hydrated Mg/Fe/Al silicate) and serpenlenite (hydrated Kg/Fe 
silicate) fits the observed astronomical positions correctly. 
They fall almost exactly at the center of the astronomical 
features. 

To the author’s knowledge only two spectra of o~orionis 
in 8~l^ym exist in the literature. Gillett et al (I 968 ) 
(figure S-'^a) obtained results in the wavelength region from 
2.8 to l^vm, "With a resolution X/ii=50* Treffers and Cohen 
( 1973 ) (figure 6-4 b) obtained a spectrum from 8-l4vm, and in 
the 20)jm region with resolution 1000. Gillett ’s and Treffers 
and Cohen’s spectra are shown in figure (6-4). Compared with 

I 

our result, all show the lOym emission feature if Gillett ’s 
black body curve is lov?ered instead of being drawn tangent to 
the observed data at 10 ym. Our spectriun shows a rather rapid 
dip at wavelengths beyond lOym while Treffers and Cohen show 
a slower nearly constant decline. Further high resolution 
observations should clear this matter up. Our own instrument 
now operates at a resolution similar to that of Treffers and 
Cohen, and if used on a telescope as large as the 120 inch 
Lick Observatory reflector that they used, sufficiently high 
signal to noise ratios should be obtained. 
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Figure 6-4 . (a) Low resolution spectrum taken by Glllett 

§ 1 ( 1969 ). Different symbols represent 
spectra taken on different nights, (b) High 
resolution spectrum taken by Treffers and 
Cohen (1973). 
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(C) Observations of Jupiter 

The observations of Jupiter were also carried out with 
the 50" infrared telescope at Kitt Peak. The beam size is 
l8" x' 47". The Kitt Peak bolometer assigned to the 50" tele- 
scope with band pass 8-l^ym was used with the Hadamard spec- 
trometer. For our Jupiter observations, the spectrometer 
operated in the 10.8-13.'4ym region with an effective resolu- 
tion x/AX around 250. The chopping frequency was 10 cycles 
per second. 

Two runs were taken. Each run consisted of twelve scans 
of the whole Jovian disk. Two lunar spectra were taken on the 
same day for correction purposes. For the lunar temperature 
we used a value of' 383°K. Figure (6-5) shows the raw spectrum 
of Jupiter and the Moon. Both the Jovian and lunar spectra 
are the sums of two Independent runs. Figure (6-6) shows the 
ratio spectrum of Jupiter corrected for lunar temperature. 

The arrows labeled by H 2 O shovr the position of telluric water 
vapor features. Most of the telluric features have been can- 
celled properly. 

Jupiter is covered by clouds which in the visible part of 
the, spectrum are seen from earth. Current models show three 
distinct cloud layers (figure 6-7, Ingersoll). The lowest 
layer is water ice, with maximum density at about 2?0°K. The 
middle cloud is solid ammonium hydrosulfide '(NH|jSH) at about 
200°K. Lewis and Prinn (1970) suggested that ultraviolet ra- 
diation from 2200 to 2700 1 is not absorbed by H 2 , He, CH|j and 
absorbed a little by Therefore, the radiation in this 
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Figure 6 - 5 . The raw spectrum of Jupiter and the Moon. 
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Figure 6-6 . The ratio spectrum of Jupiter to the Moon corrected 
for lunar temperature. 
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Figure 6-7 . 


The atmospheric profile of Jupiter for a 
solar-composition model. 
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region may reach the ammonium hydrogen cloud and photolyze hy- 
drogen sulfide there into hydrogen polysulfides (H^S ) , ele- 
mental sulfur, and ammonium polysulfides [TNHi^) 2 S^J. All of 
these species are yellow,' orange, or hrown and may explain the 
color of zones. However, Sagan and Salpeter (1976) suggested 
that even under the most optimistic assumption that every H 2 S 
photo dissociation event leads to polymeries, the implied op- 
tical depth falls short by two orders of magnitude from mat- 
ching the observed values. Moreover, pure polymeric sulfur 
fits the observed optical properties of the Jovian red chro- 
mophores only poorly (Rages and Sagan, 1977 ). The upper cloud 
is solid ammonia at around 150°K. Solid ammonia is whitish and 
probably forms the white zones on the Jovian disk. The color 
of the Great Red Spot may be due to high altitude ultraviolet 
photolysis of phosphine into ^ 2^4 amorphous red phos- 

phorus. The total depth of the cloud system is about 70 Km, and 
the pressure range probably runs from about 0.5 bar at the top 
to 4.5 bar at the cloud base. It is the cloud tops and above 
where the lOym infrared radiation originates. Our spectrum 
measures ,a ,color temperature of 135°K which is consistent with 
Ingersoll’s picture. The radiation should come from the cloud 
tops because the Jovian atmosphere at 10.5 to 13vm has appre- 
ciable opacity, as discussed in the next paragraph. 

\ 

The main constituents of the Jovian atmosphere, are hydro- 


gen molecules, helium molecules, methane and ammonia, with mi- 
nor constituents hydrogen sulfide, water, ethane and acetylene 
Table 6-1 shows their observed abundance ratio by number. 
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table 6-1 


Solar composition atmosphere (fraction by number) 


Species 

(1) 

(2) 


0,686 

0.870 

He 

0.112 

0.128 

H^O 

1.05 X 10 ”^ 

8.80 X 10 “^ 


6.30 X 10 “'*^ 

6.17 X 10"^ 

NH 3 

1.52 X 10“^ 

I.H9 X 10"^ 

H 2 S 

, 2.90 xlO"^ 

2.56 X 10"^ 


(1) Weidenschilling and Lewis (1973). 

(2) Podolak and Cameron (197^); Cameron (1973). 
The table is taken from Ingrosell. 
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The abundances of hydrogen, methane, ammonia, and helium seem 
consistent Judged from solar atomic abundances. For detailed 
information one can refer to McElroy (1973) • 

The opacity due to hydrogen molecules is caused by pre- 
ssure induced dipole absorption. The hydrogen molecule has no 
permanent dipole moment, and consequently, no permanent dipole 
spectrum. Gaseous H 2 , however, has a weak pressure - induced 
dipole spectrum which absorbs significantly over the long path 
lengths and low pressure of the Jovian atmosphere. The induced 
dipole moment results from two distinct physical process 
(Kranendonk and Kiss, 1959). The first takes place when the 
permanent quadrapole moment of one molecule Induces a dipole 
moment in another molecule by virtue of the neighbor’s pola- 
rizability. This is a long range interaction. The second 
physical process takes place when the overlap forces of the 
two adjacent molecules cause an asymmetrical distortion of 
their electronic charge clouds. The net induced dipole moment 
is modulated by the relative translational and rotational mo- 
tion of the colliding pair and this modulation produces the 
absorption of infrared radiation. The translational spectrum 
is predominant at long wavelengths with its peak at lOOym at 
100°K (Trafton and Munch, 1969). In our wavelength region 
(10.5 - 13yra) its contribution to the opacity is negligible. 

The rotational hydrogen collislonal spectrum, however, has its 
peak at 17ym and contributes a continuous opacity in our wave- 
length region (Th. Encrenaz, 1972). 

The helium molecule also has no permanent dipole moment. 
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Its opacity comes from the collision with the hydrogen mole- 
cules and resembles the collision process. The colli- 

sion is less important due to the smaller abundances of helium. 

Ammonia is an important source of opacity at lOym under 
Jovian atmospheric conditions. The lOum band of ammonia arises 
from transitions through the V 2 mode. In the mode of vibra- 
tion, the nitrogen atom oscillates vertically relative to the 
plane of the hydrogen atoms (figure 6-8a). The nitrogen atom 
is able to penetrate through the potential barrier to the other 
side of the hydrogen plane. This inverted position leads to 
inversion splitting of the levels of the ammonia molecules. 

The splitting generates both symmetric and antisymmetric energy 
levels with a given vibrational quantum number. The 10pm band 
of ammonia arises from transitions from the ground vibrational 
state to the first excited state in the V 2 mode (figure 6-8b). 
Another transition, from the first excited symmetric vibra- 
tional state to the second excited asymmetric state is also in 
10pm range, but the contribution due to this "hot band" is small 
for the low temperature in the Jovian atmosphere. 

The ammonia is clearly seen in abosrption in the spectrum. 
The centers of the bands are shown by the arrows labeled 
The ammonia absorption has been observed by different groups. 

Gillett ( 1969 ) (figure 6-9) observed Jupiter from 

2.8-l4pm with low resolution X/AX=50. Briefly, their results 
show the following: The spectrum has a depression at 3.3pm 

caused by CHij. Solar heating of the upper atmosphere via this 
band results in warming of the upper atmospheric layers. This 
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Figure 6-8 . (a) Schematic representation of the NH^ 

molecule. The components of angular momentum and 
the motion in the vg vibrational mode are also shown. 
"(b) Energy levels of 'the vg vibrational’ mode of am- 
monia. Superscripts a and s refer to the antisymme- 
tric and symmetric levels which arise due to inver- 
sion splitting. 
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proceeds imtil energy is radiated at the same rate via the 7.7viti 
band of CHij. Ammonia absorption around lOym was also detected. 
Judging from the CH^i emission, these authors were the first to 
suggest a temperature inversion on Jupiter caused by solar 
heating of the 3.3vm band of CH|j. They also showed that the 
NH^ band at lOpm is saturated, and calculated that the H 2 abun- 
dance at 12.5ym, assuming a temperature of 125°K, is 12 km-atm. 
with a oressure 'X/l/^atm. 

Aitken and Jones (1972) obtained a Jovian spectrum from 
8 - 13iim at a resolution x/AA'bli}3 (figure 6-10). The ammonia 
absorption band is again seen. They estimated that the ammo- 
nia abundance in the band is about 2.7 cm-atm. and a lapse rate 

r= at I3ym given by H=-30K, where- H is -the scale height 
dh 

-\i20km. 

The most recent published Infrared spectrum is by Lacy 
et al ( 1975 ) who used the Lick Observatory 120" telescope. 

High resolution spectra were obtained at 890 cm~^(11.2^vm) 
with X/AX=1780. Medium resolution data were observed from 
1000 cm~^ to 850 cm”^(10 'x- 12.75vm) with X/AX from 250 to 333 
(figure 6-11 a,b). The authors also calculated synthetic 
spectra, assuming that and are tthe only sources of opa- 
city. Their conclusions from comparison between observed and 
computed spectra follow: All of the prominent lines in their 
observed spectrum are saturation NH^ hands broadened to a vridth 
many’ times the pressure - broadened line width. The observed 
135°K continuum is primarily formed by the wings of the NH^ 
line. The opacity may be important if NH^ unsaturated 
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Figure 6-10. 



(a) Room temperature absorption spectrum of 
ammonia, p=0.06 atmos, w=0.6 cm atmos. 

(b) Brightness temperature; (c) Surface 
brightness of the central region of Jupiter 
from 8 to 13.5vtm. (Taken from Aitken and 
Jones ) . 




112 


near 135°K. A pressure of 0.125 atm. at 135°K is required to 
form the continuum. The minimum temperature in their synthe- 
tic model is 118±5°K while the observed minimum temperattire is 
123°K, about 5°K larger than the derived temperature due to 
incomplete resolution of the features. The lapse rate at 135 
is 7.5±2.5°K/SH. Gillett ^ aX (1969) estimated the lapse rate 
at the NH^ saturation level is 4°K/SH. A discrepancy occurs 
in the comparision of the medium resolution data between 870 
and 890 cm""^. In this region the Jovian spectrum seems to be 
depressed by about 2°K relative to the calculated curve. The 
authors suggested that it may be due to an as yet vmidentified 
minor constituent of the Jovian atmosphere. 

Our spectrum has about the same resolution as the medium 
spectra of Lacy ^ s and so the two spectra can be compared. 
The line positions match well. The vertical matching shows a 
drift towards longer wavelength. Figure (6-lla) is obtained 
by matching points A and B. The short wavelength side matches 
but the end of the long wavelength side is about 1.5 times 
higher. Also our spectrum seems to match the theoretical curve 
better at 87 O to 89 O cm“^. Figure (6-llb) is obtained by mat- 
ching A’ and B'. Then the long wavelength side matches better 
than before but the short wavelength side is different. Also 
now our spectrum matches better with Lacy ^ observed 

result at 870 - 89-0 cm“^. A conclusion about the discrepancy 
at 870 cm”^ between Lacy’s observed and calculated spectra can 
not be reached at present until there is a better way for mat- 
ching our and Lacy’s e^ ^'s spectrum. Also it is not clear 
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Figure 6-lla . Spectrum of the N and S polar regions of Jupiter at 3-i| cm”'^ 

resolution divided by the spectrum of the Moon. Data points 
are shown as solid circles, and the solid line represents 
the best fitting synthetic spectrum calculated from the 
model calculated by Lacy ^ al^ (1975HThe graph is taken 
from Lacy e;b . The dotted curve is our observed spec- 
trum by matching Lacy's spectrum at points A and B. 
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Same as 6-lla except matching our spectrum with Lacvts at 
A ' and B ' . 


Figure 6- lib. 
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whether the difference in matching of our spectrum and Lacy 
et al* 5 between long wavelength and short wavelength is real 
or not. There could be several reasons for the difference. 

It could be due to the inaccuracy of the end of the spectra, 
because the short wavelength side is the end of our spectrum 
and the long wavelength side is the end of Lacy’s spectriun; 
or it may just be due to the matching technique. More effort 
is needed to clearify this point. 

The absorption due to much less important beyond 

12p, so one may be able to use the H 2 opacity to estimate the 
lapse rate in that region. The lapse rate can be estimated 
by the equation 


Tg^i) - 


Tb(^2^ 


PH 

2 


a( X.. ) 

lu ^ 

aCx^) 


where T_(X_) is the brightness temperature at X_ 

JD X X 

is the brightness temperature at X 2 
r is the lapse rate 

a(X^) is the absorption coefficient at X^ 

a(X 2 ) is the absorption coefficient at X 2 

H is the scale height 

If one chooses X^=11.95, X^=12 . 3^ v^ith measured bright- 
ness temperature T^=133.9^, T2=133.6^, a(X^), a(X 2 ) are taken 
from Calpa and Ketebaar (1957) s one obtains a result rH=-^3*2°K 
for an H 2 opacity dominated atmosphere. Gillett e^ al (1969) 
calculated the adiabatic lapse rate (rH)ad=--42°K for an H 2 at- 
mosphere, while Aitken and Jones (1972) measured a value of 
PH=-30°K from their spectrum. Our value seems closer to the 
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value calculated by Gillett rather than to Altken’s. It is 
emphasized here that the estimate is based on the assumption 
that the opacity is the dominant opacity at wavelengths 
11.95v and 12.3^)i, which may not be true. 

Methane has a strong emission band at 1.7]i but does not 
have any band structure in our region. Methane has two im- 
portant contribution to the overall thermal structure of the 
Jovian atmosphere. The first one is that methane has an ab- 
sorption band at 3- 3pm which absorbs solar energy and rera- 
diates at 7.7pm producing an inversion temperature layer with 
a maximum temperature of 150 °K at an altitude 160 - 200 km. 
Secondly, methane is photo dissociated by ultraviolet light 
in the upper atmosphere. This results in products such as 
ethane (C^Hg), acetylene and ethylene Strobel 

(1973) estimated that column densities of C 2 H 2 , 

above the cloud top are approximately 10^^, 3 x 10^^, 

15 -2 

3 X 10 -'em respectively. C 2 Hg and C 2 H 2 were first observed 
by Ridgway (1973) using the 60” Kitt Peak solar telescope in 
the 750 - 875 cm ^ (ll.i|2 - 13. 33um)^ange with resolution 
X/AX=770, (figure 6-12a). The lines are shown in strong emi- 
ssion at the l40*^K temperature. The apparent lines are super- 
positions of many lines, each group corresponding to a sub- 
band. Ridgway calculated that the mixing ratios are N(C 2 Hg)/ 
N(H 2)74 X 10 ”^ and N(C 2 H 2 )/N(H 2)=8 x 10“^. The ratio N(C 2 Hg)/ 
N(C2H2)=50 where Strobel predicts about 200, Combes ^ aT 
(1974) ’s (figure 6-12b) observation confirms the presence of 
very strong emission lines of C 2 H 2 and 


The abundance of 
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Figure 6-12 . (a) Thermal emission spectrum of Jupiter cor- 

rected for absorption in the earth ' s .atmosphere observed by 
Ridgway (1973). The dashed line is the predicted form of 
the H 2 continuum. (b) The ratio of the Jovian spectrum to 

the atmospheric absorption spectrum observed by Combes ^ 
at (1974). The solid and dashed lines are the blackbody 
curves at 135°K and 120°K respectively. 
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ethane estimated from Ridgway*s spectra depends strongly on the 
distribution of gas temperature, whicH is not well-determinded. 

If the temperature in the mesospheric inversion layer turns out 

to have a maximum value of 150 °K, the observations indicate 'u 

—2 -1 

(Ridgway 197^) 2 x 10 gm cm of ethane in this high tem- 

perature region. Sagan and Salpeter (1976) estimate the column 

_-5 ^2 

density of ethane molecules to be 3 x 10 gm cm by assuming 
that ethane is produced by photolysis' of methane by solar ul- 
traviolet photon and destroyed mainly by eddy diffusion into 
the troposphere, followed by pyrolysis in deeper, hotter layers. 
This would be in very serious conflict with the observations, es- 
pecially since only a small fraction of the theoretical column 
density refer to the hotter inversion region. 

The ethane emission band is also shown in our spectrum 
(figure 6-6). In the figure the indicated emission line posi- 
tion was extrapolated from Ridgway ’s spectrum, and the laboratory 
observed position by Smith (19^9) are also shown for comparison. 
Our positions agree with Ridgway ’s reasonable well, while Smith's 
seem displaced from ours by 0 . Olym, possibly due to a uncertain- 
ty in position calibration. Only a portion of the <^ 2^2 
can be seen in our spectral coverage. The abundance of is 

not estimated here because the absolute amplitude of our spec- 
trum is not well calibrated. 

Our spec.trum contains both ammonia and ethane features 
while other observers have not shown both. This will be use- 
ful because one can compute the synthetic spectra including 
both. Ammonia will provide us with information about the top 
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of the cloud layer while ethane provides us with the informa- 
tion about the inversion layer. A synthetic spectrum including 
both ethane and acetylene would be interesting because the jn- 
clusion of these new gases would affect the models, especially 
around the inversion. Acetylene would appear^ around 13vm. 

The absorption of solar radiation by at 3.3vni used to be 
thought to be radiated solely by CH|^ at 7*7vm. The 7.7vni emi- 
ssion intensity is a critical test of a temperature inversion 
model and the emission intensity calculated by Wallace et ^ 
(197^) is within 25/5 of the value observed by Gillett e;^ ^ 
( 1969 ). If ethane and acetylene do contribute to emission in 
the thermal infrared, there must be some additional source of 
solar absorption in order to produce the observed inversion tem- 
perature. Additional absorption at this altitude, perhaps due 
to particles, is suggested by the low ultraviolet albedo of 
Jupiter in the wavelength region 2100 to 36 OO A (Wallace, 

Caldwell and Savage, 1972). 

Terrile and W^estphal (1976) had imaged Jupiter at high 
spatial resolution at 8 - l4ym. All images reveal a belt and 
zone structure similar to visible photographs. In the 8 - l^yra 
broad-band data, belts appear to be about 2°K hotter than the 
zone. The lowest belt-zone contrast is found in the hydrogen 
opacity dominated region at 12.5ym, while images at 9.5ym have 
the greatest contrast. This is consistent with the dynamic 
picture that zones are rising columns of air and belts are 
sinking columns of air. Ammonia gas being carried upward in 
zones will freeze' out and form a thick cloud on top of the zones. 
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giving a low infrared temperature to the zones, and the crys- 
talized NH^ particle will be carried down to the deep atraos- 
phers in the belt where they will get sublimated. One can look 
deeper into clouds in the belt because of the lack of the am- 
monia cloud on top of it and therefore see a higher infrared 
temperature. Although there are a number of high resolution 
spectra of good spectra for methane. The observation of methane 
will be interesting not only because it provides us with in- 
formation about the inversion layer, it is also useful to find 
out the temperature profile. If the temperature profile of 
the Jovian atmosphere is known, one can use it to find the am- 
monia abundance profile.. Ammonia itself is not a very good 
tool for probing the temperature profile because it has a low 
vapor pressure and its variation is very sensitive to tempera- 
ture changes. Observations of methane at 3*3ym and 7.7ym, 
should be able to accomplish this. 

The Hadamard transform spectrometer described in this 
thesis would be able to make these observations, with small 
modifications that would permit observations to be made at these 
wavelength. In addition, observations should be undertaken at 
8.0 to 9.5 microns where neither ammonia nor methane have strong 
absorption features. At these wavelengths one would be ob- 
serving the clouds. In this 8.0 to 9.5 micron region our in- 
strument should be able to image the bands and zones of Jupiter, 
to probe for spectral differences and cloud features. Such ob- 
servations should increase our understanding of Jupiter’s cloud 


structure. 
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(D) Observations of Mercury 

The observations of Mercury were' made with the newly built 
Cornell 25" telescope at Mount Pleasant, Ithaca, New York. 

The telescope has a focal ratio f/13.5 • The spectrometer’s 
acceptance beam size is "tTb" x "^8" . The dewar described in 
section III-B was used with the spectrometer. The spectrometer 
operated in the 10.5<\/13vm region with a resolution A/AX around 
300. The chopping frequency was 10 cycles per second. 

The observational procedure was carried out a little dif- 
ferently from the observations of a-orionis and Jupiter. 

First, Sun spectra were used for correction spectra rather than 
lunar spectra. There are no known molecular lines in this re- 
gion. Its temperature at ll.lOym is 5030°K(Saildy & Goody). 

The 'observations were made on August 3, 1976. At that time the 
Sun was about an hour away from Mercury, and was observed 
through roughly the same air mass as Mercury, Secondly, Mer- 
cury is so faint in broad day light that we were imable to see 
it in visible light. The way to find Mercury was the following: 
We pointed the telescope in the correct region and scanned for 
the infrared signal. The signal is so strong that one can see 
it go off scale on' the synchronous demodulator. The pointing 
accuracy of the telescope is 6 sec. of time in right ascension 
and 25 sec. of arc in declination. Thirdly, since we could not 
see Mercury visually for tracking, we adapted a different me- 
thod for tracking Mercury. Since our computer programming is 

r 

set up in such a way that the data taken when the mask is 
moving forward and moving backward are stored in different 
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areas, only took data when the mask was moving forward. 

VJhen the mask was moving backward we maximized the signal to 
assure correct pointing and waited for the next forward pass. 
Any noise introduced when moving the telescope for maximizing 
the signal would have gone into the backward-pass data bins 
and those data points were thrown away anyway. Since each pass 
takes 51 sec. only, Mercury remained at essentially the same 
position during the forward data taking pass. 

Two runs of Mercury and two runs of the Sun were taken. 
Each run consisted of ten scans of the sources. Figure (6-13) 
shows the raw spectra of Mercury and the Sun. Since Ithaca 
has a lot of moisture in the air during the summer time, the 
correction for atmospheric features is more difficult than at 
Kitt Peak and is done in a different way. A constant was added 
to the Mercury spectra such that atmospheric features in the 
Mercury /Sun ratio spectrum were minimised. This step ensures 
that the atmospheric features are largely corrected. The Mer- 
cury spectrum which has a constant added to it was then multi- 
plied by another constant to make its amplitude as close to 
that of the solar spectrum as possible. The solar spectrum 
was then subtracted from the modified Mercury spectrum. The 
multiplication of the Mercury spectrum by a constant assured 
that atmospheric features in the two spectra had similar am- 
plitudes before the subtraction step. This difference spec- 
trum was then added to a ’’perfect” solar spectrum which is 
calculated according to the, blackbody function appropriate to 
the solar temperature. What one gets from these procedures is: 
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Final Mercury Spectrvmi 

= observed corrected Mercury spectrum ~ observed solar 
spectrum + perfect solar spectrum 
= (’'perfect'* Mercury spectrum + noise in the Mercury spec- 
trum) - ("perfect" solar spectrum + noise in the solar 
spectrum) + "perfect" solar spectrum 
= "perfect" Mercury spectrum + (noise in the Mercury spec- 
trum - noise in the solar spectrum) 

Any systematic noise such as emission and absorption due to the 
sky or to the telescope will be subtracted away. The advantage 
of this method is that the final spectrum is obtained through 
subtraction rather than by division. Division, in the low sig- 
nal portion of the spectrum, produces deceptive high noise 
spikes in the ratio spectrum. The method we have used tends 
to eliminate these. 

Since this will be the first high resolution Mercury spec- 
trum obtained, we will calculate Mercury's disk integrated in- 
frared temperature and compare this temperature with the observed 
one. Morrison and Sagan (1967) had calculated the infrared 
brightness temperature of the center-of-disk as a function of 
phase angle and heliocentric longitude, but there is no disk 
integrated infrared temperature available in the literature. 

In the following we will discuss the factors that may affect 
the brightness temperature, and then- present a method of cal- 
culating the disk integrated infrared brightness and compare it 
with our observation. A good review of thermophysics of Mer- 
cury is given by Morrison (1970). 
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In 1965 , Pettinglll and Dyce (1965) used radar to discover 
that Mercury has a rotation period of ’ 59 days, two t-hirds of 
the orbital period, instead of an 88 day synchronous rotation 
around the Sun. This implies that Mercury has a solar day, 
on the planet, about 176 terrestrial days long, equal to two 
orbital revolutions in three rotations. Mercury’s non-synchro- 
nous rotation leads to time-dependent thermal emission of the 
planet due to the diurnal variation of the insolation. This 
diurnal variatipn would not happen for a synchronously rotating 
Mercury. The diurnal variation changes the brightness tempera- 
ture as a function of both phase angle and heliocentric longi- 
tude. It also allows a measurement of the thermal properties 
of Mercury's surface. 

Because of the high eccentricity of Mercury’s orbit, 
(e=0.2), the diurnal cycle of Insolation is markedly different 
from longitude to longitude, and can differ by a factor of 2.5. 
The eccentricity enters in two ways. First, the variation in 
distance from the Sun produces a solar constant that varies by 
more than a factor of 2 from perihelion to aphelion. Second, 
the changing orbital angular velocity causes the apparent speed 
of the Sun across the sky to vary; near perihelion the angular 
velocity of revolution actually slightly exceeds the angular 
velocity of rotation, and the apparent planetocentric solar mo- 
tion is retrograde (figure 6-li{, Soter and Ulrichs, I 967 ). 

The two effects of the eccentricity reinforce one another, with 
the larger flux coming at a time when the angular rate of the 
Sun across the sky is largest. The two longitudes (180° apart) 



Figure 6 - 1 ^. 



Diurnal path of the Sun about Mercury, drawn to scale. 
The relative positions of the Sun are marked at 11 day 
intervals with the planet held as a fixed reference. 
Planeto-graphic longitude are Indicated for Mercury. 
(Taken from Soter and Ulrichs, I967). 
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that see the Sun overhead at perihelion receive more than two 
and a half as much energy per period as the longitude 90° away, 
where the Sun is always small and rapidly moving while near the 
zenith. 

Besides the insolation geometry, a possible atmosphere on 
Mercury will also affect the thermal emission. CO^ is a major 
product for a possible secondary atmosphere, furthermore, since 
CO 2 could have been photodissociated and reduced to CO by pre- 
ferential loss of oxygen, Fink ^ ^ (1973) had set up a search 

for a possible Mercury atmosphere of CO^ and CO. They set up 

-it 

an upper limit 1.0 x 10 mb surface pressure for CO 2 and 2.0 x 
-5 

10 -^mb for CO. Mariner 10 ’s results also suggest that Mercury 
has no atmosphere although it may have a thin layer of He and 
other inert gas trapped by Mercury's magnetic field. 

The optical observations of Mercury show that the Integral 
spectral reflectivity of Mercury is quite similar to that for 
the integral moon (McCord and Adams, 1972). The Bond albedo 
of Mercury (de Vancouleurs, 196^}) is 0.058. 

In our model calculation we assume the following things. 
This model has been described by Murdock (197^^): 

1. The emission from the dark side at the phase angle we 

observed is negligible. On the day we made our observations, 

/ 

August 3, 1976 , the illuminated portion- amounted to 0.973 of 
the total disk and the dark portion was O.OI 7 . The dark side 
temperature is about 110°K, which at 10pm has a flux 6 . 718 I x 
10~®watt-cm'*^-u*‘^-sr'"“. The flux for 500°K at 10pm is 7.1007 x 

_'5 _2 —1 —1 

10 watt-cm -p -sr , so the contribution from the dark side 
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is negligible. 

2. At infrared wavelengths, the radiation that reaches the 
observer originates very near the surface, so the infrared tem- 
perature is assumed equal to the insolation temperature. Soter 
and Ultichs’ (196?) results show that the day time temperature 
is independent of the thermal properties of the surface material 
and determined largely by the insolation temperature. 

3. The infrared emisslvity we assumed was 0.9> which is the 
lunar value. We choose this value since Mercury’s surface may 
be similar to the lunar surface. 

In'" figure (6-15) we choose two coordinate systems on Mer- 
cury surface. The unprimed system is the ”sPlar system” with 
the z-axis pointing towards the Sun. A is the subsolar point. 
The hemisphere above the plane BCD facing the Sun is the illu- 
minated part. The primed system is the ’’earth system’ with the 
z’-axis pointing towards the earth. A’ is the subearth point. 
The hemisphere above the plane B’C’D' facing the earth is the 
portion that is being seen from earth. The two systems are 
different by an angle a with the x-axis as the common axis . 

If the surface is in equilibrium with sunlight and cannot 
conduct heat away, the subsolar point temperature is: 


S(l-A) 1/4 


( 6 - 1 ) 


6 ••1 

where S is solar constant at earth, equal to 1.360x10 erg cm s 

A is the Bond albedo for Mercury, assigned to be 0.058. 

'—5 

a is the stefan-Boltzmann constant equal to 5.67x10' 

-2 -1 
erg cm s 
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Figure 5-15. 



!Two coordinate systems on the surface of 
Mercury. The unprimed system is the "solar 
system" with the 2-axis pointing towards the 
Sun. A is the subsolar point. The primed 
system is the "earth system" with the Z'-axis 
pointing towards the earth. A’ is the sub- 
earth point. 
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e is the infrared emisslvlty at lO^m from the surface, 
assumed to be 0.9 

R is the distance between Mercury and the Sun in astro- 
nomical units 

Mercury’s subsolar point temperature varies with distance from 
—1 /2 

the Sun as R and therefore is a function of heliocentric 

longitude due to the eccentricity we discuss above. 

In the unprimed system the temperature distribution on 
the surface will be concentric isothermal bands around the sub- 
solar point. The temperature of a band at colatitude e is given 
hy 

T(e) = T^ cos^'^^e (6-2) 

The total intensity at any wavelength region will be composed 
of the contributions at that wavelength from many isothermal 
regions each with its own apparent area. 

Since one is interested in the flux coming to the earth 
one will have. 


dF- 

= 1^(9’, (j.’) 

dA’ 

(6-3) 



r^sinS ’ d6 ' d4> ' 

(6-il) 



cose ’r^sinQ ’ d0 ’ d4> ’ 

(6-5) 


where dP^ is the flux at wavelength, X coining from the area dA’ 
is the intensity at (6’,ij)’) 
dA’ is the differential area on Mercury 
r is the radius of Mercury 
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= 1^(0 ' ,4 * )/cos6 ’ is the component of flux 

that radiates toward the 
earth. 

is a complipated function of (6’, 4’’) hacause the 
temperature distribution is a complicated function of (9 
However, one can convert the system to the ”sun coordinate" 
where 5^(9,(})) is a simple function. 

Since ' dA’ = dA (6-6) 

cose’ = Gosecosa - sinBsin^sina (6-7) 

equation(6-5) becomes 


dP^ = [^cosQcosa - sinBsintfisinc^ r^sin6d6d(j> 



where Ci= 1.1909 x lo\att cm”^y“^sr”^ 

C2= 1.^388 X lO^uK 
and T(e) is given by (6-2) 
so 

= r^/de/d(}) H^(9 ) (sinecosecosa - sin^esin^sina) 

(6-9) 

The integral (6-9) can be separated into two parts, for 9<0_ 
where ^ = I" " “ (6-10) 

is the limit of the cap shared by both the Sun and the earth, 
i.e. the limit of integration over de can 4> go from 0 to 2ir. 
For e>0, the limit of the integral over d is constrained by 
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the physical condition that some part on that is illuminated 
by the Sun can not be seen from the earth. 

So (6~9) becomes 


P 


X 



2^ 

d'Qf d <}> 
0 


'X 




ir 


de/ 

e 


^ (e) 

2 

<l>l( 0) 


( e) (sinecos 0COS a - sin^a sin(|isina) 


d(j) E^( e) (sin 6COS 6COS a - sin^esinrjisina) 


( 6 - 11 ) 


To find if>i(0) and <(.2(0) one notices that ^ is given when 6’=90°. 


From 

(6-7) 




COS0' = 

cos0cosa - sinesin^sina 



e» 

^ =4 sinij) = eotacoto 

(6-12) 

and 


sin~^(cotacot e) 

(6-13) 


Since (fi will also be symmetric about the y-axis, equation (6-11) 
can be rewritten as; 


IT 

•2 2it 

~ (0 ) (sin0cos0cosa - sin^esin4>sino) 


^ 2Tr+sin”^(cotacot0) 

+ / d0/ E (0)(sin0cos0cosa - 

_ ir-sin (cotacots) ^ 

2 sin^esin^ysina) 

(6-li|) 


After, evaluating the integral ..one obtains the followingi 


2nr^ cosa 


H,(0) 


sin^-0 


2 “^+ irr^cosa 


5,(6)^ 


IO|=l 
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TT^ 

+ de sin~*^(cotacot 0 ) sinec'osecosa 2^(0) 



+ I S ( e ) (cose /sin 


— ^ — — 

a- cos 0 + sin asln 


cose 

sina 


where 


=,(e) 


sin^e 





(6-15) 


is the mean of ( H^'(e) in the interval of 6' from 0 to 

-a . The same meaning applies to the third term of (6-15) • 
As a check of equation (6-15) s if ct=0, that means when 
subsolar point and subearth point coinside, ^let H( 9)= constant 
evaluating (6-15) gives: 

F = irr^H 

If a= , that means the subsolar point and subearth point are 
90° apart. With h(’B) assumed to be constant, equation (6-15) 
gives 



which is as one expects since one is seeing half of Mercury. 

Equation (6-15) can be readily integrated on a computer. 

It is applied to our case with the following physical parameters: 
Date: August 3, 1976 . 

Phase Angle: 53° 


roj=< ro(=i 
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Radius vector: 0.^l4 A,U. 

Orbital longitude; 198.09*^ 

Mercury perihelion point: 0.3075 A.U. 

Subsolar temperature at perihelion point: 700°K 

Subsolar temperature at cc=53°: 603°K 
Equation (6-15) was computed on a LSI mini computer at each 
wavelength, from 10.6 to 13.2pm. The program and the result are 
shown in the Appendix D. The integral was divided into twenty 
steps. E(0) was evaluated from equation (6-8) 


^ X5 ^C2/kT(e)_^ 

1 /u 

where T(e) = T^cos 6 

The calculated spectrum is a measure of color temperature 
and is used to compare with the observed spectrum. Figure (6-l6) 
shows the final Mercury spectrum conected for solar temperature, 
with a number of blackbody slopes shown to match. The calculated 
spectrum (cross) matches the blackbody temperature 525 °K, which 
also matches the observed spectrum. We concluded that the best 
fit lies in the 500°K region. Murdock (197^) measured a effec- 
tive brightness temperature at 10.8pm at the same phase angle 
to be around 650°K. Our results disagree with his results. 



Relative Flux 


in 

1 — I 



11.0 11 5 12.0 ' 12.5 13.0 


Figure 6~l6 . The final Mercury spectrum, corrected for solar tem- 
perature, with a number of blackbody slopes shown to 
match . 


APPENDIX A 


ESTIMTE 0 ? CODING ERROH FOR FOURIER 
TRANSFORM SPECTROMETRY 

Let 5 be the path difference in a two beam interferometer. 
?(v) is the power at wavelength v» (i.e. the spectral density 
function) v here is taken to.be 1 /x. S(g) is the power received 
for path difference 5. Then (p.96 Stewart) 

S(c) = /^PCv) cos^(2it5v) dv A~1 

0 

= 1/2P„ t l/2/”P(v) C0s(l|,rSv) dv A-2 

o 0 

The reciprocal Fourier property is (Morse and Feshback P.^ 5 ^) 

that if 

F(s) = /"cos(5v )f (v)dv 

IT 0 

Then f(v) = /“cos (£v )F( 0<^5 

0 

which implies, neglecting the constant term, that 

P(v) = 16 /“cos(^tt 5 v )s (5 )dg A-3 

0 

Suppose we take measurement at (N+ 1 ) equally seperated 
steps in the variable 5 . Let the step length be t , then 
C = nx 

In general x is chosen such that 
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T 


1 


iKv 


max 


V . ) 

nun' 


A-4 


which, 

cycle 


if V . <<v , 

min max ^ 

(Stewart, p.303) 


effectively implies sampling twice per 


One can now write the integral for P(v) as 


N 

P(v) = l6z TS(nx) cos(4irvnT) A-5 

n=0 


Now consider frequency v= + m6 


where 6 


V — \) J 

max min 
N 


a-6 


and m is an integer. m=0,l**»* N 

The 5P(v) is the power in the resolved spectral band width 
6 about frequency v 


N 

6P(v) = 16t5 1 S(nT) cos(4irvnT) A-7 

n=0 


but t5 = 

^^'’max ~ ’^min^ 

= -L. 

4N5 


4N 


A-8 
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Therefore 

<5P(v) 


H 

N 


N 

Z 

n=0 


S(nT) 


cos(^TtvnT) 


A-9 


Writing this out in matrix fornij with ® ) 
and with have 




7 

h 

> 

' — ' 

j 


r 

1 COS6i COSN01 

S(0) 

« 


• 

1} 

1 COS02 C0SN92 

• 

s(l) 

« 

• 

5 5 

• 

• 

• 

P(v -6) 

max ' 

P(v ) 

max 

N 

« 

• 

1 COSSjj^^^**' COsN9j^__^ 

1 COS0 • • ♦ • cosN0„ 
n N 

S(N) 
~ - 


or ir = W~^S A-10 

From (2-8) 


lo 


(1 + co5^e. + 

n2 i 


+ oos^Ne^ ) 


1/2 


I ( Z cosine 
n=Q 


I 


n 

2 


/I 

^ W 


for N>>! 


A-11 


The SNR improvement is the reciprocal of this quantity. 
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APPENDIX B 
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IN VS 

0038 

0019 

F9 00 


JST 

ERR 

0039 

OOlA 

0110 

MON 

EAR 


0040 

0013 

0210 


CAR 


0041 

OOlC 

9AAG 


STA 

TETR 

0042 

OOlD 

0110 


EAR 


0043 

ODIE 

F20D 


JMP 

NEXT 

0044 

OOIF 

0110 

TAPE 

EAR 


0045 

0020 

9A9C 


STA 

TETR 

0046 

0021 

F202 


JMP 

T4 

004-7 

0022 

C60I 

IN VS 

LAP 

1 

00^8' 

0023 

9A99 


STA 

TETR 

0049 

0024 

0108 

T4 

EXR 


0050 

0025 

F9 00 


JST 

IKB 

0051 

0026 

COAB 


CAI 

*+ » 

0052 

0027 

0408 


CXR 


0053 

0025 

COAD 


CAI 


0054 

0029 

0105 


EXR 




0055 

002A 

EASE 


STX 

0056 

002B 

0110 


2AR 

0057 

0020 

E29 2 

NEXT 

LDX 

005S 

002D 

0128 


IXR 

0059 

002£ 

0210 


CAR 

0060 

002F 

9A7F 


STA 

006 1 

0030 ' 

00?? 


CAI 

0062 

0031 

F203 


JM? 

0063 

0032 

0030 


TXA 

0064 

0033 

8A7A 


ADD 

0065 

0034 

0048 


TAX 

0066 

0035 

EA7 6 

F0UD 

STX 

0067 

0036 

C7F5 


LAM 

0068 

0037 

9A79 


STA 

0069 

0038 

07 FF 


LAM 

007 0 

0039 

9A76 


STA 

0071 

003A 

FABl 


JST 

0072 


* IN PUT 

L00P 

007 3 

003B 

FABB 

ILP 

JST 

0074 

0030 

9A60 

* 

STA 

0075 

003D 

EA6C 


STX 

0076 

003E 

07FF 


LA-I 

007 7 

003? 

9A72 

• 

STA 

0078 

0040. 

B26B 


LDA 

0079 

004 1 

9A6B 


STA 

008 0 

0042 

327 1 


LDA 

008 1 

0043 

9A0F 


STA 

008 2 



*TRAI'IS?0R:<I 

008 3 

0044 

E265 

TLP 

LDX 

0084 

0045 

Olio 


2AR 

003 5 

0046 

D366 


Q«3 5 

003 6 

0047 

F203 


JMP 

008 7 

0048 

B274 


LDA 

0088 

0049 

218 3 


JAL 

0089 

004A 

0110 


EAR 

009 0 

004B 

0108 


ZXR 

009 1 

0040 

F204 


JMP 

009 2 

004D 

B25B 


LDA 

009 3 

004E 

F9 00 


JST 

009 4 

004? 

F-201 


JMP 

009 5 

0050 

B253 

TIA 

LDA 

009 6 

0051 

DA53 

T1 

IMS 

009 7 

0052 

F900 


JST 

0098 

0053 

0000 

IBP 

DATA 

0099 

0054 

0202 


a;<i 

0100 

0055 

EE02 


STX 

0101 

0056 

DA53 


IMS 

0102 

0057 

F613 


JMP 

0103 

0058 

E2S3 


LESC 

0104 

0059 

0123 


IXR 

0105 

005A' 

B254 


LDA 

0106 

005B 

0000 


CAI 

0107 

0050 

0302 


SXI 

0103 

005D 

EA4E 


STX 


TRH'l, 


MASK 

CHA^IGE DIRECTION 

DIR 
; FF 
F0UD 


H255 


MSKl 

245' 

DCT2 

255 

START data C0UNT 

DOT 

TURI'J0N 

WAIT FOR LIGHT 

M0NS 

INPUT 

CNT 

aSAVE 

CMT+ 1 

255 

SOT 

SPECTRUM C0UOT 

MSK 1 
nSK2 

RESET MASK POINTER 

13 

GET BUFFER POINTER 

IBP 

C0LU-IN T0 

INPUT BUFFER 

□'JT+ 1 


*MSK2 
TIA 
TETR 
S+ 4 



T1 

CNT 


DR^i 
T1 
CNT 
MSK 2 


DPACC 

0 

ADD 

2 

IBP 

BUMP P0IiNTER 

SOT 

D0NE? 

TT-P 

N0 

MSK 1 

YES.. MOVE C0LUI-1N 

DIR 

O'- 

+ IIF FWD 

2 

MSKl 

-IIF REV. 



0109 

005E 

DAS 2 


IMS 

DOTS 


0110 

005F 

F201 


JMP 

£+2 


0111 

0060 

4007 


SEL 

24r 7 


0112 

0061 

DA4E 


IMS 

DOT 

END 0F ROW? 

0113 

0062 

F627 


JMP 

ILP 

N0 

01 M 




0F R0V/ 

ST0P 

TEST 

0115 

0063 

B259 


LDA 

TETR 


0116 

0064 

208 1 


JAl'l 

£+2 


0117 

0065 

F209 


JMP 

STOP 


01 18 

0066 

3403 


JSS 

ST0P 

SS D0WN? 

0119 

0067 

43 07 


SSN 

: 07 

N0^ ST0P D0WI>I? 

0120 

0068 

F206 


JMP 

STOP 

YES 

0121 

0069 

FA56 


JST 

ACC 

N0^ SAVE DATA 

0122 

006A 

E24D 


hISK 

N0H-I 

• 

0123 

0063 

0110 


ZAR 



0124 

0060 

FiOO 


JMP 

D1 SP 


0125 

006D 

B241 

RTTR 

LDA 

DIR 


0126 

006S 

F642 


JMP 

NEXT 


0127 

006F 

E243 

STOP 

LDK 

NI 


01 25 

0070 

0601' 


LAP 

I ‘ 


0129 

0071 

FIOO 


JMP 

DI SP 


0130 

0072 

C6BF 

RT0P 

LAP 

‘ * 


0131 

0073 

F9 00 


JST 

0TT 


0132 

0074 

01 IF 


DATA 

MSG 


0133 

0075 

F900 


JST 

IKB 


0134 

0076 

COCE 


CAI 

’N ' 

AB0RT LAST PASS? 

0135 

0077 

F203 


JMP 


N0 

0136 

0078 

C0D9 


CAI 

•Y* 

YES 

0137 

0079 

-F203 


JMP 

YES 

YES 

0133 

007A 

F608 


JMP 

RT0P 

WR0NG H'lTRY 

0139 

007B 

FA44 

N0 

JST 

ACC 

SAVE DATA 

0140 

0070 

F203 


JMP 

Y2 


0141 

007D 

06AC 

YES 

LAP 

’ 

PRINT AB0RT 

0142 

007E 

F9 00 


JST 

on: 

• 

0143 

0077 

0128 


DATA 

ABT 


0144 

OOBO 

B234 

Y2 

LDA 

CT 

SET COU-'ITER 

0145 

003 1 

9A30 


STA 

SOT 


0146 

008 2 

0350 


ARP 


RESET SUBSCRIPTS 

0147 

003 3 

9900 


STA 

XB 

. 

0148 

008 4 

9900 


STA 

XC 


0149 

003 5 

E230 


LDJC 

NPEM 


0150 



*ADD 

T0 BEAM ARRAY 

0151 

003 6 

EE33 

ALP 

STX 

IBP 


0152 

003 7 

B400 


LDA 

§0 

GET DATA 

0153 

0038 

E401 


LD>C 

§1 


0154 

0039 

F9 00 


JST 

DPFLT 

FLOAT IT 

0155 

008A 

9A21 


STA 

MS1^ 1 

SAVE IT 

Oi'56 

OOS.B 

EA2-1 


STX 

M-SK2' 


01*57 

003 0 

F9 00 


JST 

FAD 

ADD T0 BEAM 

0158 

008 D 

OOAC 


DATA M SK 1 


0159 

003 2 

0000 

BEAM 

RES 

2* O' 


0160 

009 0 

D9 00 


IMS 

XB 

BUMP SUBSCRIPTS 

0161 

009 1 

D900 


IMS 

XC 

* 

0162 

009 2 

E63F 


LDX 

IBP 




1^6 


0163 

009 3 

0202 


AXI 

2 

0164} 

0094 

DAID 


IMS 

SCT 

0165 

009 5 

F60F 


JMP 

ALP 

0166 



* PRINT 

PASS 

C0UNT 

0167 

009 6 

F900 


JST 

CRLF 

0168 

009 7 

0702 


LAM 

2 

0169 

0098 

9A23 


STA 

C0N1 

0170 

0099 

E21C 


LDK 

NPE4 

0171 

009A 

0110 

MK 

EAR 


0172 

009 B 

E401 


LEK 

©1 

0173 

009 C 

F9 00 


JST 

DPFLT 

0174 

009 D 

9A0E 


STA 

MSKl 

0175 

009 E 

EAOE 


STX 

MSK2 

0176 

009F 

F9 00 


JST 

0FPA 

0177 

OOAO 

OOAC 


DATA 

MSKl 

0178 

OOAl 

Olio 


ZAR 

• 

0179 

00A2 

F900 


JST 

0TL 

018 0 

OOAO 

0130 


DATA 

RUNS 

018 1 

00A4 

E212 


LDX 

NMHl 

0182 

00A5 

DAl6' 


IMS 

C0N 1 

0183 

00A6 

F60C 


JMP 

MK ■ 

0184 

00A7 

F7A7 

T2 

RTN 

TRANS 

0185 



*DATA 

ST0RAGE 

018 6 

00 AS 

FAOO 

Z CT 

DATA 

- 1536 

0187 

00A9 

0000 

el'll 

DATA 

o; 0 


OOAA 

0000 




0188 

OOAB 

1400. 

NT 

DATA 

: 1400 

0159 

OOAC 

0000 

N SK i 

DATA 

O' 

019 0 

OOAD 

0000 

MSK2 

DATA 

0 

019 1 

OOAE 

OOFF 

H255 

DATA 

255 

019 2 

OOAF 

0000 

DIR 

DATA 

0 

019 3 

0030 

0000 

DOT 

DATA 

0 

0194 

OOBl 

0000 

DCT2 

DATA 

0 

0195 

0032 

0000 

SOT 

DATA 

0 

0196 

OOB3 

1800 

N1 

DATA 

: 18 00 

0197 

00B4 

1802 

13 

DATA 

: 18 02 

0198 

0035 

FEOO 

CT 

DATA 

-512 

0199 

0036 

1400 

N PK-I 

DATA 

V1400 

0200 

0037 

1600 

NM S-I 

DATA 

; 1600 

0201 

00B8 

0000 

N0H>1 

DATA 

0* 

0202 

0039 

0000 

trhm 

DATA 

0 

0203 

OOBA 

9 002 

FLS 

DATA 

: 9002 

0204 

003B 

9 202 

MINS 

DATA 

:9202 

0205 

ODBC 

0000 

C0N 1 

DATA 

0 

0206 

003D 

0000 

TETR 

DATA 

0 

0207 

OOBE 

1000 

BR'JT 

DATA 

I 1000 

0208 

OOBr 


MASK 

REF 


0209 



* 



02*10 



INPUT ARRAY : 

^0211 



* 



0212 

0000 

08 00 

ACC 

ENT 


0213 

OOCl 

3604 


LDA 

TETR 

0214 

00C2 

208 3 


JAM 

S+4 

0215 

0003 

360A 


LDA 

trem 


D0NE? 

N0 
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0216 

00C4 

0210 


CAR 



0217 

00C5 

F201 


JMP 

S+2 


0218 

D0C6 

3617 


LDA 

DIR 


0219 

00C7 

COOO 


CAI 

0 


0220 

0008 

F202 


JM-P 

S+3 


0221 

0009 

E613 


LDK 

NPai 


0222 

OOCA 

F201 


JMP 

S+2 


0223 

OOCB 

E6I4 


LEK 

NMBl 


0224 

0000 

EEl4 


S7K 

N0BI 


0225 

OOCD 

DCOl 


IMS 

el 

3l?MP C0U>ITER 

0226 

OOCE 

0202 


AXI 

2' 


0227 

OOCF 

EA12 


STX 

TBP 

SAVE TBIP P0 INTER 

0228 

OODD 

C7FF 


LAM 

255 

DATA C01MT 

0229 

OODl 

9E28 


STA 

CNT 


0230 

00D2 

E61F 


LEK 

NI 

INPUT P0INTER 

0231 

00D3 

0202 

A1 

AXI 

2 

BUMP 00 UN TER 

0232 

00D4 

EES 1 

* 

STX 

IBP 

+ SAVS 

0233 

00D5 

B618 


LDA 

TETR 


0234 

00D6 

2188 


JAL 

T5 


0235 

00D7 

3400 


LDA 

©0 


023 6 

OODS 

E401 


LDX 

@1 


0237 

00D9 

1328 


LLX 

1 ' 


0238 

OODA 

1BS6 


LLR 

7 


0239 

00D3 

13A3 


LRX 

1 


0240 

OODO 

1356 


LLA 

7 


0241 

OODD 

10D6 


ARA 

7 


0242 

OODE 

F202 


JMP 

T6 


0243 

OODF 

B400 

75 

LDA ■ 

eo 


0244 

OOEO 

E401 


LDX 

61 


0245 

OOEl 

F9 00 

T6 

JST 

DPACC 


0246 

00E2 

0000 

TBP 

DATA 

0 


0247 

00E3 

0202 


AXI 

2 

BL^IP TH1P 

0248 

00E4 

EE02 


STX 

TBP 


0249 

00E5 

E69 2 


LDX 

IBP 


025 0 

00E6 

0110 


EAR 



0251 

00E7 

9000 


STA 

©0 


0252 

00E8 

9001 


STA 

©1 


0253 

00E9 

DE40 


IMS 

□vJT 

D0NE? 

0254 

OOEA 

F617 


JMP 

A1 

N0 ■ 

0255 

OOEB 

F72B 


RTW 

ACC 


0256 



* 




0257 



*VAIT. 

F0R ALIGNMENT 

PULSE 

0258 



* 




0259 

OOEO 

08 00 

TVKQ0N 

BIT 



0260 

OOED 

4006 


01 D 


DISABLE AUT0 

0261 

OOEE 

48 01 


SSN ■ 

01 

LIGHT 0FF? 

0262 

-00 EF 

F601 


JH-P - 

1 

N0 • - 

02.63 

OOFO 

OEOO 


SK-I 

* • 

YES BYTE 0N 

0264 

OOFl 

4901 


SEN 

: 01 

LIGHT 0N? 

0265 

00F2 

F60I 


JMP 

£- 1 

N0 

0266 

OOF3 

OFOO 


sm 

• • 

YES, BYTE 0FF 

0267 

00F4 

4004 


SEL 

: 04 

CLEAR FLAG 

0268 

00F5 

4005 


Cl E 


HvJABLE AUT0 

0269 

00F6 

F70A 


rtn 

TURN0N 
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0270 



* 




027 1 



*INPUT 

FR0M 

M0N SANTO 


0272 



CONVERT BCD T0 BINARY 

027 3 



* 




0274 

00F7 

05 00 

M0NS 

HNT 



■0275 

ODFB 

B633 


EDA 

TETR 


027 6 

00F9 

20SF 



T3 


0277 

OOFA 

C6FF 


LAP 

255 


0278 

00F3 

8E4B 


ADD 

DCT 


0279 

OOFC 

1050 


ALA 

1 


028 0 

OOFD 

E644 


LIK 

THEM 


025 1 

OOFE 

38 02 


JXN 

S+ 3 


025 2 

OOFF 

5A0S 


ADD 

mTR 


025 3 

0100 

F201 


JMP 

£+2 


028 4 

0101 

8A05 


ADD 

DPTR 


028 5 

0102 

0048 


TAX 



026 6 

0103 

B400 


LDA 

@0 


028 7 

0104 

E401 


LDlX 

@1 


0235 

0105 

F9 00 


JST 

DPFIX 


0289 

0106 

F70F 

- 

RTN 

M0NS 


029 0 

0107 

1002 

DPTR 

DATA 

: 1002 


029 1 

OIOS 

1202 

H-ITR 

DATA 

; 1202 


029 2 

0109 

49 C6 

T3 

SEN 

: C5 

FLAG 

029 3 

OlOA 

F601 


JMP 

S- 1 

N0 

0294 

0103 

5AC6 


I NX 

; C6 

YES^ INPUT T0 X 

029 5 

01 DC 

2843 


JX2 

£-3 

I GN0RE ZER0ES 

029 6 

OlOD 

C7 04 


LAM 

4" 

SET DIGIT C0UNT 

029 7 

OlOE 

9E65 


sta 

• CNT 


0293 

01 OF 

0110 


ZAR 


CLEAR TALLY 

0299 

0110 

9E63 


STA 

MSK2 


0300 

0111 

F206 


JMP 

M2 


0301 

0112 

1350 

Ml 

LLA 

1 

XIO 

0302 

0113 

9E66 


STA 

MSK2 


0303 

0114 

1351 


LLA 

2 


0304 

0115 

8E68 


ADD 

MSX2 


0305 

0116 

9E69 


sta 

MSX2 


0306 

0117 

0110 


ZAR 



0307 

01 18 

1303 

M2 

LLL 

4 

GET NEXT DIGIT 

0303 

0119 

8E6C 


ADD 

MSK2 

ADD T0 TALLY 

0309 

01 lA 

DE71 


IMS 

CNT 

LAST DIGIT? 

0310 

0113 

F609 


JMP 

M 1 

N0 

031 1 

01 1C 

0045 


TAX 



0312 

OllD 

0110 


ZAR 



0313 

01 IE 

F727 


RTI^ 

MOMS 

YES 

0314 




ST0RAGE 


0315 

01 IF 

SD3A 

MSG 

DATA 

: 5D3A 


0316 

0120 

C1C2 


TEXT 

•AB0RT LAST RUN? ’ 

• 

0121 

CFD2 






0122 

D4A0 






0123 

CCCl 






0124 

D3D4 






0125 

A0D2 






0126 

D5CE 






0127 

BFAD 







1^9 


0317 012S 8DGA ABT DATA :5D5A 

0318 0129 BlAO TEXT M RUN ABORTED,' 

0l2A D2D5 

012B CEAO 
0 12C C.IC2 
01 2D CFD2 
012E D41CS 
0 1 2F CAAC 

0319 0130 A0D2 RUNS TEXT ‘ RUNS KEPT* 

0131 D5CE 

0132 D3A0 

0133 CBC5 

0134 D0D4 

0320 0135 8IBA DATA :8D8A,0 

0136 0000 

0321 0137 8DSA MODE DATA :SDBA 

0322 0138 CDCF TEXT 'MONSANTO, TAPE OR IN VERSE? • 

0139 CED3 
013A CICE 
013B D4CF ' 

013C ACDA 
013D ClDO 
013E C5A0 
013F CFD2 

0140 A0C9 

0141 CED6 

0142 C5D2 

0143 D3C5 

0 1 44 3FA0 

0323 

0000 ERRORS 


END 
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0001 




MAI-I 

IKBj iefc bipt^xeq 

0002 




NAi'l 

0TT^ ERR» CHLF 

0003 




WAI-I 

0TL^ 0FPA 


0004 

0000 



REL 

0 


0005 






- 

0006 



¥ PANIC SWITCH 

ORIGINAL PAGE IS 

0007 



* 



OP POOR QUALITY 

0008 

0000 

FBOO 


JST 

* CMN D 


0009 

0001 


CMN D 

REF 



OOiO 







0011 



^{^KETBOARD INPUT 


OOlS 



* 

* 



0013 

0002 

08 00 

IKB 

ENT 



OOU 

0003 

4038 


SEL 

7j 0 

AUT0-ECH0 

0015 

0004 

4039 


SEL 

7^ 1 

KBD H0DE 

0016 

0005 

59 39 


RDA 

7 , 1 

READ 0W FLAG 

0017 

0006 

4030 


SEL 

7 , 4 

RESET 

0018 

0007 

F705 


RTN 

IKB 


0019 

0003 

03 00 

I ER 

H'JT 



0020 

0009 

F£07 


J ST 

IKB 


0021 

OOOA 

CODF 


CAI 

DF 

BACK ARR0U? 

0022 

0003 

FFOA 


JST 

*CMND 

YES 

0023 

0000 

003 A 


CAI 

;8A 

LINE FEED? 

0024 

OOOD 

FFOO 


JST 

* CMN D 

YES 

0025 

OOOE 

F706 


RTN 

lER 

NEl THERtRETURIvI 

0026 



* 


■ 

* 

0027 



^ PAPER 

TAPE 

; INPUT 


0028 



DSO 

= 0 

FOR TTY 


0029 



* 

1 

FOR HSR 


0030 







0031 

OOOF 

03 00 

BI FT 

H'JT 



0032 

0010 

53 01 

BIP2 

I SA 


READ SWITCHES 

0033 

0011 

13D0 


LRA 

1 

DSO UP FOR TTY 

0034 

0012 

220A 


J0S 

HSR 

D0WW F0R HSR 

0035 

0013 

49 3B 


SHJ 

7 , 3 

TTY busy? 

0036 

0014 

F604 


JMP 

BIP2 

YES 

0037 

0015 

403A 


SEL 

7, 2 

N0, STEP READER 

0038 

0016 

48 39 

UT 

SSN 

7^ 1 

FLAG? 

0039 

0017 

F203 


JMP 

IT' 

YES 

0040 

0018 

0150 


I AR 


NO^BUMP COUNT 

0041 

0019 

2149 


JAZ 

BIP2 

RESTART IF TIME UP 

0042 

OOlA 

F604 


JMP 

UT 

ELSE CHECK FLAG AGAIN 

0043 

OOlB 

58 38 

IT 

INA 

7 , 0 

INPUT FR0M TTY 

0044 

0010 

F70D 


RTN 

BIPT 

& RETURIvI 

0045 

OOID 

49 33 

HSR 

SH'J 

6 , 3 

HSR BUSY? 

0046 

OOIE 

F60E 


JMP 

BIP2 

yes 

0047 

OOlF 

4032 


SEL 

6; 2 

N0^ STEP READER 

0048 

• 

0020 

48 35 

I'JH 

SSN 

6^ 5 

FLAG? 

0049 

0021 

F203 


JMP 

IH 

YES 

0050 

0022 

0150 


lAR 


N0J BUMP C0UI^T 

0051 

0023 

2153 


JAZ 

.BIP2 

RESTART IF TIME UP 

00 52 

0024 

F604 


JMP 

WH 

ELSE CHECK FLAG AGAIN 

0053 

0025 

53 35 

IH 

INA 

6^ 5 

IN PUT' FROM HSR 
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0054 

0026 

F717 


RTN 

BIPT 

& RETURIvI 

0055 



* 




0056 



H^UAI t 

FOR EXECUTE SIGNAL 

0057 



* 




0053 

0027 

08 00 

XEQ 

PJT 



0059 

0028 

FE20- 


JST 

lER 

INPUT 

0060 

0029 

C03D 


CAI 

■ ;8D 

CARRIAGE RETURi^I? 

0061 

002A 

F703 


RTN 

XEQ 

YES, RETURN 

0062 

002B 

F603 


JMP 

X EQ+ 1 , 

N0, GET MORE 

0063 







0064 



Sic OUTPUT T0 

TTY 


0065 



SfC 




0066 

0020 

08 00 

0TT 

ENT 



0067 

002D 

4030 


SEL 

7, 4 

RESET INTERFACE 

0063 

002E 

6D3B 


URA 

7, 3 

WRITE ON NOT BUSY 

0069 

002F 

49 3B 


SEN 

7, 3 

DON E? 

007 0 

0030 

F601 


JMP 

$- 1 

N0 

0071 

0031 

F705 


RTN 

0TT 

YES 

007 2 



* 




007 3 



*O0MMAND ERROR EXIT 


0074 







0075 

0032 

03 00 

ERR 

ENT 



007 6 

0033 

C6DF 


LAP 

; DF 

PRINT ARROW 

0077 

0034 

FEOS 


o ST 

0TT 


0078 

0035 

FF34 


JST 

4=Ci’lND 

RESTART COMMAND 

0079 



* 




003 0 



=ic CARRI AGE- RETURN, L IN E 

FEED 

008 1 




■ 



003 2 

0036 

08 00 

CRLF 

ENT 



008 3 

0037 

C68D 


LAP 

: 8D 

CR 

003 4 

0038 

FEOC 


JST 

0TT 


008 5 

0039 

068 A 


LAP 

• 8 A 

LF 

008 6 

003A 

FEOE 


JST 

0TT 


008 7 

003B 

F705 


RTN 

CRLF 


0088 



* 




0089 



* OUTPUT TEXT FROM BUFFER 

009 0 



* 




009 1 

0030 

08 00 

0TE 

E'3T 



009 2 

003D 

8A0F 


ADD 

CAI 

MAKE COMPARE INSTRUCTION 

009 3 

003E 

9A06 


STA 

0T2 

a SAVE IT 

009 4 

003F 

9A09 


STA 

0T3 


009 5 

0040 

E704 


LEX 

*0TL 

GET TEXT POINTER 

009 6 

0041 

DE05 


IMS 

6TL 

SET RETURN ADDRESS 

009 7 

0042 

B400 

0T1 

LDA 

®0 

GET WORD 

0098 

0043 

1 1D7 


RRA 

8 

PRINT FIRST BYTE 

0099 

0044 

FE18 


JST 

0TT 


0100 

0045 

0000 

0T2 

CAI 

0 

last ONE 

0101 

0046 

F70A 


RTN 

OfL 

YES, return' 

0f02 

0047 

1157 


RLA 

8 

PRINT SECOND BYTE 

0103 

0048 

FEIO 


JST 

0TT 


0104 

0049 

0000 

0T3 

CAI 

0 

LAST ONE? 

0105 

004 A 

F70E 


RTN 

■0TL 

YES RETURN 

0106 

004B 

0128 


IX R 


BUMP POINTER 



0107 

0040 

F60A 


JMP 

0T1 

L00P 

0108 

004D 

GOOO 

CAI 

CAI 

0 


0109 



•T* 




0110 



•+ OUTPUT FLOATING POINT 

NUMBER 

01 1 1 







0112 

004E 

03 00 

0FPA 

ENT 



0113 

004F 

E701 


LD>1 

*0FPA 

GET P0 INTER 

01 14 

0050 

DE02 


IMS 

0FPA 

SET RETURJ ADDRESS 

01 15 

0051 

EAOl 


STX 

OPT 

SAVE POINTER 

01 16 

0052 

FBOE 


JST 

^i-.FAS 

CONVERT T0 ASCII 

01 17 

0053 

0000 

OPT 

DATA 

0 


01 18 

0054 

0059 


DATA 

BUF 


01 19 

0055 

0110 


Z AR 


SET E\ID FLAG 

0120 

0056 

FElA 


JST 

0TL 

PRINT NUMBER 

0121 

0057 

0059 


DATA 

BUF 


0122 

0058 

F70A 


RTIn] 

0FPA 


0123 

0059 

0000 

BUF 

RES 

8, 0 


0124 

0061 


FAS 

REF 



0125 

0000 

ERRORS 



n\ID 
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0001 





READ, CL EAR 

0002 





PUMCH, 

GRAPH 

0003 




SXTR 

IKB,0TT, BIPT 

OOOA 




DCTR 

0TL, CPEF, FPL0T, DFPA 

0005 




EXTR 

XA,XB,XC. 

, FAD 

0006 

0007 

0000 



REL 

0 


ORIGINAL PAGE Ib 

0008 



READ PAPER 

TAPE & 

ADD T0 BUFFER Qp pOQR QUALITY 

0009 



* 





0010 

0000 

08 00 

READ 

H\5T 




001 1 



vIWITIALIZE 

VARI ABL ES 


0012 

0001 

1 128 


■ RLX' 

1 


SET INDEX BIT 

0013 

0002 

1400 


S8V 




0014 

0003 

1 IAS 


RPX 

1 



0015 

0004 

EA33 


STX 

R5 


& SAVE POINTER 

0016 

0005 

EA34 


ST>C 

R5+2 



0017 

0006 

0528 


XRP 



RESET SUBSCRIPTS 

CO 18 

0007 

E9 00 
0000 


STX 

XA 



0019 

0008 

E9 00 


STX 

XB 





0000 






0020 

0009 

E9 00 
0000 


STX 

xc 



0021 

OOOA 

B236 


LDA 

TP 


SET INPUT BUFFER POINTER 

0022 

OOOB 

9A38 


STA 

MPT 



0023 

OOOC 

B235 


LDA 

CT 


SET INPUT C0UNT 

0024 

OOOD 

9A37 


STA 

CNT 



0025 



*SKIP 

L EADER 

& LABEL 



0026 

OOOE 

C63A 


LAP 

: 8A 


LINE FEED 

0027 

OCOF 

F9 00 
0000 


JST 

0TT 



0028 

0010 

F9 00 
0000 


J ST 

BIPT 



0029 

0011 

2141 


JA2 

S- 1 


SKIP LEADER 

0030 

0012 

F9 00 
0000 

R1 

JST 

0TT 


ECH0 LABEL 

0031 

0013 

F9 00 

« 

JST 

BIPT 


READ TAPE 



0000 






0032 

0014 

C09 2 


CAI 

: 92 


END 0F LABEL? 

0033 

0015 

F201 


JMP 

R2 


YES 

0034 

0016 

F604 


JMP 

R1 


N0 

0035 

0017 

F9 00 
0000 

R2 

JST 

BIPT 


READ TAPE 

0036 

0018 

GOFF 


CAI 

: FF 


FILE MARK? 

0037 

0019 

F201 


JMP 

R3 


YES' 

00 38 

OOIA 

F603 


JMP 

R2 


N0 

0039 



H«READ 

L00P 




0040 

00TB 

F9 00 

R3 

J'ST 

Bi'pr 


LOWER' BYTE 

• 


0000 






0041 

OOIC 

1BS7 


LLR 

8 


SAVE 

0042 

OOID 

F9 00 


JST 

BIPT 


UPPER BYTE 



0000 






0043 

OOIE 

1B07 


LLL 

8 


REST0RE 
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0044 

OOIF 

E224 


LDK 

MPT 

SAVE LOW BITS 

0045 

0020 

9001 


STA 

ei 


0046 

0021 

F9 00 


JST 

BIPT 

LOWER BYTE 



0000 





0047 

0022 

1B8 7 


LLR 

8' 


0048 

0023 

F9 00 


J ST 

BIPT 

UPPER BYTE 



0000 





0049 

0024 

1B07 


LLL 

8 


0050 



* CON VERT BASI C- F. P. T0 

CAI-F. P. 

0051 

0025 

2109 


JAZ 

r i ' ‘ 

0= 0= 0 ■ ‘ 

005S 

0026 

0048 


TA)C 


SAVE SI©I 

0053 

0027 

303 1 


JAP 

S+ 2 

ABS. VAL. 

0054 

0028 

0310 


WAR 



0055 

0029 

1^7 


LLR 

8 

REMOVE MSB 

0056 

002A 

1328 


LLX 

1 


0057 

002B 

SA17 


ADD 

D64 

FIX CHARACTERI STI C 

CD 58 

0020 

1B07 


LLL 

8 


0059 

002D 

1300 


LAO 


RECOVER SIGN 

0060 

002E 

1 IDO 


ERA 

1 


0061 

002F 

9BU 

Z 1 

STA 

j^MPT 

SAVE HI BITS 

0062 

0030 

DA13 


IMS 

MPT 

BWMP POINTER 

0063 

0031 

DAI 2 


IMS 

MPT 

TWICE! 

0064 

0032 

DAl2 


IMS 

CNT 

MORE? 

0065 

0033 

F618 


JMP 

R3 

YES 

0066 


' 

*ADD INPUT BUFFER TO BEAi-l ARRAY 

0067 



:^=PANI C 

SWITCH DISABLED 

FOR DURATION 

0065 

0034 

4006 


CID 


DISABLE AUTO 

0069 

0035 

B20C 


LDA 

CT 

RESET COUNTER 

0070 

0036 

9A0E 


STA 

CNT 


0071 

0037 

F9 00 

R4 

JST 

FAD 

F. P. ADD 

' 


0000 




• * 

0072 

0038 

0000 

R5 

DATA 

0^ ; 9400^ 0 


0039 

9400 






003A 

0000 





007 3 

003B 

D900 


IMS 

XA 

BUMP SUBSCRIPTS 



0000 


' 



0074 

0030 

D900 


IMS 

XB 




0000 





0075 

003D 

D9 00 


IMS 

xc 




0000 


■ 



0076 

003E 

DA 06 


IMS 

CNT 

MORE? 

0077 

003F 

F60S 


JMP 

R4 

YES 

0078 

0040 

F740 


RTW 

READ 


0079 


- 

STORAGE 



005 0 

0041 

1400 

TP 

DATA 

: 1400 


008 1 

0042 

ffoo 

CT 

DATA 

-256 


008 2 

0043 

0040 

D64 

DATA 

64 


008 3 

0044 

0000 

MPT 

DATA 

0 


003 4 

0045 

0000 

CNT 

DATA 

0 


008 5 



* 




008 6 



CL EAR 

BEAM 

ARRAY 


008 7 



* 


• 


0038 

0046 

08 00 

CLEAR 

EMT 
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0089 

0047 

B60S 


LDA 

CT 

SET C0UNT 

009 0 

0048 

9E03 


STA 

CNT 


009 i 

0049 

0110 


Z AR 


CLEAR 

009 2 

004 A 

9C00 

Cl 

STA 

GO 

HI BITS 

009 3 

004E 

9 C-0.1 


STA 

ei 

L0 BITS 

009 4 

004C 

C202 


A>CI 

2‘ 

bump P0 inter 

009 5 

004D 

DE08 


IMS 

CNT 

DON E? 

009 6 

004E 

F604 


JMP 

Cl 

N0 

009 7 

004F 

F709 


RTN 

CL EAR 


0098 



* 




0099 



* PUN CH 

BEAM 

ARRAY - 

BASIC F0PMAT 

0100 



* 


• 


0101 

0050 

08 00 

PUNCH 

HIT 



0102 

0051 

FA22 


JST 

LEAD 

L EADER 

0103 

0052 

B610 


LDA 

CT 

START C0U1>ITER 

0104 

0053 

9E0E 


STA 

CNT 


0105 

0054 

F9 00 

PI 

JST 

IK3 

ECHO LABEL 



0000 




- 

0106 

0055 

009 2 


CAI 

: 92 

CTRL /TAPE? 

0107 

0056 

F201 


JMP 

P2 

YES ' 

0108 

0057 

F603 


JilF 

PI 

N0 

0109 

0058 

C6FF 

P2 

LAP 

: FF 

PUI-JCH FILE MARK 

Olio 

0059 

F9 00 


JST 

0TT 

■ 



0000 





01 1 1 

005A 

B401 

P3 

LDA 

§1 

PUNCH L0 BITS 

01 12 

005B 

F9 00 


J ST 

0TT 




0000 





01 13 

005C 

13D7 


LRA 

8 


01 14 

005D 

F9 00 


JST 

0TT 




0000 





01 15 

005E 

B400 


LDA 

©0 

GET HI BITS 

01 16 

005F 

C202 


AXI 

2 

BUMP pointer 

0117 

0060 

EEIC 


STK 

MPT 

AND SAVE 

0118 



* CONVERT CAl 

“F, P. T0 

BASIC- F. P. 

01 19 

0061 

210A 


JAS 

' Z2' 

0= 0= 0 ■ ■ 

0120 

0062 

0048 


TA/C 


SAVE SIGN 

0121 

0063 

1350 


LLA 

1 

CLEAR A15 

0122 

0064 

1B8 7 


LLR 

8 

SPLIT 

0123 

0065 

9 622 


SUB 

D64 

FIX CHARACTERISTIC 

01 24 

0066 

1400 


SDV 


INSERT MSB 

0125 

0067 

1150 


RLA 

1 


Ol26 

0068 

lB06 


LLL 

7 

RE-F0 HI AT 

0127 

0069 

1329 


LLX 

2 

RECOVER SIGN 

0l28 

006A 

3201 


J0R 

$+2 

C0RRECT F0R IT 

0129 

006B 

0310 


WAR 



0130 

006C 

F9 0Q 

Z2 

JST 

0TT 

PUI'ICH HI BITS 



0000 




“ 

0131 

006D 

1B8 7 


LLR 

8 


0152 

006E 

F9 00 


JST 

0TT 




0000 





0133 

006F 

E62B 


LCX 

MPT 

REC0VER POINTER 

0134 

0070 

DE2B 


IMS 

'CNT 

DON E? ■ 

0135 

0071 

F617 


JMP 

P3 

N0 
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0136 

0072 

FAOl 


J ST 

LEAD 

YES^ PUNCH LEADER 

0137 

0073 

F723 


RTN 

PUTnICH 


0133 



* PUNCH 

5" 0F 

L EADER 


0139 

0074 

08 00 

LEAD 

ENT 



0140 

0075 

07 32 


LAM 

50 


0141 

0076 

9E31 


STA 

CNT 


0143 

0077 

0110 


EAR 



0143 

0078 

F9 00 

L2 

JST 

0TT 




0000 





0144 

0079 

DE34 


IMS 

CNT 


0145 

007A 

F602 


JMP 

L2 


0146 

0147 

007B 

F707 

'i' 

RTN 

LEAD 


0143 



=}=PL0T 

DATA ARRAY 


0149 



* 




0150 

007C 

03 00 

GRAPH 

ENT 



0151 

007D 

EA06 


sue 

T 

SAVE C0UNT POINTER 

0152 

007E 

0202 


AY I 

2 


0153 

007F 

EA03 


srx 

PT 

SAVE DATA POINTER 

0154 

008 0 

0110 


EAR 


PRINT COUNT 

0155 

008 1 

F9 00 


JST 

0TL 




0000 





0156 

003 2 

008A 


DATA 

CTY 


0157 

008 3 

F9 00 
0000 


J ST 

0FPA 


0158 

008 4 

0000 

T 

DATA 

0 


0159 

003 5 

F9 00 
0000 


J ST 

CRLF 


0160 

003 6 

C7FF 


LAM 

255 

PLOT DATA 

0161 

003 7 

F9 00 


JST 

FPL0T 




0000 





0162 

0038 

0000 

PT 

DATA 

0 


0163 

008 9 

F70D 


RTN 

GRAPH 


0164 

003 A 

8i^A 

CTX 

DATA 

; 8 AS A 


0165 

008 B 

C3CF 


TE>IT 

’COUNT = 

t 


008 C 

U5CE 






003 D 

D4A0 






008 E 

BDAO 





0166 

008 F 

0000 


DATA 

0 


0167 

0000 

ERRORS 



Hv D 





0001 







0002 



* FH-0T - FL0ATIWG POINT 

H-0TTER 

0003 



* 




000-^ 



* 

CALLING SEQUENCE; 


0005 



5)- 



UDA MDIM 

0006 






JST *FPL0T 

0007 






DATA ARRAY 

0008 






C RETURN) 

0009 



* 




0010 




NAM 

FH.0T 


001 1 

0000 



REL 

0 


0012 

0000 

08 00 

FPL0T 

HIT 



0013 

0001 

9 A4C 


STA 

NR0U 


0014 

0002 

9A4C 


STA 

CNT 


0015 

0003 

B703 


LDA 

xFPLOT 


0016 

0004 

9A1F 


STA 

RP3 


0017 

0005 

9A24 


STA 

RP4 


0018 

0006 

A249 


I0R 

B15 


0019 

0007 

9A05 


STA 

RPl 


0020 

0008 

9 ACS 


STA 

RP2 


0021 

0009 

0350 


ARP 



0022 

OOOA 

9B4F 


STA 

*XA 


0023 

0003 

F204 


JMP 

MOVE 


0024 

OOOC 

FB4E 

LPl 

JST 

*FCP 


0025 

OOOD 

0000 

RPl 

DATA 

O..MA)C 



OOOE 

0051 





0026 

OOOF 

2183 


JAL 

INC 


0027 

^ 0010 

FB4B 

HOVE 

JST 

*FMV 


0023 

0011 

0000 

RP2 

DATA 

0>MAK 



0012 

0051 





0029^ 

--0013 

D546 

INC 

IMS 

siiXA 


0030 

0014 

DA3A 


I'M S 

CNT 


0031 

0015 

F609 


JMP 

LPl 


'00.32 

0016 

FB41 


JST 

*CRLF 


0033 

0017 

FB40 


JST 

^CRLF 


0034 

" 0018 

C6BD 


LAP 

«— » 


0035 

0019 

FB3F 


JST 

*0TL 


0036 

00 1 A 

0061 


DATA 

SF 


0037 

OOlB 

FB43 


JST 

>i:FDV 


0033 

OOlC 

0055 


DATA 

F50..MAX., SCALE 


OOID 

0051 






OOIE 

0053 





0039 

OOIF 

FB40 


JST 

*0FPA 


0040 

0020 

0053 


DATA 

SCALE 


0041 

0021 

FB36 


JST 

*CRLF 


0042 

0022 

FB35 


JST 

*CRLF 


0043 

0023 

FB3C 

LP2 

JST 

*0FPA 


0044 

0024 

0000' 

RP3 

DATA 

0 


0045 

0025 

DEOl 


IMS 

RP3 


0046 

0026 

DE02 


IMS 

RP3 


0047 

0027 

C6A0 


LAP 

t f 


0048 

0028 

FB2E 


JST 

^0TT 


0049 

0029 

FB34 


JST 

*FMP 
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0050 

002A 

0000 

RP4 

DATA 

Oj SCALED max 


002B 

0053 





002C 

0051 




0051 

002D 

DE03 


IMS 

RP4 

0052 

002E 

DE04 


IMS 

RP4 

0053 

002F 

FB2D 


JST 

*FIX 

0054 

0030 

0051 


DATA 

MAXj CNT 


0031 

004F 




0055 

0032 

32 1C 


LDA 

CNT 

0056 

0033 

0048 


TAX 


0057 

0034 

303 5 


JAP 

P0S 

0058 

0035 

C6AD 


LAP 

» — » 

0059 

0036 

F320 


JST 

^^0TT 

0060 

0037 

C6B0 


LAP 

•O' 

0061 

00 33 

FBIE 


J ST 

^0TT 

0062 

0039 

F20E 


JMP 

OUT 

0063 

003A 

C6A0 

P0S 

LAP 

1 t 

0064 

0033 

F31B 


JST 

^0TT 

0065 

003C 

23 09 


JXZ 

MAJIK 

0066 

003D 

C6B0 


LAP 

•0’ 

0067 

003E 

F318 


JST 

V0TT 

0063 

003F 

0503 


WXR 


0069 

0040 

EAOE 


STX 

C'JT 

007 0 

0041 

C6AA 


LAP 

» 

0071 

0042 

F201 


JMP 

£+2 

0072 

0043 

FBI 3 


JST 

^0TT 

007 3 

0044 

DAOA 


IMS 

CNT 

0074 

0045 

F602 


JMP 

S-2 

0075 

0046 

C6AA 

MARK 

LAP 


0076 

0047 

F30F 


JST 

»0TT 

0077 

0048 

FBOF 

0UT 

JST 

4: CRLF 

0078 

0049 

DAO 4 


IMS 

NR0¥ 

0079 

004A 

F627 


JMP 

LP2 

008 0 

004B 

FBOC 


JST 

*CRLF 

008 1 

004C 

DE4C 


IMS 

FPL0T 

003 2 

004D 

F74D 


RTM 

FPL0T 

008 3 

004E 

08 00 

MR9U 

KLT 


008 4 

004F 

03 00 

CNT 

HLT 


008 5 

0050 

8000 

B15 

DATA 

:8000 

008 6 

0051 

0000 

MAK 

RES 

2, 0 

008 7 

0053 

0000 

SCALE 

RES 

2, 0 

0038 

0055 

4 343 

F50 

DATA 

; 4348^ 0 


0056 

0000 




0089 

0057 


0TT 

REF 


009 0 

0058 


CRLF 

REF 


009 1 

0059 


0Tu 

REF 


009 2 

005A 


XA 

REF 


009 3 

0053 


FCP 

REF 


009 4 

005C 


RIV 

REF 


009 5 

005D 


FIX 

REF 


009 6 

005E 


Fi'IP 

REF 


009 7 

005F 


FDV 

REF 


0093 

0060 


0FPA 

REF 
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0099 

0061 

D3C3 SF 

TEXT 

•sea. FACTOR 


0062 

ClCC 




0063 

8506 




0064 

C1C3 




0065 

D4CF 




0066 

D2A0 




0067 

BDAO 



0100 



Hsj D 




FL0ATIWG POINT PLOTTER 


0001 

0002 

0003 

0004 


:t:CRT - 

CALLING SEQUEMC 


0005 






0006 






0007 






0008 




NAM 

CRT 

0009 

0000 



REL 

0 

0010 

0000 

08 00 

CRT 

H'lT 


00 1 1 

0001 

C7FF 


LAM 

255 

0012 

0002 

9A5F 


STA 

NR0U 

0013 

0003 

9 ASF 


STA 

CNT 

0014 

0004 

49 C4 


SEN 

24:» 4 

0015 

0005 

F601 


JMP 

S” 1 

0016 

0006 

40C7 


SEL 

24; 7 

0017 

0007 

C202 


AXI 

2 • 

0018 

0008 

0030 


TXA 


0019 

0009 

A261 


I0R 

B15 

0020 

OOOA 

9A0B 


STA 

RPl 

0021 

OOOB 

9A0E 


STA 

RP2 

0022 

OOOC 

9A10 


STA 

RP3 

0023 

OOOD 

9A13 


STA 

RP4 

0024 

OOOE 

9A36 


STA 

RP5 

0025 

OOOF 

0110 


2 AR 


0026 

0010 

9A54 


STA 

MIN 

0027 

001 1 

9A54 


STA 

MIN+ 1 

0023 

0012 

0350 


ARP 


0029 

0013 

9B5B 


STA 


0030 

0014 

F204 


JMP 

MOVE 

0031 

0015 

FB5A 

LPl 

J ST 

*FCP 

0032 

0016 

0000 

RPl 

DATA 

O^MAX 


0017 

0067 




0033 

0018 

218 3 


J AL 

IN Cl 

0034 

0019 

FB57 

MOVE 

JST 

=iJFMV 

0035 

OOlA 

0000 

RP2 

DATA 

O^MAX 


OOlB 

0067 




0036 

OOlC 

FB53 

INCl 

JST 

*FCP 

0037 

OOID 

0000 

RP3 

DATA 

0^ MIN 


OOIE 

0065 




0038 

OOIF 

308 3 


JAP 

INC2 

0039 

0020 

FB50 


JST 

*FI'IV 

0040 

0021 

0000 

RP4 

DATA 

O^MIN 


0022 

0065 




0041 

0023 

DB43 

IWC2 

IMS 

*XA 

0042 

0024 

DA3E 

■ 

IMS 

CNT 

0043 

0025 

F610 


JMP 

LPl 

0044 



* 



00*45 

0026 

FB4F 


JST 

*CRLF 

0046 

0027 

FB4E 


JST 

H^CRLF 

0047 

0028 

C6BD 


LAP 

t ^ t 

0048 

0029 

FB4D 


JST 

■jROTL 

0049 

002A 

0079 


DATA 

SF 


d ST CRT 
c return}) 



0050 



* 



0051 

002B 

FB46 


JST 

^FSB 

0052 

002C 

0067 


DATA 

HAX,MIN^MAX 


002D 

0065 





002E 

0067 




0053 

002F 

FB43 


JST 

=!cFDV 

0054 

0030 

0060 


DATA 

F250jMAX^ scal: 


0031 

0067 





0032 

0069 




0055 

0033 

FB44 


JST 

=5=0 FPA 

0056 

0034 

0069 


DATA 

SCALE 

0057 

0035 

FB40 


J ST 

*CRLF 

0058 



❖ 



0059 

0036 

FB3D 


JST 

*H-1P 

0060 

0037 

0069 


DATA 

SCALE>MIN^MAX 


0038 

0065 





0039 

0067 




0061 

003A 

F33A 


JST 

* FIX 

0062 

003B 

0067 


DATA 

MA}<^MIN 


003C 

0065 




0063 

003D 

B227 


LDA 

MIM 

0064 

003E 

0310 


WAR 


0065 

003F 

308 1 


JAP 

S+2 

0066 

0040 

0110 


EAR 


0067 

0041 

9A22 


STA 

X 

0068 

0042 

0350 


ARP 


0069 

0043 

9B2B 


STA 

iHXA 

0070 






0071 

0044 

FB2F 

LP2 

J ST 

* FMP 

0072 

0045 

0000 

RP5 

DATA 

O^ SCALELrMAX 


0046 

0069 





0047 

0067 




0073 

0048 

DB26 


IMS 

*XA 

0074 

0049 

FB2B 


J ST 

*FIX 

0075 

004A 

0067 


DATA 

MA3C, CNT 


004B 

0063 




007 6 

004C 

E217 


LDK 

X 

0077 

004D 

B215 


LDA 

CNT 

0078 

004E 

49 04 

DOT 

SH'3 

24^ 4 

0079 

004F 

F601 


JMP 

$- 1 

008 0 

0050 

49 02 


SH\I 

24» 2 

008 1 

0051 

F601 


JMP 

S- 1 

008 2 

0052 

6EC2 


OTJC 

24, 2 

008 3 

0053 

2107 


JAZ 

DONE 

003 4 

0054 

318 3 


JAG 

P0S 

008 5 

0055 

00 AS 


EKR 


008 6 

0056 

. 0150 


I AH 


0087 

0057 

F609 


JMP 

DOT 

0038 

0058 

0128 

P0S 

IXR 


0089 

0059 

OODO 


DAR 


009 0 

005A 

F60C 


JMP 

DOT 

009 1 

005B 

3208 

DOME 

LDA 


009 2 

0050 

8A1 r 


ADD 

H256 



l62 


009 3 

005D 

9A06 


STA 



009 4 

005E 

DAO 3 


IMS 

NR0W 


009 5 

005F 

F6 IB 


JMP 

LP2 


009 6 

0060 

FB15 


J ST 

>i:CRLF 


009 7 
0098 

0061 

F761 


RTW 

CRT 


0099 

0062 

0000 

I'lROU 

DAiTA 

0 


0100 

0063 

0000 

C\!T 

DATA 

0 


0101 

0064 

0000 

X 

DATA 

0 


0102 

0065 

0000 

MIN 

RES 

2^ 0 


0103 

0067 

0000 

max 

RES 

2^ 0 


0104 

0069 

0000 

so ALE 

RES 

2^ 0 


0105 

006B 

8000 

B15 

DATA 

;8 000 

of JS 

^ QUALHY 

0106 

0060 

006D 

447A 

0000 

F250 

DATA 

: 44 7A, 0 

0107 

006E 

0100 

K256 

DATA 

256 


0103 







0109 

006F 


XA 

REF 



01 iO 

0070 


FOP 

REF 



01 1 1 

007 1 


R-IV 

REF 



0112 

0072 


FSB 

REF 



0113 

0073 


FDV 

REF 



01 U 

0074 


RIP 

REF 



0115 

0075 


FIX 

REF 



0116 

0076 


CRLF 

REF 



0117 

0077 


ETL 

REF 



01 13 

0073 


0FPA 

REF 



01 19 

0120 

0079 

D3C3 

SF 

TEXT 

• SCALE 

FACTOR = ' 


007A 

ClCC 






0073 

05 AO 






0070 

0601 






007D 

03D4 






007E 

CFD2 






007F 

AOBD 





0121 

0122 

0000 

ERRORS 


* 

END 
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0001 
0002 
00 0 3 

0004 

0005 

0000 

03 0 0 

VI Ei-J 

NA'I 

DCTR 

NAl 

DCTR 

D'JT 

0006 

0001 

9A7 1 

DI SP 

STA 

0007 

0002 

G202 


AXI 

0008 

0003 

EA66 


STM 

0009 

0004 

EA66 


SIX 

0010 

0005 

C7FF 


LAT-l 

0011 

0006 

9A6A 


STA 

0012 

0007 

9A64 


STA 

0013 

0008 

0110 


ZAR 

0014 

Uuu9 

9A63 


^‘lA 

0015 

OOOA 

9A63 


STA 

0016 

0003 

9A63 


STA 

0017 

OOOC 

9A63 


STA 

0018 

GOOD 

B400 

ST CM 

LDA 

0019 

OOOS 

D260 


CMS 

0020 

OOOF 

F20A 


JMP 

0021 

0010 

F205 


JMP 

0022 

0011 

B401 


LDA 

0023 

0012 

D25D 


CMS 

0024 

0013 

F206 


JMP 

00 2 5 

0014 

9A53 


STA 

0026 

0015 

F21 i 


JMP 

0027 

0016 

9A58 

MAC 1 

STA 

0023 

0017 

3401 


LDA 

0029 

0018 

9A57 


STA' 

0030 

0019 

F20D 


JMP 

0031 

OOIA 

3400 

MIN 1 

LDA 

0032 

OOIB 

D251 


C'IS 

0033 

OOlC 

F206 


JMP 

0034 

OOID 

F209 


JMP 

0035 

OOIE 

3401 


LDA 

0036 

OOIF 

D24E 


CMS 

0037 

0020 

F202 


JMP 

0038 

0021 

F205 


JMP 

0039 

0022 

F204 


JMP 

0040 

0023 

3400 

MI M2 

LDA 

0041 

0024 . 

9 A48 


STA 

0042 

0025 

3401 


LDA 

0043 

0026 

9A47 


STA 

0044 

0027 

E242 

FI CM 

LDX 

0045 

0028 

C202 


AXI 

0046 

0029 

EA40 


STK 

0047 

002A 

DA46 


IMS 

0048 

0023 

F61E 


JMP 

0049 

0050 

002C 

3242 

* 

LDA 

0051 

002D 

E242 


LDC 

0052 

002E 

F9 00 


JST 



0000 




VI El-' 

DPi-J:, DPSUE:> DPDI V: ^ DPSI: 

DI SP 

RTTP^ RT0P 

VEVR 
2 ■ 

RPID P0INTER 

RP2D 

255 

CNT 

NR0W 

XivITIPLIZE MISC i5-ID MIW 

MIN 
MIN+ 1 
MiK 
Mi^X+ 1 

eo 

M/SC 

MINI AKMP?: 

Mi«i Al>M/X 

Q1 Al=M^5>: 

MAX+ 1 

MINI" A1=MAC A2<MiX+l 

MAK+1 Ai=MAC A2>MPX+1 

FI CM ■ 

MA>I 

§1 

M^+ 1 
FI 0^1 ■ 

go 

MIN 

MIN2 A1<MIN 

FI CM A1>MIN 

§1 Al=MIN 

MIN+1 

MIN2 Ai=MIN A2<MIN+1 

FlCr^ Al=MIN A2>MIN+1 

FI CM A1=M1N A2=MIN+1 

eo 

MIN 

@1 

MlN+ 1 
RPID 
2 ■ 

RPID 
CNT 
ST CM 

M/X 
MAX+ 1 
DPSUB; 
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0053 

002F 

006D 


DATA 

MIN 

0054 

0030 

1050 


ALA 

1 

0055 

0031 

1329 


LLX 

2 

0056 

0032 

3201 


J0R 

S+2 

0057 

0033 

0150 


I AR 


005S 

0034 

13A3 


LRX 

1 

0059 

0035 

9A39 


STA 

MAX 

0060 

0036 

EA39 


STX 

MAX+ 1 

0061 






0062 

0037 

3235 


LDA 

MIN 

0063 

0038 

E235 


LDK 

MIN+ 1 

0064 

0039 

F9 00 


JST 

DPDI V; 



0000 




0065 

003A 

006F 


DATA 

MAX 

0066 

0033 

9A2D 


STA 

TEM 

0067 

0030 

lODl 


ARA 

2 

0068 

003D 

0310 


WAR 


0069 

003E 

8 A2A 


ADD 

TE'I 

0070 

003F 

10D5 


ARA 

6 

0071 



5k 



007 2 

0040 

0310 


WAR 


007 3 

0041 

303 1 


JAP 

S+2 

0074 

0042 

0110 


EAR 


0075 

0043 

9A2E 


STA 

X 

007 6 






0077 

0044 

E226 


L3X 

RP2D 

0078 

0045 

B400 

LPD 

L DA 

eo 

0079 

0046 

E401 


LDK 

©1 

003 0 

0047 

F9 00 


J ST 

DPDI V; 



0000 




003 1 

0048 

006F 


DATA 

MAX 

003 2 

0049 

9A1F 


STA 

TEM 

003 3 

004A 

lODl 


ARA 

2" 

003 4 

0043 

0310 


WAR 


003 5 

0040 

SAIC 


ADD 

TH'l 

003 6 

004D 

10D5 


ARA 

6 

003 7 



ik 



0038 

004E 

E223 


LEK 

X 

0039 

004F 

49 04 

DOT 

SH'] 

24j 4 

009 0 

0050 

F601 


JMP 

S- 1 

009 1 

0051 

49 02 


SHvJ 

241 2 

009 2 

0052 

F601 


JMP 

$-1 

009 3 

0053 

6EC2 


0TX 

241 2 

009 4 

0054 

2107 


JAZ 

DONE 

009 5 

0055 

318 3 


JAG 

P0S 

009 6 

0056 

00 AS 


DXR 


009 7 

0057 

0150 


lAR 


0098 

0053 

F609 


JMP 

DOT 

009 9 

0059 

0123 

P0S 

IXR 


0100 

005A 

OODO 


DAR 


0101 

0053 

F600 


JMP 

DOT 

0102 





• 

0103 

0050 

3215 

DONE 

LDA 

X 
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0104 

005D 

8A16 

0105 

005E 

9A13 

0106 

005F 

E20B 

0107 

0060 

C202 

0103 

0061 

EA09 

0109 

00‘62 

DA09 

0110 

0063 

F61E 

01 1 1 

0064 

B20E 

0112 

0065 

2101 

0113 

0066 

FlOO 

0000 

0114 

0067 

FlOO 

■* ' 


0000 

0115 

0068 

F768 

0116 

0069 

0000 

0117 

006A 

0000 

01 18 

006B 

0000 

0119 

006C 

0000 

0120 

006D 

0000 

0121 

006F 

0000 

0122 

0071 

0000 

0123 

0072 

0000 

0124 

0073 

0000 

0125 

0074 

0100 

0126 



0000 

ERRORS 




ADD 

H256 


STA 



LDX 

RP2D 


AXI 

2 


SIX 

RP2D 


TM S 

N-R0U 


JMP 

LPD 


LDA 

VEVR 


J AE 

S+2 


JMP 

RT0P 


JMP 

RTTR 


RTN 

VIEW 

TEM 

DATA 

0 

RPID 

DATA 

0 

RP2D 

DATA 

0 

NR0¥ 

DATA 

0 

MIN 

RES 

2/ 0 

MAX 

RES 

2 , 0 

CNT 

DATA 

0 

X 

DATA 

0 

VEVR 

DATA 

0 

H256 

DATA 

EJD 

256 


ORIGINAi; PAGE IS 
OP POOH QUALITY 



0001 




WAI4 

DPAOO 

0002 




NAii 

DPFLT 

0003 




NAl'l 

DPFIX 

0004 




EKTR 

DPM: 

0005 




EKTR 

DPM 

0006 

0000 



REL 

0 

0007 






0003 



'^=D0UELE PRECISI0W ACCUi'lULATE 

0009 






0010 

0000 

03 00 

DPACC 

EWT 


0011 

0001 

9a10 


STA 

TEMP 

0012 

0002 

0030 


rxA 


0013 

0003 

E703 


LDX 

* DPAOO 

00 1 4 

0004 

1200 


R0V 


0015 

0005 

SCOl 


ADD 

@I 

0016 

0006 

8200 


AUD 

MASK 

0017 

0007 

9 COl 


STA 

©1 

0018 

0008 

B209. 


LDA 

TE-IP 

0019 

0009 

3204 


J0R 

DCl 

0020 

OOOA 

1200 


R0V 


0021 

OOOB 

0150 


I AR 


0022 

0000 

3201 


J0R 

DCl 

0023 

OOOD 

0110 


Z AR 


0024 

OOOh 

8 000 

DCl 

ADD 

©0 

0025 

OOOF 

9 000 


STA 

90 

0026 

0010 

DEIO 


IMS 

DPACC 

0027 

0011 

F71 1 


RTw 

DPACC 

0028 

0012 

0000 

TRIP 

DATA 

0 

0029 
% • 

0013 

7FFF 

nASK 

DATA 

; 7FFp 

sj ^ 

0031 



DOUBLE PREOISI0N T0 OAI-F. 

0032 



* 


•• • • 

0033 

0014 

08 00 

DPFLT 

EMT 


0034 

0015 

33 01 


JXN 

DPFl 

0035 

0016 

2113 


JAZ 

DPF4 

0036 

0017 

BA13 

DPFl 

EI'IA 

BITS 

0037 

0018 

0110 


ZAR 


0038 

0019 

BAl 1 


E'lA 

BITS 

0039 

OOIA 

9Al 1 


STA 

SI GN 

0040 

OOlB 

308 1 


JAP 

DPF2 

0041 

OOlC 

F9 00 


JST 

DPN: 

■ 

■ 

0000 




0042 

OOID 

1328 

DPF2 

LLX 

1 

0043 

OOIE 

IBOO 

DPF3 

LLL 

1 

0044 

OOIF 

DAOB 


IMS 

BITS 

0045 

0020 

3242 


J0R 

DPF3 

0046 

0021 

1BS8 


LLR 

9 

0047 

0022 

BA 08 


E-IA 

BITS 

00-48 

0023 

0310 


MAR 


0049 

0024 

8 AOS 


ADD 

D160 

0050 

0025 

1356 


LLA 

7' 

0051 

0026 

A204 


10R 

BITS 

0052 

0027 

BA04 


EMA 

SIGN 



005 3 

0028 

0054 

0029 

0055 

002A 

0056 

002B 

0057 

00.2C 

0053 

002D 

0059 

D02E 

0060 

002F 

0061 

0030 

0062 

0031 

0063 

0032 

0064 

0033 

0065 

0034 

0066 

0035 

0067 

0036 

0068 

0037 

0069 

00 38 

007 0 

0039 

0071 

00 3A 

0072 

003B 

007 3 

0030 

0074 

003D 

0075 

003E 

0076 

003F 

0077 

0040 

0073 

0041 

0079 

0042 

008 0 

0043 

003 1 

0044 

003 2 

0045 

008 3 

0046 

003 4 

0047 

003 5 

0048 

008 6 

0 049 

008 7 

004A 

0038 

004B 

0039 

0040 

009 0 

004D 

009 1 

004E 

009 2 

004F 

009 3 

0050 

009 4 
0000 

ERRORS 


8205 


A202 


F716 

DPF4 

0000 

BITS 

O.O.OJ) 

SI GN 

00 AO 

D160 

8 000 

DPM 

08 00 

DPFIX 

9E04 


8210 


13D6 


9606 


9E09 


B609 


3 213 


A218 


BEOD 


2109 


319 0 


EEIO 


1B07 


1330 

DPFX 1 

DE13 


F602 


IBOO 


13i^ 


F201 


3E18 

DPFX2 

BE 13 


303 3 


BEIA 


F9 00 


0000 


F201 


3E1D 

DPFX3 

F71B 


01 IS 

DPFX4 

1400 


F71E 


7FFF 

DPP 

007F 

MATIT 

008 0 

SBIT 


AND 

DPM 

I0R 

SIGI^ 

RTW 

DPFLT 

DATA 

0 

DATA 

0 

DATA 

160 

DATA 

;8000 

EJT 

STA 

SIGN 

AND 

DPP 

LRA 

7 

SUB 

D160 

STA 

BITS 

LDA 

SIGN 

AND 

MAxiT 

I0R 

SBIT 

EM A 

BITS 

JAZ 

DPFX2 

JAG 

DPFX4 

H'lA 

BITS 

LLL 

5 

LLR 

1 

IMS 

BITS 

JMP 

DPFX 1 

LLL 

1 

LPJC 

1 

JMP 

S+2 

E/IA 

BITS 

E-IA 

SIGN 

JAP 

DPFX3 

H'l A 

SI GM 

JST 

DR-J: 

JMP 

S+2 

HI A 

SIGN 

RTI'J 

DPFIX 

ZAX 


S0V 

RTW 

DPFIX 

DATA 

: 7FFF 

DATA 

: 7F 

DATA 

:80 

H\JD 
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APPENDIX C 
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0001 




NAM 

CMN D 


0002 




EXTR 

I EFL 0TT^ 0TL^ CRL F 

0003 




E<TR 

ERR,XEQ 


0004 

0100 



ABS 

: 100 


0005 



H.-PRI-N T 

N A>1 £ 



0006 

0100 

0110 


EAR 



0007 

0101 

F9 00 


JST 

0TL 




0000 





0008 

0102 

012B 


DATA 

NA4E 


0009 

0103 

F201 


JMP 

S+2 


0010 



* 




0011 



* COMMAND INTERPRETER 


0012 



* 


• 


0013 

0104 

08 00 

a-IN D 

E\3T 



001 4 

0105 

OAOO 


EIN 



0015 

0106 

4005 


Cl E 


EvIABLE PAnJIC BUTTON 

C3016 

0107 

F9 00 


JST 

CRLF 

PRINT PROMPT CriARACTER 



0000 

• 




0017 

0108 

C687 


LAP 

;87 


0018 

0109 

F9 00 


J ST 

0TT 




0000 





0019 

OlOA 

C6AA 


LAP 

' 


0020 

0103 

F9 00 


JST 

0TT 




0000 





0021 

OlOC 

0108 


2XR 



0022 

OlOD 

F9 00 


J ST 

I ER 

INPUT COMMAND 



0000 





0023 

OlOE 

C0D4 


CAI 

• > 


0024 

01 OF 

E213 


LD>C 

TRANS 


0025 

0110 

C0D2 


CAI 

'R' 


0026 

Oil 1 

E212 


LDX 

READ 


0027 

0112 

C0C3 


CAI 

• C' 


0028 

0113 

E21 1 


LD)C 

clear 


0029 

0114 

38 OA 


JXN 

0KF 


0030 

0115 

CODO 


CAI 

.pi 


0031 

0 116 

E20F 


L ESC 

Pm CH 


0032 

0117 

C0C7 


CAI 

’ G’ 


0033 

01 13 

E20F 


LDX 

FPC0T 


0034 

0119 

C0D6 


CAI • 

•V 


0035 

01 lA 

E20E 


LD>C 

CRT 


0035 

01 IB 

C0C4 


CAI 

"D’ 


0037 

01 1C 

E20A 


LEK 

DATUM 


0033 

01 ID 

38 01 


JXN 

0KF 


0039 

01 IE 

F9 00 

BAD 

JST 

ERR 




0000 



. 


0040 

01 IF 

EAOA 

0KF 

STX 

FUNC 

SAVE CALL ADDRESS 

0041 

0120 . 

F9 0.0, 


JST . 

XEQ 

WAIT F0.R. G0 . 

• 


0000 





0042 

0121 

FB08 


JST 

* FUN C 

CALL FUJJCTI0N 

004 3 

0122 

F61D 


JMP 

CMN Df 1 


0044 

0123 


TRANS 

REF 



0045 

0124 


READ 

REF 



0046 

0125 


CLEAR 

REF 
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0047 

0126 


PUI'J Crf 

REF 


0043 

0127 


DATUl-I 

REF 


0049 

0123 


FPL0T 

REF 


0050 

0129 


CRT 

REF 


0051 

012A 

0000 

FUI'JC 

DATA 

0 

0052 

0 12B 

B1B5 

NAI'lE 

TEXT 

’ 15>k255 foul HTS’ 


012C 

AAB2 





012D 

B5B5 





012E 

A0C6 





012F 

D5CC 





0130 

CCAO 





0131 

C8D4 





0132 

D3A0 




0053 

0133 

8D3A 


DATA 

;8D3Aj 0 


0134 

0000 




0054 




EMD 



0000 ERRORS 
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0001 




MPI-l 

TRACIS 

0002 




NAtl 

RTTR 

0003 




EXTR 

DI SP 

0004 

0000 



REL 

0 

0005 

0000 

08 00 

TRAN'S 

EI'JT 


0006 

0001 

C6BF 


LAP 

1 2 1 

0007 

0002 

F9 00 


JST 

*0TL 



8 127 




0008 

0003 

013E 


DATA 

M0TA 

0009 

0004 

F9 00 


JST 

*IKB 



8 128 




0010 

0005 

1200 


R0V 


001 1 

0006 

GOOD 


GAI 

*M ’ 

0012 

0007 

F203 


JNP 

M 0N 

0013 

0008 

C0D4 


CAI 

T t 

0014 

0009 

F204 


JMP 

TAP 

0015 

OOOA 

F9 00 


JST 

*ERR 

* 


8 12A 




0016 

OOOB 

F9 00 

M0N 

JST 

*CRLF 

■* 


8 129 




0017 

OOOC 

0108 


ZXR 


0013 

OOOD 

F20F 


JMP 

BEG 

0019 

OOOE 

F9 00 

TAP 

JST 

^ CRL F 



8 129 




0020 

OOOF 

C63A 


LAP 

: 8A 

0021 

0010 

F9 00 


JST 

5t:0TT 



8 126 




0022 

0011 

F9 00 


JST 

:t=BIPT 


” " 

8 125 




0023 

0012 

2l4l 


JAZ 

S- 1 

0024 

0013 

F9 00 

T1 

JST 

'^077 



8 126 




0025 

0014 

F9 00 


JST 

*BI PT 



8 125 




0026 

0015 

C09 2 


CAI 

:92 

0027 

0016 

F201 


JMP 

T2 
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8A27 


ADD 

MSK2 

009 5 
009 6 

0050 

005D 

9A26 

0110 


STA 

ZAR 

MSK2 

009 7 

005E 

1B03 

M2 

LLL 

4 

0098 

005F 

8A23 


ADD 

MSK2 

0099 

0060 

DA2-1 


IM-S 

OWT 

0100 

0061 

F609 


JMP 

M 1 

01*04 

0062 

F714 


RTN 

DOWS 

0102 

0063 

DAIA 

DLP2 

IMS 

DDT 

0103 

0064 

F63A 


JMP 

DLP 

0104 

0105 

0065 

0066 

F629 

0TL 

JMP 

REF 

’CLPl 



j-/y 


0106 0067 8C8A DIRE DATA :8D3A 

0107 0068 C6CF TEXT ’ FORWARD^ BACKWARD 0R BOTH?' 

0069 D2D7 
006A C1D2 
006B C4AC 
006C C2C1 
006D C3CB 
006E D7CI 
006F D2G4 

0070 AOCF 

0071 D2A0 

0072 C2CF 
007 3 D4C8 
0074 BFAO 


0108 

0075 


CRLF 

REF 


0109 

0076 


IKB 

REF 


0110 

0077 


IBC 

REF 


01 1 1 

0078 


ERR 

REF 


0112 

0079 

0000 

DIRl 

DATA 

0 

01 13 

007A 

FOOD 

Z CTl 

DATA 

“408 5 

01 14 

007B 

0000 

CTl 

DATA 

0 

0115 

007C 

0000 

IB 

DATA 

0 

0116 

007D 

0000 

DIRO 

DATA 

0 

01 17 

007E 

0000 

DDT 

DATA 

0 

oi is 

007F 

0000 

DDT2 

DATA 

0 

01 19 

003 0 

0000 

DAT 

DATA 

0 

0120 

008 1 

0000 

nJT2 

DATA 

0 

0121 

003 2 

0000 

CHT 

DATA 

0 

0122 

003 3 

0000 

MSK2 

DATA 

0 

0123 

008 4 

0000 

R0W 

DATA 

0 

0124 




EID 


0000 

ERRORS 









lOU 


QOOl 




NAl-1 

IKB^ lER^ BIPT,XEQ 


0002 




WAl^ 

0TT^ ERR, CRLF 


0003 




N AT'l 

0TLJ0FPA. 0DEC 


0004 

0000 



REL 

0 



0005 








0006 



* PA^JIC SWITCH 



0007 








00 OS 

0000 

FBOO 


J ST 

:t:CMND 



0009 

0001 


CMND 

REF 




0010 



* 





0011 



5SKE7B0ARD INPUT 



0012 



* 





0013 

0002 

03 00 

IKB 

ENT 




0014 

0003 

403 


SEL 

7 ^ 0 

AUT0-ECH0 


0015 

0004 

4039 


SEL 

7 , 1 

KBD M0DE 


0016 

0005 

59 39 


RDA 

7 , 1 

READ 0N FLAG 


0017 

0006 

4030 


SEL 

7 , 4 

RESET 


0018 

0007 

F705 


RTN 

1KB 



0019 

0003 

08 00 

lER 

ENT 




0020 

0009 

FE07 


JST 

IKB 



0021 

OOOA 

CODF 


CAI 

: DF 

BACK ARROW? 


0022 

OOOB 

FFOA 


JST 

^ CMN D 

YES 


0023 

0000 

COS A 


CAI 

; 8 A 

LINE FEED? 


0024 

GOOD 

FFOC 


J ST 

*CMND 

YES 


0025 

OOOE 

F706 


RTM 

lER 

NEITHER, RETURN 


0026 



* 





0027 



* PAPER 

TAPE 

INPUT 



0028 



^ DSO = 

= 0 F0R TTY 



0029 



« 1 FOR HSR 




0030 



* 





0031 

OOOF 

03 00 

DIPT 

EfJT 




0032 

0010 

53 01 

BI P2 

ISA 


READ SWITCHES 


0033 

0011 

13D0 


LRA 

1 

DSO UP F0R TTY 


0034 

0012 

220A 


J0S 

HSR 

D3UN F0R HSR 


0035 

0013 

49 3B 


SEN 

7 j 3 

TTY BUSY? 


0036 

0014 

F604 


JNP 

BIP2 

YES 


0037 

0015 

403A 


SEL 

7^ 2 

N0^ STEP READER 


0038 

0016 

48 39 

UT 

SSN 

7 , 1 

FLAG? 


0039 

0017 

F203 


JMP 

IT 

YES 


0040 

0018 

0150 


lAR 


N0^ BU^IP C0U\IT 


0041 

0019 

2149 


JAZ 

BIP2 

RESTART IF TIME 

UP 

0042 

OOlA 

F604 


JMP 

UT 

ELSE CHECK FLAG 

AGAIN 

0043 

OOIB 

58 38 

IT 

INA 

7 , 0 

INPUT FR0M TTY 


0044 

OOlC 

F70D 


RTN 

BIPT 

& RETUH'J 


0045 

OOID 

49 33 

HSR 

SHxJ 

6; 3 

HSR BUSY? 


0046 

OOIE 

F60E 


JIIP 

BIP2 

YES 


0047 

OOIF 

^4032 


SEL 

6 , 2 

N0^ STEP READER 


0048 

0020 

48 35 

vrH 

SSN 

6^ 5 

FLAG? 


0049 

0021 

F203 


JMP 

IH 

YES 


0050 

0022 

0150 


I AR 


N0^ BU-IP COUNT 

4 

0051 

0023 

2153 


JAZ 

BIP2 

RESTART IF TIME 

UP 

0052 

0024 

F604 


JtiP 

■ l-TH 

ELSE CHECK FLAG 

AGAIN 

00 5 3 

0025 

53 35 

IH 

INA 

6 , 5 

INPUT FR0M HSR 




0054 

0055 

0056 

0057 
00 5B 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 

0067 

0068 
00 69 

007 0 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 
003 0 

008 1 
008 2 
008 3 
008 4 
008 5 
003 6 

008 7 
0038 
0089 

009 0 
009 1 
009 2 
G09 3 
009 4 
009 5 
009 6 
009 7 
0093 
009 9 
0100 
0101 

0102 

0103 

0104 

0105 

0106 


of & PAGi; IS 
QVAJjjjyi 


l8l 


0026 

F717 

'■a 

^ UAI T 
in 

XEQ 

RTM 

BIPT 

& RETURN 



F0R EXECUTE SIGNAL 

0027 

08 00 

ENT 



0028 

FE20 


JST 

lER 

INPUT 

0029 

003 D 


CAI 

:8D 

OARRIAGE RETURN? 

002A 

F703 


RTW 

XEQ 

YES> RETURN 

002B 

F603 


JMP 

X EQ+ 1 

N0.» GET MORE 



n-* 

OUTPUT T0 

TTY 


0020 

08 00 

'T* 

0TT 

E-JT 



002D 

4030 


SEL 

7, 4 

RESET INTERFACE 

002E 

6D3B 


URA 

7^ 3 

WRITE 0W NOT BUSY 

002F 

49 3B 


SE\I 

7, 3 

D2NE? 

0030 

F601 


JMP 

$- 1 

NO 

0031 

F705 


RTi'I 

DTT 

YES 








^icCOilMAI'JD ERROR EKIT 


0032 

08 00 

ERR 

H\JT 



0033 

C6DF 


LAP 

: DF 

PRINT ARR0U 

0034 

FE03 


JST 

0TT 


0035 

FF34 


JST 

❖ CMN D 

RESTART COMMAND 



^5= CARRI AG E- RETURN^ LIME 

FEED 

0036 

08 00 

CRLF 

ET^T 



0037 

068 D 


LAP 

:8D 

OR 

0033 

FEOC 


J ST 

0TT 


0039 

068 A 


LAP 

:8A 

LF 

003A 

FEOE 


JST 

0TT 


003B 

F705 


RTM 

CPJLF 




❖OUTPUT TEXT 

FR0M BUFFER 

0030 

03 00 

8TL 

e:>it 



003D 

8A0F 


ADD 

CAI 

MAICE COMPARE INSTRUCTION 

003E 

9A06 


STA 

0T2 

a SAVE IT 

003F 

9A09 


STA 

0T3 


0040 

E7 04 


LDIC 

❖ 0TL 

GET TEXT POINTER 

0041 

DE05 


IMS 

0TL 

SET RETURIxI ADDRESS 

0042 

B400 

0T1 

LDA 

60 

GET WORD 

0043 

1 1D7 


RRA 

8 

PRINT FIRST BYTE 

0044 

FE13 


JST 

0TT 


0045 

0000 

0T2 

CAI 

0 

LAST ONE 

0046 

F70A 


RTM 

0TL 

YES^ RETURN 

0047 

1157 


RLA 

8 

PRINT SECOND BYTE 

0048 

FEIC 


JST 

0TT 


0049 

0000 

0T3 

CAI 

0 

LAST ONE? 

004A 

F70E 


RTN 

■ 0TL 

YES RETURIvJ 

004B 

0128 


IXR 


BUMP POINTER 



182 


0107 

0C4C 

F60A 


JMP 

0T1 . 

LOOP 

0108 

0C4D 

COOO 

CAI 

CAI 

0 ' 


0109 



* 




0110 



* OUTPUT FLOATING POINT 

NUMBER 

0-1 1 1 






- 

0112 

0C4E 

03 00 

0FPA 

BIT 



0113 

004F 

E701 


LD)C 

*0FPA 

GET POINTER 

01 U 

0050 

DE02 


IMS 

0FPA 

SET RETUH'J ADDRESS 

01 15 

0051 

EAOl 


STIC 

0PT 

SAVE POINTER 

0116 

0052 

FBOE 


JST 

*FAS 

CONVERT T0 ASCII 

01 17 

0053 

0000 

0PT 

DATA 

0 


01 IS 

0054 

0059 


DATA 

BUF 


01 19 

0055 

0110 


ZAR 


SET END flag 

0120 

0056 

FElA 


JST 

0TL 

PRINT NIE-IBER 

0121 

0057 

0059 


DATA 

BUF 


0122 

0058 

F70A 


RTN 

0FPA 


0123 

0059 

0000 

BUF 

RES 

Sj 0 


0124 

0061 


FAS 

REF 



0125 

■ 


* 




0126 



0DEC 

0UTPUT DECIMAL ( 

C+/~ DDDDD) 

0127 



** 0DEC CONVERTS THE BINARY VALUE IN THE 

0128 



*A REG 

AI^D PRINTS IT AS A SIGNED 5 

0129 



=1; DIGIT DECIMAL NUMBER 

0N THE TELETYPE. 

0130 



*■ AR AND OV 

ARE DESTROYED. 

0131 






*■ 

0132 



^ L DA 

VAL AR 

: = VALUE 


0133 



* SM'l 

MUST 

BE IN U0RD 

MODE 

0134 



* JST 

❖ 0 DEC 

CALL ROUTINE 

0135 



* 

RETURI>I 

XR UNCHANGED 

0136 

0062 

08 00 

0DEC 

E\JT 



0137 

0063 

EA19 


STX 

S 

SAVE XR 

0138 

0064 

C4AB 


LXP 

I.}. » 


0139 

0065 

303 2 


JAP 

S+3 


0140 

0066 

0310 


NAR 


MAKE VALUE + 

0141 

0067 

C202 


AXI 

2 

MAKE SIGN - 

0142 

0068 

9A15 


STA 

V 

SAVE VALUE 

0143 

0069 

0030 


TXA 



0144 

006A 

FE3E 


JST 

0TT 

PRINT SIGN 

0145 

006B 

B215 


LDA 

STRT 


0146 

0060 

9A13 


STA 

PTR 

INITIALIZE TSL PTR 

0147 

006D 

C705 


LAt-I 

5 


0148 

006E 

9A10 


STA 

T 

SET FOR 5 DIGITS 

0149 

006F 

B20E 

0 1 

LDA 

V 


0150 

0070 

C4AF 


LXP 

: AF 

Z ER0 T0 - 1 

0151 

0071 

9 30E 


SUB 

*PTR 


0152 

0072 

0128 


IX R 



01-53 

007 3 

30C2 


JAP 

$- 2 . 


0^54 

0074 

8 BOB 


ADD 

*PTR 


0155 

0075 

9A0S 


STA 

V 


0156 

0076 

0030 


TXA 



0157 

0077 

FE4B 


JST 

0TT 

PRINT DIGIT 

0158 

0078 

DA07 


IMS 

■ PTR 


0159 

0079 

DAOS 


IMS 

T 




0160 

007A 

F60B 


JMP 

0 1 


0161 

007B 

E201 


LDK 

S 

RESTORE XR 

0162 

0070 

F7 lA 


RTW 

0DEC 

RETUHnI 

0163 

007D 

OSOO 

S 

HLT 


TB4P FOR XR 

016A 

007E 

08 00 

V 

HLT 


VALUE 

0165 

007F 

08 00 

T 

HLT 


COUTnIT 

0166 

003 0 

03 0 0 

PTR 

HLT 


POINTER 

0167 

003 1 

003 2 

STRT 

DATA 

TBL 

TABLE ADDR 

0168 

003 2 

27 10 

T3L 

DATA 

10000, 

1000, 100 

■ 

008 3 

03ES 






008 4 

0064 





0169 

008 5 

OOOA 


DATA 

10, 1 



008 6 

0001 





0170 




EMD 



0000 

ERRORS 
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0001 




NAl-I 

READ, PUMCri, CLEAR 

0002 




EKTR 

I KB, 0TT, 

BI PT, 0 TL, CRL F, I ER, ERR 

0303 

0000 



REL 

0 


0004 

0000 

03 00 

READ 

Ei\]T 



0005 

0 00-1 

F9-00 


J SJ 

CRLF 




0000 





0006 

0002 

F9 00 


J ST 

I ER 




0000 





0007 

0003 

0030 


CAI 

'O’ 


0008 

0004 

F203 


Ji-lP 

REZE 


0009 

0005 

C0C7 


CAI 

•G’ 


0010 

0006 

F204 


JMP 

RESM 


001 1 

0007 

f'9 00 


J ST 

ERR 




0000 





0012 

0008 

E21F 

REZE 

LEK 

IBC 


0013 

0009 

32 IF 


LDA 

ZCTO 


0014 

OOOA 

F204 


JMP 

RELP 


0015 

OOOB 

32 1C 

RES.J 

LDA 

IBC 


0016 

0000 

3A1E 


ADD 

SNO 


0017 

OOOD 

0043 


TAX 


ORIGINAL PAGE IS 

0018 

OOOE 

3213 


LDA 

Z CTl 

OP POOR QUALITT 

0019 

OOOF 

9A1C 

RELP 

STA 

CMT 


0020 

0010 

C68A 


LAP 

; 8A 


0021 

001 1 

F9 00 


JST 

0TT 




0000 





0022 

0012 

F9 00 


JST 

BIPT 




0000 





0023 

0013 

2141 


JAZ 

$- 1 

SKIP LEADER 

0024 

0014 

F9 00 

R1 

JST 

0TT 




0000 





0025 

0015 

F9 00 


JST 

BIPT 

READ TAPE 



0000 





0026 

0016 

009 2 


CAI 

:92 

EvID 0F LABE. 

0027 

0017 

F201 


JMP 

R2 

YES 

0023 

0018 

F604 


JMP 

R1 

m 

0029 

0019 

F9 00 

R2 

JST 

BIPT 

READ TAPE 



0000 





0030 

OOIA 

COFF 


CAI 

; FF 

FILE MARK 

0031 

OOIB 

F201 


JMP 

R3 

YES 

0032 

0010 

F603 


JMP 

R2 

N0 

0033 

OOID 

F9 00 

R3 

JST 

BIPT 




0000 





0034 

OOIE 

1357 


LLA 

8 


0035 

OOIF 

9A0E 


STA 

Tmp 


0036 

0020 

F9 00 


JST 

BIPT 




0000 





0037 

0021 

A20C 


10R 

TEMP 


0033 

0022 

3 000 


ADD 

eo 


0039 

0023 

9000 


STA 

60 


0040 

0024 

0128 


IXR 



0041 

0025 

DAO 6 


IMS 

CMT 


0042 

0026 

F609 


JMP 

• R‘3 


0043 

0027 

F727 


RTN 

READ 




185 


0044 

0028 


IBC 

REF 


004 5 

0029 

FOlO 

ZCTO 

DATA 

-403 0 

0046 

002A 

FFOl 

ZCTl 

DATA 

-255 

0047 

002B 

CEFl 

SMO 

DATA 

33 25 

0043 

002C 

0000 

CMT 

DATA 

0 

0049 

002D 

0000 

PNT 

DATA 

0 

0050 

002E 

0000 

TE'IP 

DATA 

0 

0051 



* 



0052 



* PUMCrf 


0053 

002F 

03 00 

PUMCri 

ENT 


0054 

0030 

F9 00 
0000 


JST 

CRLF 

0055 

0031 

F9 0C 
0000 


JST 

lER 

005 6 

0032 

COBO 


CAI 

•O' 

0057 

0033 

F203 


JMP 

pnze 

0053 

0034 

C0C7 


CAI 

' G ’ 

0059 

0035 

F204 


JMP 

pnsn 

0060 

0036 

F9 00 
0000 


J$T 

ERR 

0061 

0037 

E60F 

PMZE 

LEK 

IBC 

0062 

0038 

B60F 


LDA 

ZCTO 

0063 

0039 

F204 


JMP 

PNL P 

0064 

003A 

B612 

FM SW 

LDA 

IBC 

0065 

003B 

8E10 


ADD 

SI\JO 

0066 

0030 

0043 


TAK 


0067 

003D 

B613 


LDA 

Z CT 1 

0068 

003E 

9E12 

PMLP 

STA 

CNT 

0069 

003F 

FAOD 


J ST 

LEAD 

007 0 

0040 

F9 00 
0000 

PI 

J ST 

IKB 

0071 

0041 

C09 2 


CAI 

;92 

0 07 2 

0042 

F201 


JMP 

P2 

0073 

0043 

F603 


JMP 

PI 

0074 

0044 

C6FF 

P2 

LAP 

: FF 

007 5 

0 04 5 

F9 00 
0000 


J ST 

8TT 

007 6 

0046 

B400 

P3 

L DA 

§0 

0077 

0047 

FAOD 


JST 

0TU 

0078 

0048 

0123 


IXR 


0079 

0 049 

DEID 


IMS 

CMT 

008 0 

004A 

F604 


JMP 

P3 

003 1 

004B 

FAOl 


JST 

LEAD 

008 2 

0040 

F71D 


RTW 

PUN CH 

008 3 



* PUN OH 

•5" 0F 

L EADER 

008 4 

004D 

08 00 

LEAD 

EivIT ' 


008 5 

004E 

0732 


LAJ4 

50 

003 6 

004F 

9E22 


STA 

PNT 

008 7 

0050 

01 10 


Z AR 


0038 

0051 

F9 00 
0000 

L2 

JST 

0TT 

008 9 

0052 

DE25 


IMS 

PNT 

009 0 

0053 

F602 


JMP 

L2 


ORIGBSTAU PAGE IS 
OP POOR QUALEY 


ECH0 LABEL 

CTRL/TAPE 

YES 

W0 

PUiMCH FILE MARK 


PU\ICH BITS 
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009 1 
009 2 

0054 

F707 

009 3 

0055 

03 00 

009/1 

0056 

1 1D7 

009-5 

005’7 

F9 00 
0000 

009 6 

0058 

1 157 

009 7 

0059 

F9 00 
0000 

0098 

0099 

005A 

F7 05 

OiOO 

005B 

08 00 

0101 

005C 

F9 00 
0000 

0102 

005D 

F9 00 
0000 

0103 

005E 

COBO 

0104 

005F 

F203 

0105 

0060 

CCC7 

0106 

0061 

F204 

0107 

0062 

F9 00 
0000 

OIOS 

0063 

E63B 

0109 

0064 

B63B 

01 10 

0065 

F204 

01 1 1 

0066 

B63E 

01 12 

0067 

8E3C 

0113 

0068 

0048 

01 14 

0069 

B63F 

0115 

006A 

9E3E 

0116 

0C6B 

0110 

0117 

006C 

9000 

01 18 

006D 

0128 

01 19 

006E 

DE42 

0120 

006F 

F603 

0121 

0122 

007 0 

F715 

0000 

ERRORS 



5k 

RTW 

LEAD 

0TW 

EMT 



RRA 

3 


J ST 

0TT 


RLA 

8 


JST 

0TT 


RTW 

0TW 

* CL EAR 

PROGRAM 

CLEAR 

HaIT 



JST 

CRLF 


JST 

I ER 


CAI 

'0* 


JMP 

CLZ E 


CAI 

' G’ 


JMP 

CLSN 


JST 

ERR 

CLZE 

LDK 

IBC 


LDA 

Z CTO 


JMP 

CLLP 

CLSN 

LDA 

IBC 


ADD 

TAX 

SNO 


LDA 

Z CTl 

CLLP 

STA 

EAR 

CiMT 

CR 

STA 

IXR 

eo 


IMS 

CNT 


JMP 

CR 


RTW 

HMD 

CLEAR 
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0001 

0002 

0003 

0004 

0005 

0006 
0007 
0003 

0009 

0010 
0011 
0012 
0013 
00 U 

0015 

0016 
0017 
00 IS 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 
0027 
0023 

0029 

0030 
00 31 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 
00.48 

0049 

0050 

0051 

0052 

0053 





NAi'i 

FFL0T 


0000 



REL 

0 


0000 

03 00 

FPL0T 

ETJT 



0001 

C70F 


LAM 

15 


0002 

9 ABB 


STA 

NR0W 


0003 

B2D0 


LDA 

IBC 


0004 

9 ABA 


STA 

IB 


0005 

0108 


ZXR 



0006 

FBCF 


J ST 

*CRLF 


0007 

FBDl 


JST 

* ! ER 


0008 

COBO 


OAl 

•O' 


0009 

F221 


JMP 

FZR0 


GOOA 

GOBI 


CAI 

’ 1 ' 


OOOB 

F222 


JMP 

ForJE 


OOOC 

C0B2 


OAI 

*2* 


OOOD 

7222 


JMP 

FT1?0 


OOOE 

COBS 


OAI 

’3’ 


OOOF 

F222 


JMP 

FTHRE 


0010 

C0B4 


OAI 

’4' 


0011 

F222 


JMP 

FF0UR 


0012 

C0B5 


OAI 

'5' 


0013 

F222 


JMP 

fffve 


0014 

C0B6 


CAI 

’6’ 


0015 

F222 


JMP 

Fsrx 


0016 

C0B7 


CAI 

’7’ 


0017 

F222 


JMP 

fseve 


0018 

COBS 


CAI 

'8 * 


0019 

F222 


JMP 

FEI GH 


OOIA 

C0B9 


CAI 

• 9 . 


OOIB 

F222 


JMP 

FM rw E 


OOlC 

COCl 


CAI 

’A' 


OOID 

F222 


JMP 

FT£M 


OOIE 

C0C2 


CAI 

’B’ 


OOIF 

F222 


JMP 

FECE 


0020 

C0C3 


CAI 

'O' 


0021 

F222 


JMP 

FTWVE 


0022 

0004 


CAI 

*D’ 


0023 

F222 


JMP 

FrHD 


0024 

C0C5 


CAI 

' E ' 


0025 

F222 


JMP 

FFRTN 


0026 

0006 


CAI 

’F’ 


0027 

F222 


JMP 

ff^tn 


0028 

0007 


CAI 

’G' 


0029 

F222 


JMP 

FSXTN 


002A 

FBAF 


JST 

s):ERR 


002B 

0110 

FZR0 

ZAR 



002C 

9A9 0 


STA 

OPEL 


002D 

F22A 


JMP 

F00P1 


002E 

0601 

F0NE 

LAP 

1 

ORIGINAL PAGE IS 

002F 

0030 

F21E 

0602 

FTW0 

JMP 

LAP 

FST0 

2 

OF POOR QUALITY] 

0031 

F21C 


JMP 

FST0 


0032 

0603 

FTHRE 

LAP 

3 
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0054 

0033 

F21A 


JMP 

FST0 


0055 

0034 

0604 

FF0UR 

LAP 

4 


0056 

0035 

F218 


JMP 

FST0 


0057 

0036 

0605 

FFl VE 

LAP 

5 


0053 

00'3‘7 

F216 


JMP 

FST0 


0059 

0038 

0606 

FSIX 

LAP 

6 


0060 

0039 

F214 


JMP 

FST0 


0061 

003A 

0607 

FSEVE 

LAP 

7 


0062 

003B 

F212 


JMP 

FST0 


0063 

003C 

0603 

FEI GH 

LAP 

8 


0064 

003D 

F210 


JMP 

FST0 


0065 

003E 

0609 

FNINE 

LAP 

9 


0066 

003F 

F20E 


JMP 

FST0 


0067 

0040 

060A 

FTH^ 

LAP 

10 


0063 

0041 

F200 


JMP 

FST0 

okigmaepag:^ 
OF POOR OUALm 

0069 

0042 

060B 

FEu E. 

LAP 

1 1 

0070 

0043 

F20A 


JMP 

FST0 

0071 

0044 

0600 

FTWVE 

LAP 

12 

0072 

0045 

F208 


JMP 

FST0 


007 3 

0046 

060D 

fthd 

LAP 

13 


0074 

0047 

F206 


JMP 

FST0 


007 5 

0043 

060E 

FFRTO 

LAP 

14 


0076 

0049 

F204 


JMP 

FST0 


0077 

004A 

060F 

FFVTN 

LAP 

15 


0078 

004B 

F202 


JMP 

FST0 


0079 

004C 

0610 

FSXTO 

LAP 

16 


003 0 

004D 

F200 


JMP 

FST0 


003 1 

004E 

9A7A 

FST0 

STA 

OALL 


008 2 

004F 

0310 


MAR 



008 3 

0050 

9A79 


STA 

FR0¥ 


003 4 

0051 

326D 


LDA 

IE 


008 5 

0052 

DA77 

FR0N0 

IMS 

FR0W 


003 6 

0053 

F201 


JMP 

S+ 2 


003 7 

0054 

F202 


JMP 

S+ 3 


0088 

0055 

8A75 


ADD 

H255 


0089 

0056 

F604 


JMP 

FR0N0 


009 0 

0057 

9A67 


STA 

IB 


009 1 

0058 

B270 

F00P1 

LDA 

OALL 


009 2 

0059 

2101 


JAZ 

•S+2 


009 3 

005A 

F204 


JMP 

S+ 5 


009 4 

005B 

B260 


LDA 

ZCT 


009 5 

005C 

0150 


lAR 



009 6 

005D 

9A63 


STA 

CMT 


009 7 

005E 

F203 


JMP 

$+ 4 


0098 

005F 

07FF 


LAI<1 

255 


0099 

0060 

0150 


lAR 



-0100 

0061 

9 ASF 

- 

STA 

ONT 


0101 

0062 

B25C 


LDA 

IB 


01*02 

0063 

9A5E 


STA 

IPBl 


0103 

0064 

B35D 


LDA 

*IPB1 


0104 

0065 

9A5E 


STA 

MAX 


0105 

0066 

DA5B 

L00P2 

IMS 

IPBl 


0106 

0067 

B35A 


LDA 

>IPB1 




01 

07 

0068 

ol 

08 

0069 

01 

09 

006A 

01 

10 

006B 

01 

1 1 

0060 

01 

12 

006D 

01 

13 

006E 

01 

14 

006F 

01 

15 

0070 

01 

16 

0071 

01 

17 

0072 

01 

IS 

0073 

01 

19 

0074 

01 

20 

0075 

01 

21 

0076 

01 

22 

0077 

01 

23 

0078 

01 

24 

0079 

01 

25 

007A 

01 

26 

007B 

01 

27 

0070 

0128 

007D 

0129 

007E 

01 

30 

007F 

01 

31 

008 0 

01 

32 

OGS 1 

01 

33 

003 2 

01 

34 

003 3 

01 

35 

008 4 

01 

36 

008 5 

01 

37 

008 6 

01 

38 

003 7 

01 

39 

0038 

-0140 

0089 

0141 

003 A 

01 

42 

003 B 

01 

43 

008 0 

0144 

003 D 

01 

45 

008 E 

01 

46 

008 F 

01 

47 

009 0 

01 

48 

009 1 

01 

49 

009 2 

01 

50 

009 3 

01 

51 

0094 

01 

52 

009 5 

01 

53 

009 6 

01 

54 

009 7 

01 

55 

0098 

01 

56 

0099 

01 

57 

009A 

01 

58 

009B 

01 

59 

009 C 


D25B 
F201 
9A59 
DA55 
F606 
0350 
9A56 
9A56 
B256 
D252 
F202 
F20C 
F20B 
B24E 
0103 
9 24F 
0128 
31C2 
EA4A 
1200 
1 iPB 
3201 
0128 
EA46 

FB55 L08P 

FB54 

C6ED 

FB54 

OOCD 

C6A0 

FB50 

B23D 

FB4C 

C601 

9A41 

FB4A 

F349 

B23E 

FB4S 

C7FF 

9A2F 

B22D 

9A30 

FB42 L00P1 

B32E 

FB3F 

C6A0 

FB3F 

C6B0 

FB3D 

B323 

2192 

0103 


CMS 

MAX 

JMP 

5+2 

STA 

MAX 

IMS 

CNT 

JdP 

L00P2 

ARP 


STA 

SCALE 

STA 

R0UND 

LDA 

H50 

CMS 

MAX 

JMP 

5+ 3 

JMP 

L00P 

JMP 

L00P 

LDA 

MAX 

ZXR 


SUB 

H50 

IXR 


JAG 

5-2 

STX 

SCALE 

R0V 


RRX 

1 

J0R 

$+2 

IXR 


STX 

ROUI-JD 

JST 

^}=CRLF 

JST 

=}=CRLF 

LAP 

t ~ r 

JST 

*0TL 

DATA 

SF 

LAP 

f r 

JST 

=R0TT 

LDA 

SCALE 

J ST 

sfi0DEC 

LAP 

1 

STA 

R0N0 

JST 

*CRLF 

J ST 

CRLF 

LDA 

R0W0 

JST 

*0TT 

LA'-I 

255 

STA 

NC0L 

LDA 

IB 

STA 

IPB2 

JST 

=!«CRLF 

LDA 

*IPB2 

JST 

^0DEC 

LAP 

f t 

JST 

>R0TT 

LAP 

• 0’ 

JST 

^0TT 

LDA 

* I PB2 

AL» 

CLOSE 

ZXR 
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0160 

009 D 

9 227 


SUB 

SCALE 

0161 

009 E 

0128 


IXR 


0162 

009 F 

3102 


JAG 

S-2 

0163 

OOAO 

8A24 


ADD 

SCALE 

0164 

OO'Al 

00A3 


n?CR 


0165 

00A2 

D223 


CMS 

R0U1'JD 

0166 

00 A3 

F202 


JMP 

S+ 3 

0167 

00A4 

0000 


NOP 


016B 

00A5 

0128 


IXR 


0169 

00A6 

0030 


TXA 


0170 

00A7 

218 6 


JAL 

CLOSE 

0171 

OOAS 

C6AA 


LAP 

• 

0172 

00A9 

0503 


MXR 


0173 

OOAA 

EA16 


STX 

CNT 

0174 

OOAB 

FB2B 


JST 

5}:0TT 

0175 

OOAC 

DAI 4 


IMS 

CNT 

0176 

00 AD 

F602 


JMP 

S- 2 

0177 

OOAE 

DA14 

CLOSE 

IMS 

IPB2 

0178 

OOAF 

0000 


MOP 


0179 

OOBO 

DAOF 


IMS 

WOOL 

OISO 

OOBl 

F61E 


JMP 

L00P1 

018 1 

0 0B2 

FB23 


JST 

>r=CRLF 

018 2 

00B3 

B215 


LDA 

CALL 

0133 

00B4 

2101 


JAZ 

S+ 2 

0184 

00B5 

F206 


JMP 

FSH 

0185 

OOB6 

C6FF 


LAP 

255 

0186 

00B7 

3A07 


ADD 

IB 

018 7 

00B8 

9A06 


sta 

IB 

0188 

00B9 

DAI 2 


IMS 

R8W0 

0139 

OOBA 

DA03 


IMS 

WR0W 

0190 

OOBB 

F623 


JMP 

L00P1 

019 1 

OOBC 

FE19 

FSH 

JST 

*CRLF 

0192 

OOBD 

F7BD 


RTN 

FPL0T 

019 3 

OOBE 

0000 

NR0W 

DATA 

0 

0194 

OOBF 

0000 

IB 

DATA 

0 

0195 

OOCO 

0000 

WOOL 

DATA 

0 

0196 

0001 

0000 

CMT 

DATA 

0 

0197 

0002 

0000 

IPBl 

DATA 

0 

0198 

0003 

0000 

IPB2 

data 

0 

0199 

0004 

0000 

MAX 

DATA 

0 

0200 

0005 

0000 

SCALE 

DATA 

0 

0201 

0006 

0000 

ROUND 

DATA 

0 

0202 

0007 

0032 

H50 

DATA 

50 

0203 

0003 

FI OF 

Z CT 

DATA 

-33 25 

02 04 

0009 

0000 

CALL 

DATA 

0 

0205 

OOCA 

0000 

FR0U 

DATA 

0 

0206 

OOCB 

OOFF 

H255 

DATA 

255 

0^07 

OOCO 

0000 

R0N0 

DATA 

0 

0203 

OOCD 

D3C3 

SF 

TDCT 

’SCALE FACTO R= 


OOCE 

ClCC 





OOCF 

C5A0 





OODO 

0601 





OODl 

C3D4 






0209 

00D2 

00D3 

00D4 

CFD2 

BDAO 

IBC 

REF 

0210 

00D5 


0DEC 

REF 

021 1 

OOD6 


CRLF 

REF 

0212 

00D7 


0TT 

REF 

0213 

00D3 


0TL 

REF 

0214 

00D9 


lER 

REF 

0215 

OODA 


ERR 

REF 

0216 

0000 

ERRORS 



HSID 



0001 




NAI«1 

CRT 


0002 




NA’l 

DISP 


0003 




NAM 

DISD 


0004 




BCTR 

RTDA 


0005 




EXTR 

RTTR 


0006 

0000 



REL 

0 


0007 

0000 

03 00 

CRT 

EMT 



0008 

0001 

C70F 


la:-i 

15 


0009 

0002 

9900 
01 12 


STA 

NP^U 


0010 

0003 

49 C 4 


SE'J 

24, 4 


0011 

0004 

F601 


JMP 

£- 1 


0012 

0005 

4007 


SEL 

24, 7 


0013 

0006 

0110 


ZAR 



0014 

0007 

9900 

0113 


STA 

OFFSET 


0015 

0003 

0601 


LAP 

1 


00.16 

0009 

9900 
01 16 

DI SP 

STA 

TRi 


0017 

OOOA 

3106 


JAI\! 

S+7 


0018 

OOOB 

BlOO 


LDA 

ADI 6 


0019 

OOOC 

0126 

8900 

0133 


ADD 

IBC 

0020 

GOOD 

9900 

0114 


STA 

IB 


0021 

OOOE 

0110 


ZAR 



0022 

OOOF 

9900 

0117 


STA 

CALL 


0023 

0010 

F25F 


JMP 

MXMI 


0024 

001 1 

BiOO 


LDA 

IBC 



' 

0133 





0025 

0012 

9900 

0114 


STA 

IB 


0026 

0013 

0103 


ZXR 



0027 

0014 

F9 00 
8 133 


JST 

M^CRLF 


0023 

0015 

F9 00 
8 134 


JST 

^-lER 


0029 

0016 

COBO 


CAI 

• 0’ 


0030 

0017 

^F221 


JMP 

Zero 


0031 

0018 

GOBI 


CAI 

• 1 • 


0032 

0019 

F222 


JMP 

0ME 


0033 

OOIA 

G0B2 


CAI 

*2’ 


0034 

OOIB 

F222 


JMP 

TWO 


0035 

OOlC 

0033 


CAI 

’3’ 


0036 

OOID 

F222 


JMP 

THREE 


0037 

OOlE 

C0B4 


CAI . 

•4' 


00 S8 

OOIF 

F222 


JMP 

F0Dr 


0039 

0020 

C0B5 


CAI 

’5* 


0040 

0021 

F222 


JMP 

FI WE 


0041 

0022 

C0B6 


CAI 

•’6' ‘ 


0042 

0023 

F222 


JMP 

SIX 
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0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

007 2 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

008 0 
003 1 
003 2 
003 3 
003 4 
003 5 
003 6 

008 7 
0038 
0089 

009 0 
009 1 
009 2 
009 3 
009 4 
009 5 


0024 

C0B7 


CAI 

. 7 . 


0025 

F222 


JIIP 

SE\/EM 

obigwalpageK 
OE POOR QPALIi’‘ 

0026 

0027 

0028 

COBS 

F222 

0039 


CAI 

JMP 

CAI 

’8 ' 

El GHT 

.9 • 

0029 

F222 


JMP 

NINE 


002 A 

COCi 


CAI 

’A’ 


002B 

F222 


■jMP 

TEN 


002C 

0002 


CAI 

’B' 


002D 

F222 


JMP 

ELE 


002E 

0003 


CAI 

’ 0 * 


002 F 

F222 


JMP 

tuve 


0030 

0004 


CAI 

’D’ 


0031 

F222 


JMP 

THD 


0032 

C0C5 


CAI 

•E’ 


0033 

F222 


JMP 

frtn 


0034 

0006 


CAI 

•F* 


0035 

F222 


JMP 

FVTN 


0036 

0007 


CAI 

»G’ 


0037 

F222 


JMP 

SXTIn! 


0038 

FBFC 


JST 

=t=ERR 


0039 

0110 

Z ER0 

EAR 



003A 

9 ABC 


STA 

CALL 


003B 

F22E 


JMP 

L00P1 


003C 

0601 

ONE 

LAP 

1 


003D 

F222 


JMP 

STORE 


003E 

0602 

TU 0 

LAP 

2 


003F 

F220 


JMP 

STORE 


0040 

0603 

THREE 

LAP 

3 


0041 

F21E 


JM P 

STORE 


0042 

0604 

FOUR 

LAP 

4 


0043 

F21C 


JMP 

STORE 


0044 

0605 

FI VE 

LAP 

5 


0045 

F21A 


JMP 

STORE 


0046 

0606 

SIX 

LAP 

6 


0047 

F218 


JMP 

STORE 


0048 

0607 

SEVEl 

LAP 

7 


0049 

F216 


JMP 

STORE 


004A 

0603 

EIGHT 

LAP 

8 


0043 

F214 


JMP 

STORE 


0040 

0609 

NINE 

LAP 

9 


004D 

F212 


JMP 

STORE 


0C4E 

C60A 

TEN 

LAP 

10 


004F 

F210 


JMP 

STORE 


0050 

C60B 

ELE 

LAP 

1 1 


0051 

F20E 


JMP 

STORE 


0052 

0600 

TUVE 

LAP 

12 ■ 


0053 

F20C 


JMP 

STORE 


0054 

C60D 

THD 

LAP 

13 


0055 

F20A 


JMP 

STORE 


0056 

C60E 

FRTN 

LAP 

14 


0057 

F208 


JMP 

- STORE 


0058 

C60F 

FVTN 

LAP 

15 




194 


009 6 
009 7 

0098 

0099 
01-00 
0101 
0102 

0103 

0104 

0105 

0106 
0107 
0103 
0109 
01 10 
01 1 1 
0112 
0113 
01 14 

0115 

0116 
01 17 
01 18 
01 19 
0120 
0121 
0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 
of 44 

0145 

0146 

0147 

0148 


0059 

F206 


JMP 

ST0RE 

005A 

0610 

SXTN 

LAP 

16 

005B 

F204 


JMP 

ST0RE 

005C 

0108 

D1 SD 

ZXR 


005-D 

EAB8 


STX 

TRl 

005E 

E2D4 


LIK 

I BO 

005F 

EAB4 


STX 

IB 

0060 

9AB6 

STORE 

STA 

CALL 

0061 

0310 


MAR 


0062 

9AB5 


STA 

R0 UVA 

0063 

B2B0 


LDA 

IB 

0064 

DABS 

R0VJN0 

IMS 

R0UVA 

0065 

F201 


JMP 

S+2 

0066 

F202 


JMP 

$-}■ 3 

0067 

8 ABB 


ADD 

H255 

0068 

F604 


JMP 

R0V.JN0 

0069 

9 AAA 


STA 

13 

006A 

B2A0 

L00P1 

LDA 

CALL 

006B 

2101 


JAZ 

S+2 

006C 

F203 


JMP 

$+ 4 

006D 

B2A0 


LDA 

ZCT 

006E 

9AA0 


STA 

CNT 

006F 

F202 


JMP 

$+3 

0070 

07FF 

MXMI 

LAfl 

255 

0071 

9AA9 


STA 

CNT 

0072 

B2A1 


LDA 

IB 

0073 

9AAS 


STA 

IBPl 

0074 

0110 


Z AR 


0075 

9AA3 


STA 

MIN 

0076 

9-AA8 


STA 

MAX 

0077 

E3A4 

L00P2 

LDA 

^i^IBPl 

0078 

D2A6 


CMS 

MAX ' 

0079 

F202 


JMP 

$+ 3 

007A 

9AA4 


STA 

MAX 

007B 

F203 


JMP 

£+4 

007C 

D2A1 


CMS 

MIN 

007D 

9AA0 


STA 

MIN 

007E 

0000 


NOP 

■ 

007F 

DA9C 


IMS 

IBPl 

008 0 

DA9A 


IMS 

CNT 

003 1 

F60A 


JMP 

L00P2 

008 2 

0350 


ARP 


008 3 

9A9C 


STA 

SCALE 

008 4 

9A9C 


STA 

R0 UN D 

008 5 

3299 


LDA 

MAX 

008 6 

9 29 7 


SUB 

MIN 

008 7- 

9-A9--7 

- 

STA 

'MAX 

0088 

B28D 


LDA 

TRl 

0089 

2106 


JAZ 

SCA 

008A 

B28C 


LDA 

CALL 

008 B 

2101 


JAZ 

$+2 

008 0 

F203 


JMP 

■ $+4 

008 D 

B29 4 


LDA 

HI 00 
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0U9 

008 E 

9A98 


STA 

T0TD0T 


0150 

008 F 

F202 


JMP 

S+3 


0151 

009 0 

B29 3 

SCA 

LDA 

H200 


0152 

009 1 

9A9 5 


STA 

T0TD0T 


0153 

009 2 

B29 4 


LDA 

T0TD0T 


0154 

009 3 

D23B 


aMS 

MAX 


0155 

009 4 

F202 


JMP 

S+ 3 


0156 

009 5 

F20C 


JMP 

L00P 


0157 

009 6 

F20B 


JMP 

L00P 


0153 

009 7 

B23 7 


LDA 

MAX 


0159 

0160 

0098 

0099 

0108 
9 28D 


ZXR 

SUB 

T0TD0T 

ORIGINAI] PAGE IS 

0161 

0 09 A 

0128 


IXR 


OE POOR QTJALTTy 

0162 

009 B 

31C2 


JAG 

S-2 


0163 

009 C 

EA8 3 


STX 

SCALE 


0164 

009 D 

1200 


R0V 



0165 

009 E 

1 1A3 


RPX 

1 


0166 

009 F 

3201 


J0R 

S+2 


0167 

00 AO 

0128 


IXR 



0163 

OOAi 

EA7F 


STX 

ROUND 


0169 

00A2 

B273 

L00P 

L DA 

TRl 


0170 

OOA3 

21 OB 


JAZ 

MI 


0171 

OOA4 

FB9 3 


JST 

CRLF 


0172 

00A5 

FB9 2 


JST 

:{:CRLF 


0173 

C0A6 

C6BD 


LAP 

1 — t 


0174 

00A7 

FB5F 


JST 

^0TL 


0175 

00A8 

012C 


DATA 

SF 


0176 

00A9 

C6A0 


LAP 

f t 


0177 

OOAA 

FB8B 


JST 

:P0TT 


0173 

OOAB 

B274 


LDA 

SCALE 


0179 

OOAC 

FBSC 


JST 

5ti0 DEC 


0180 

00 AD 

FBSA 


JST 

*CRLF 


013 1 

OOAE 

FB39 


JST 

*CRLF 


0182 

OOAF 

B26E 

MI 

LDA 

MIN 


018 3 

OOBO 

0310 


WAR 



013 4 

0031 

0108 


ZXR 



0185 

00B2 

926D 


SUB 

SCALE 


0186 

00B3 

0128 


IXR 



0187 

00B4 

3102 


JAG 

S-2 


0188 

0035 

8A6A 


ADD 

SCALE 


0139 

00B6 

00 AS 


DXR 

■ 


019 0 

0037 

D269 


CMS 

ROUND 


019 1 

00B3 

F202 


JMP 

S+3 


0192 

00B9 

0000 


NOP 



0193 

OOBA 

0128 


IXR 



0194 

OOBB 

EA6C 


STX 

XI 


0195 

OOBC 

B26B 

L06P4 

LDA 

X 1 


019 6 

OOBD 

9A6B 


STA 



01*9 7 

OOBE 

C7FF 


LA-I 

255 


0198 

003F 

9A59 


STA 

NC0L 


0199 

OOCO 

B253 


LDA 

IB 


0200 

OOCl 

9A5B 


STA 

■ IBP2 


0201 

00C2 

B250 

L00P3 

LDA 

OFFSET 




0202 

00C3 

SA65 


add 

X 

0203 

00C4 

9A66 


STA 

TO TOFF 

0204 

0005 

B357 


EDA 

^=IBP2 

0205 

00C6 

1200 


R0V 


0206 

00C7 

1 150 


RLA 

1 

0207 

00 OS 

2207 


JOS 

NEG 

0203 

00C9 

1 I'DO 


RRA 

1 

0209 

OOCA 

0108 


ZXR 


0210 

OOCB 

9 254 


SUB 

SOALE 

021 1 

OOCC 

0128 


IXR 


0212 

OOCD 

31C2 


JAG 

$-2 

0213 

OOCE 

8A51 


ADD 

SCALE 

0214 

OOCF 

F206 


JllP 

S+ 7 

0215 

OODO 

1400 

NEG 

S0V 


0216 

OODl 

1 IDO 


RRA 

1 

0217 

00D2 

0108 


ZKR 

Sp^LE 

0218 

OOD3 

8A4C 


ADD 

0219 

00D4 

00-^ 


EKR 


0220 

00D5 

2002 


JA'I 

S- 2 

0221 

OOD6 

9 249 


SUB 

SCi^E 

0222 

00D7 

00 AS 


DXR 


0223 

0028 

D248 


OMS 

ROUND 

0224 

OOD9 

F202 


JMP 

$+ 3 

0225 

OODA 

0000 


MOP 


0226 

OODB 

0128 


IXR 


0227 

OODC 

0030 


TXA 


0228 

OODD 

9A4C 


STA 

DO TWO 

0229 

CODE 

E24C 


LD>C 

TO TOFF 

0230 

OODF 

49 04 

DOT 

SE'l 

24^ 4 

0231 

OOEO 

F6Q1 


JMP 

S- 1 

02'32 

OOEl' 

49 02 


SEW 

24. 2 

0233 

00E2 

F601 


JMP 

S“ 1 

0234 

00E3 

6EC2 


OTJC 

24; 2 

0235 

00E4 

B231 


LDA 

TRl 

0236 

00E5 

2108 


JAZ 

S+9 

0237 

00E6 

B230 


LDA 

CALL 

0233 

00E7 

2101 


JAZ 

S+2 

0239 

00E8 

F205 


JWP 

£+ 6 

0240 

00E9 

B240 


LDA 

DO TWO 

0241 

OOEA 

8A40 


ADD 

T0T0FF 

0242 

OOEB 

0043 


TAX 


024 3 

OOEC 

6EC2 


0TX 

24> 2 

0244 

OOED 

F20B 


JMP 

DOME 

0245 

OOEE 

B23B 


LDA 

D0TN0 

0246 

OOEF 

2109 


JAZ 

DOME 

0247 

OOFO 

318 4 


JAG 

P0S 

0243 

OOFl 

00 AS 


DXR 


0249 - 

00F2 

0T50 


lAR 


0250 

OOF3 

9A36 


STA 

D0TN0 

0251 

00F4 

F615 


JMP 

DOT 

0252 

0 0F5 

0123 

P0S 

IXR 


0253 

00F6 

oOdo 


DAR 

. 

0254 

00F7 

9 A3 2 


STA 

DOTInIO 



0255 

OOFS 

F619 


JMP 

D0T 

0256 

00F9 

B22F 

D0NE 

LDA 

X 

0257 

OOFA 

8 A2A 


ADD 

H256 

0253 

OOFB 

9A2D 


STA 

yv 

0259 

OOFC 

DA2 0 


IMS 

IBP2 

0260 

OOFD 

DAIB 


IMS 

NC0L 

0261 

OOFE 

F63C 


JMP 

L00P3 

0262 

OOFF 

B216 


LDA 

TRl 

0263 

0100 

3104 


JAN 

$+5 

0264 

0101 

B215 


LDA 

CALL 

0265 

0102 

2101 


JAZ 

S+2 

0266 

0103 

FlOO 

0000 


JMP 

RTDA 

0267 

0104 

FlOO 

OOOO 


JMP 

RTTR 

0263 

0105 

B21 1 


LDA 

CALL 

0269 

0106 

2101 


JAiZ 

$+2 

027 0 

0107 

F203 


JMP 

FIN 

0271 

0103 

C60A 


LAP 

10 

0272 

0109 

8A09 


ADD 

OFFSET 

027 3 

OlOA 

9 AOS 


STA 

OFFSET 

0274 

01 OB 

B217 


LDA 

H255 

0275 

OlOC 

SA07 


ADD 

IB 

0276 

OlOD 

9A06 


STA 

IB 

0277 

OlOE 

DAO 3 


IMS 

NR0U 

0278 

01 OF 

F653 


JMP 

L00P4 

0279 

0110 

FB27 

FIN 

J ST 

>1=CRLF 

028 0 

01 1 1 

FlOO 

8000 


RTN 

CRT 

023 I 

0112 

0000 

MRSV 

DATA 

0 

023 2 

0113 

0000 

0FFSET 

DATA 

0 

023 3 

0114 

0000 

IB 

DATA 

0 

028 4 

01 15 

0000 

I3Q 

DATA 

0 

023 5 

0116 

0000 

TRl 

DATA 

0 

023 6 

0117 

0000 

GAEL 

DATA 

0 

023 7 

01 18 

0000 

R0WVA 

DATA 

0 

0238 

01 19 

0000 

NC0L 

DATA 

0 

0239 

01 lA 

FlOE 

Z CT 

DATA 

“38 26 

029 0 

01 IB 

0000 

CNT 

DATA 

0 

029 1 

01 1C 

0000 

IBPl 

DATA 

0 

029 2 

01 ID 

0000 

IBP2 

DATA 

0 

029 3 

01 IE 

0000 

MIN 

DATA 

0 

029 4 

01 IF 

0000 

MAX 

DATA 

0 

029 5 

0120 

0000 

SCALE 

DATA 

0 

029 6 

0121 

0000 

ROUiy'D 

DATA 

0 

029 7 

0122 

0064 

HI 00 

DATA 

100 

0298 

0123 

OOFF 

H255 

DATA 

255 

0299 

0124 

oocs 

H200 

DATA 

200 ' 

0300 

0125 

0100 

H256 

DATA 

25 6" 

0301 

0126 

OEFl 

AD16 

DATA 

38 25 

0302 

0127 

0000 

T0TD0T 

DATA 

0 

0303 

0128 

0000 

X 1 

DATA 

0 

0304 

0129 

0000 

X ■ 

DATA 

0 



0305 

012A 

0000 

D0TN0 

DATA 

0306 

0123 

0000 

T0T0 FF 

DATA 

0307 

012C 

D3C3 

SF 

TEKT 


0I2D 

ClCC 




0-L2-E 

C-SAG 


- 


012F 

C6C1 




0130 

C3D4 




0131 

CFD2 




0132 

BDAO 



0308 

0133 


IBC 

REF 

0309 

0134 


lER 

REF 

0310 

0135 


ERR 

REF 

0311 

0136 


0TT 

REF 

0312 

0137 


0TL 

REF 

0313 

0133 


CRLF 

REF 

0314 

0139 


0DEC 

REF 

0315 




E\]D 

0000 

ERRORS 





0 

0 

’SCALE FACT0R= 



PAGE 

0001 


0001 



0002 

0000 


0003 

0000 

0001 

OOOA 

0001 

0001 

0005 

0002 

0001 

0006 

0003 

0001 

0007 

0004 

FFFE 

0008 

0005 

0001 

0009 

0006 

FFFF 

0010 

0007 

0001 

0011 

0003 

0001 

0012 

0009 

FFFF 

0013 

OOOA 

FFFF 

0014 

OOOB 

0001 

0015 

oooc 

FFFF 

0016 

GOOD 

FFFF 

0017 

OOOE 

FFFF 

0018 

OOOF 

0001 

0019 

0010 

0001 

0020 

0011 

000 1 

0021 

0012 

0001 

0022 

0013 

FFFF 

0023 

0014 

0001 

0024 

0015 

FFFF 

0025 

0016 

0001 

0026 

0017 

0001 

0027 

00 IS 

FFFF 

0028 

0019 

FFFF 

0029 

001 A 

0001 

0030 

0013 

FFFF 

0031 

OOlG 

FFFF 

0032 



0000 

ERRORS 


NAl-i MASKO 

REL 0 

MASKO DATA 1 

DATA 1 

DATA 1 

DATA 1 

DATA - 1 

DATA I 
DATA - I 

DATA I 
DATA 1 
DATA - I 

DATA “■ 1 

DATA I 
DATA - 1 

data - 1 
DATA - 1 
DATA I 
DATA i 
DATA 1 
DATA 1 
DATA - 1 
DATA I 
DATA - 1 
DATA 1 ■ 
DATA 1 
DATA - 1 
DATA “ 1 
DATA I 
DATA “ I 
DATA - i 
Ei'JD 



200 


0001 


NAM 

IBC 

0002 

0000 

REL 

0 

0003 

0000 0000 IBC 

DATA 

0000 

0004 

0000 

ERRORS 

EHD 






APPENDIX D ORIGINAZ PAffP tq 


1 DEF FiMRCX) = l:\IT( lOO+X+0. 5)/ iOO 

10 DIM 0C 255)^ F( 255),XC 255) 

20 PRINT "SPECTRUM C0RRECTED F0R NEGETIVE DIP" 

25 PRINT 

26 PRINT 

30 CALL C20) 

40 CALL <5,0(0)^ 256^2) 

50 LET F<O) = 0(O) 

60 F0R W= 1 T0 255 

70 LET Il=N-i 

75 I F I I> 0 THD>I 3 5 

80 LET I 1=255 

85 LET 12=N+ 1 

90 I F I 2< 256 THEN IOC 

95 LET 12=12-255 

100 LET l3=N+24 

105 IF I3<256 THEs] 115 

no LET 13=13-255 
115 L ET I 4=N+ 25 

120 IF I4<256 THE'J 130 

125 LET 14=14-255 

130 LET C=(0( I 3) + 0( I4)-0CI 1)-0C 12) )/20 
140 LET F(N)=0(N) + C 
150 NEXT N 

160 PRINT "A; 1 F0R DISPLAY, 2 FOR GRAPH, 3 F0R TAPE" 

170 PRINT "B: 0 F0R ORIGINAL, 1 F0R FINAL" 

13 0 PRINT 

200 PRINT "A='U 

210 INPUT A 

220 PRINT "B=”; 

230 INPUT B 

240 PRINT "RESOLUTION D= 

250 INPUT D 

260 PRINT "INITIAL WAVE LENGTH LO="J 

270 INPUT LO 

23-0 PRINT 

29 0 PRINT 

300 IF A=3 THEM 730 

310 . LET Z=0 

320 LET M= lElO 

330 F0R N= 1 T0 255 


201 



202 


350 IF B= 1 TriEiJ 400 

360 LLT X(.']> = 0(N) 

370 G0T0 410 

400 L£T X<I4)=FCiM) 

410 IF XCN)> = i'l TUHnI 430 

420 LET M=XCiM) 

430 IFX<W)<=Z THEM 450 

440 LET Z=X(iM) 

450 NEXT N 

460 IF ri>=0 THEN 510 

470 LET SO=255/(2-N) 

^0 LET S1=43/CZ-N) 

49 0 LET Y0=-M*S 

500 G0T0 540 

510 LET S0=255/Z 

520 LET Sl=43/Z 

530 LET Y0=0 

540 PRINT "NAX=”;Z; ’MIN=’’;N 

550 IF A=2 THErsI 630 

560 PRINT "SCALE FACT0R=’*;SO 

570 F0R N= 1 T0 255 

53 0 LET E=INT( S0*XCN) + 0. 5) 

59 0 CALL C3^N^Y0>2, E) 

600 NEXT N 

620 G0T0 160 

630 PRINT "SCALE FACTO R="J SI 

635 PRINT 

636 PRINT 

640 FOR N= 1 T0 255 

650 LET E=INTC Sl>»:XCN) + 0. 5) 

660 LET L=L0+CN- 1)=(=D 

670 PRINT FNRCL);TABC8)JX(N) 

715 N EX T N 

730 IF B=0 THH'I 760 

740 IF B= 1 THEN 73 0 

750 G0T0 200 

760 CPI.L< 6^ 0C 0>^ 256^ 2) 

77 0 G0T0 79 0 

78 0 CALLC 6^ F< 0)^ 256^ 2) 

79 0 CALLC 6^ 0^ Oj 3) 

800 ST0P 


OfiajBSAj, PArp 1 
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